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SUMMARY 


Our experimental and theoretical effort continues to be aimed at the 
establishment of relationships among crystal growth parameters* materials 
properties, electronic properties and device applications of GaAs. Toward 
this goal» we undertook the development of new approaches to the preparation 
and characterization of GaAs. This extensive ground-based program constitutes 
a necessary step for insuring successful processing of GaAs under zero-gravity 
conditions. Furthermore, due to its unique scope combining crystal-growth 
characterization, and device-related properties and phenomena, this program 
bears directly on the exploitation of the potential of GaAs in device applica- 
tions. 

Our crystal growth effort includes electroepitaxy as a non-conventlonal 
and extremely promising approach to liquid phase epitajqr and melt growth of 
bulk GaAs crystal. We have substantially advanced the understanding of funda- 
mental processes involved in electroepltaxlal growth. Our quantitative model 
of electroepltaxlal growth in terms of electromlgratlon and the Peltier effect 
has been successfully extended to the growth of ternary and quarternary compounds 
predicting high growth rates and remarkable stabilization of the composition of 
electroepltaxlal layers. 

High growth rates attained by electroepitaxy confirmed out prediction 
that interface kinetic phenomena have only negligible effects on electroepl- 
taxlal growth. Thus, it is now experimentally verified that electroeplt^ucy 
enables the extension of the LPE process to high growth rates so far obtained 
only in melt growth. As clearly demonstrated by our growth-property studies, 
high growth rates result in a remarkable Improvement of the carrier recombina- 
tion characteristics of epitaxial layers. 

In our melt-growth studies we have completed the construction and testing 
of an advanced system for horizontal and/or vertical growth of GaAs. Successful 


f 
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growth of bulk GaAs under precisely controlled growth conditions has been 
Initiated. Crystals recently obtained were extensively characterised 
on a macro- and micro-scale. 

Our effort to establish a "state of the art" characterization facility 
has been 'highly successful. A modified system for Deep Level Transient 
Capacitance Spectroscopy was up, suitable for the determination of 

bulk levels and Interface states. Scanning Electron Mlcroscopy-cathodolumlnes- 
cence studies were Initiated In the temperature range of 80K to 700K. 

These newly Introduced techniques, together with our recently developed 
methods such as IR scanning absorption, derivative photovoltage and photocapacl- 
tance spectroscopies, transport techniques and SEM-electron beam-induced current 
now permit the unique characterization of all essential electronic parameters 
of GaAs on a macro- and a micro-scale. 

The Importance of micro-characterlzatlon In establishing growth-property 
relationships of GaAs was recently demonstrated through our studies of 
carrier concentration and compensation mlcropruflles in melt-grown material. 
These studies showed clearly that, unlike those of elemental semiconductors 
(SI, Ge), the electronic characteristics of GaAs are controlled by amphoteric 
doping and deviations from stoichiometry rather than by Inqiurlty segregation. 

also Initiated preliminary characterization studies on InP which 
very recently received considerable attention as a material competitive with 
GaAs In microwave and optoelectronic applications. Although such studies 
Introduce a new element to our program, we believe that they are essential for 
proper assessment of the relative potential of GaAs versus other compound semi- 
conductors in future applications. 


I 
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INTRODUCTION 

During the last three years we have made significant advancements in 

the electroepitaxial growth of GaAs» In the electronic characterization on a 

microscale, and in the study of growth-property relationships relevant to 

GaAs applications. We have also Initiated the study of melt growth of bulk 

GaAs crystals. Our program on "Crystal Growth of Device Quality GaAs in Space" 

* ^ 

has been widely exposed to the scientific and engineering comnunlty of leading 
industrial and educational institutions, and we have succeeded in establishing 
direct contacts which have proven extremely valuable in assessing the present 
status, major problems and future prospectlves of GaAs growth and applications. 

Table I summarizes the most important among the above developments. A 
more detailed discussion is given in our publications and 'earlier annual reports. 
Presently we will outline the most recert developments, i.e., since April 1, 

1979. 

CRYSTAL GROWTH 

Electroepitaxy 

Electroepitaxy represents a novel approach to liquid phase epitaxy (LPE) 

of semiconductor compounds in which growth is achieved by passing an electric 

current through the substrate solution Interface while the overall tsnperature 

of the system is maintained constant. 

The prospect of achieving precise control of the crystal growth process 

by simply controlling the current density stimulated extensive experimental 

and theoretical studies of electroepitaxy. 

In our earlier theoretical investigations of binary systems a model of 

(1-4) 

the growth kinetics and impurity segregation has been developed. This 

model defines the contribution of the Peltier effect at the solid-solution inter- 
face and that of solute electromigration to the overall growth process. Thus, 
it has been shown that electromigration of solute species to the growth inter- 
face is the primary means of mass truisport required for supersaturation and 
continued growth; the contribution of electromigration to growth is dominant 
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SEM-Cathodolunlnescence Advanced variable tenperature system has been 

set up for cathodolunlnescence microprof lllng 
of defects. Impurities & carrier concentration 
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Working Contacts Preliminary contacts have been estab- 

lished with certain industrial organi- 
zations In the area of GaAs characteri- 
zation, growth & device applications. 
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in the absence of significant convection in the liquid. The presence of 
convection in the solution enhances the contribution of the Peltier effect 
(cooling or heating at the growth interface) with pronounced concentration 
gradients which increase mass transport by diffusion. 

The ‘presence of electromigration and of the Peltier effect render 
electroepitaxial growth sensitive to the direction and density of the 
electric current, to the electrical resistivity of the solution, to the con- 
ductivity type (n or p) of the substrate and to the thickness of the substrate; 
as has been shown, these parameters provide a high flexibility in controlling 
and studying the growth process which is unattainable In standard LPE. 

Our present studies were addressed to three important aspects: 

(1) the feasibility of achieving high growth velocities comparable to sielt 
growth in GaAs; (2) the theoretical and experimental analysis of electroepi- 
taxial growth of GaAs-related multicomponent systems; (3) the development of 
new advanced systems for electroeptiaxial growth. 

Electroepitaxial Growth Velocity 

The experimental results of electroepitaxial growth velocity as a function 
of electric current density are given in Figures 1 and 2. Figure 1 presents 
the results obtained on substrates with (100) orientation commonly employed 
in LPE growth, while pig* 7 corresponds to substrates with (111) orientation. 
Only slight differences are seen between (111) and (100) orientations, and 
these differences are within variations expected from slight systematic dif- 
ferences in electron concentration of the employed substrate materials (leading 
to differences in Peltier effect contiibutions to the growth velocity). It 
is thus apparent that the electroepitaxial growth velocity is insensitive to 
substrate orientation. Furthermore, growth velocity increases linearly with 


ELECTROEPITAXIAL GROWTH VELOCITY (M"Vmln) 


“■ 8 



Fig. 1 ElecCroepltaxlal growth velocity versus electric 

current density for CaAs grown on n^type substrates 
(9.10^’cm“^) of (100) orientation. 
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•lectric current density, even for high current densities, i^lch provides 

direct evidence of the negligible role of Interfece kinetic phenooens. It 

should be noted thst the highest growth velocity obtained at 92S* for a cur- 

2 

rent density of ISO A/cn Is about 0.1 Inches/hr, l.e., comparable to tj^lcal 
growth rates of GaAs from the melt (0.3 to 0.5 Inch/hr) directional solidification. 
Electroepitaxy of Multicomponent Systems 

Our recent study of the electroepitaxy of multicomponent compounds 
has shown that the same phenomena which dominate the electroepltaxial growth 
of binary compounds are also of fundamental importance In this case. Thus, 
a theoretical electroepltaxial growth model of multicomponent systems was 
developed based on electromigration In the solution and on the Peltier 
effect at the growth interface. Quantitative relationships were derived for 
the dependence of the growth velocity and the cor^ositlon of the layers on 
the current density and other growth parameters. The model was successfully 
used for the analysis of experimental results on the electroepltaxial growth of 
Ga^_^^Al^As, and It enabled to explain the observed remarkable stabilization of 
the composition In the direction of growth under constant density. It was also 
found that the composition of the solid can be precisely varied by varying the 
current density. Experimentally determined cooposltlonal changes 
are shown In Figure 3 as a function of current density. The dashed line is 
computed from the theoretical model of multlcomponent^system electroepitaxy. 

It should be noted, In conjunction with the present results, that multicomponent 
systems play a key role In optoelectronic device structures and their tech- 
nology relies exclusively on epltsxlal growth. Standard LPE provides only 
limited flexibility in controlling the microscopic growth velocity, thus, 
precise control of the coaposltlon of ternary and quartemary cooqtound layers 
Is still a problem. The unique advantages of electroepitaxy in this respect 
are quite evident. 


Energy gap (eV) 
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Sew Experimental Sy»t«m» for El«ctro«pitaxy 

During the last year we i n 1 t i a ted the design and construction of 
two highly advanced systems for electroepitaxial growth of thick GaAs lav^rs 
and of improved quality heterostructure layers. 

The first of these systems incorporates a configuration of "rotating 
substrate" (see figure 4), originally proposed and studied in 
conjunction with electroepitaxy at the Institute of Inorganic Chemistry, 

Siberian branch of the Soviet Academy of Sciences. As shown schematically 
in figure A, the GaAs substrate is mounted on a rotating shaft and growth 
from equilibrated Ga->As solution is performed by passing electric current 
through the solution-substrate interface. The vertical motion of the graphite 
boat brings the substrate into contact with solution, makes possible pulling 
of the substrate with respect to the solution during the growth, and permits 
the separation of the substrate from the solution after growth is terminated. 

Rapid spinning of the substrate («100 RPM) provides a convenient means for an 
effective removal of the residual solution from the surface of the grown GaAs 
layer. We believe that Lhis approach (which in essence Incorporates all of the 
convenient features of Czochralski-type arrangement in epitaxial growth) will 
enable extension of the electroepitaxial growth to thick layer, and possibly 
to the growth of high quality, bulk GaAs crystals. Construction of the ap- 
paratus has already been completed. Extensive work on testing and calibrating 
the system is in progress; preliminary first growth experiments were carvied 
out and yielded promising results. 

The second electroepitaxy system, which is in the final states of con- 
struction, was designed for the growth of ultra-high quality layers and hetero- 
structures. It combines our own experience in current-contolled LPE with the advanced 
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know-how of industrial organizations (RCA, Nippon T & T, Fujitsu Labs) 
leading in the field of thermal LPE. 

Figures 5, 6, and 7 depict three essential subsystems of the new apparatus, 
i.e., the electronics system, the ambient gas system and the electroepitaxial 
boat, respectively. The electronic system utilizes a microprocessor which 
enables computer programming and controlling of the growth process. Detailed 
temperature control monitoring is considered a key factor for consistent and 

reproducible growth of a multicomponent layer. 

Significant improvements were also introduced in the ambient gas puri- 
fication system (figure 6) which, according to a new design. Includes H 2 
and N 2 purifiers and sorption pumps eliminating oil contamination of standard 
mechanical pumps. 

The multi-well boat (figure 7) was designed on the basis of the latest 
experience of RCA with modifications permitting the passage of electric cur- 
rent through the solution-substrate interface. This boat enables the growth 
electroepltaxlally, and/or thermally, of four layers of the desired composition 
and/or doping characteristics, and it will be employed for electroepitaxial 
growth of heterostructure layers and device structures. Growth experiments will 
be carried out in the very near future. 

Melt Growth 

Growth Apparatus 

During the last yet we completed the construction and testing of an 
advanced system for horizontal and/or vertical growth of GaAs. An overall 
view of this melt-growth apparatus is given in figures 8 and 9 in the horizontal 
and vertical operation mode, respectively. 
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Fig. 6 Gas purification system utilized in our new electro- 
epitaxy apparatus. 









- 20 - 


To satisfy unique thernal requirements the furnace system was constructed 
with four Independently controlled resistance windings, an Internal gas 
cooling system for enhanced heat extraction, and a heat pipe in the low 
temperature zone for Improved thermal stability and "flat" temperature profile. 
Power taps on the windings permit tailoring of the thermal profile In each 
zone. The furnace Is supported by a carriage assembly which rides on two 
parallel rails mounted on a two-piece main frame. The main frame permits 
positioning the system at any orientation from horizontal to vertical (figures 
8 and 9). Furnace motion is accomplished by a variable speed cable drive 
system pulling the carriage assembly. Drive rates range from 0.05 cm/hr to 
to 6 cm/hr. All of the support systems (power Inputs, gas feeds, sensor 
outputs and cable drive) have been designed to enable operation In any con- 
figuration without major modifications. 

A 15-channel chart recorder provides continuous monitoring of the 
thermal conditions throughout the system. Digital readouts enable visual 
monitoring of system temperatures. Four adjustable voltage sources permit 
high sensitivity thermal measurements to be made by offsetting the thermocouple 
potential. During the testing procedures en^> basis was placed on temperature 
control. All detectable sources of temperature instabilities In the sensor 
and control circuits have been systematically isolated and stabilized. Air 
circulation effects in the furnace tube which tend to flatten the temperature 
profile, were minimized by the use of zlrconla and alumina insulation between 
the quartz process tube and the furnace Interior wall. Long-term temperature 
stability Is typically better than -0.02*C. 

During the Initial growth experiments three major problems related to the 
actual growth process were Identified and solved. The first of these 
involved temperature control in the "cold zone". Since the temperature in 
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this cone determines the arsenic vapor pressure and thus the melt stoichio- 
metry and growth temperature. It Is essential that It be maintained very 
stable. It was found that heat transfer within the ampoule was acting as 
a "short circuit" path between the high and low temperature zones. The 
transfer* caused erratic temperature control In the low temperature zone leading 
to depletion of arsenic In the melt and collapse of the ampoule due to the 
decreased pressure Inside the ampoule. External control of the "cold zone" 

was abandoned In favor of control at the coldest point on the ampoule end 
wall. Implementation of this change led :o a thermal stability of -0.01 to 

io.02"c. 

The second problem was related to the accurate positioning of the Initial 
solid-liquid growth Interface without visibility through the furnace wall. 

To solve this problem a series of partial seed melting experiments were per- 
formed which led to the establishment of an experimental procedure which makes 
possible positioning of the Interface (through controlling of the exterior 
temperature profile) with a precision of -1 mm. 

A third practical problem was related to surface tension effects which 
complicate seeding of the crystal due to a tendency of GaAs melt to detach 
Itelf from the seed crystal. This problem was overcome by redesigning the 
boat configuration, which now permits reproducible seeding and growth of 
slngls crystals of GaAs. 

Finally, a new ampoule design was generated which permits the reuse of 
the most complex portion of the ampoule. This design reduces the time and 
costs of rebuilding the ampoule for each run. 

Characteristics of Melt-Grown GaAs 

The Initial melt-growth experiments were designed for testing the growth 
system rather than for studying growth-property relationships. Nevertheless, 
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these experiments were quite successful in yielding bulk GaAs with electronic 
properties (e.g. , nobilityt compensation and homogeneity) comparable to 
commercially available material of similar electron concentration. 

The electronic properties of two crystals are given in Table II. It 

should he noted that characterization of these crystals was based on the 

utilization of our advanced characterization facility, Thus, the compensation 

ratio was determined from carrier concentration and electron mobility with 

(\2 131 

our transport technique. ' * The homogeneity of the material was assessed 

with IR Scanning absorption. A typical electron concentration microprofile 

of the end portion (most inhomogeneous) of crystal 1 is presented in figure 10. 

It is seen that local changes of electron concentration are -202, i.e. , typical 

(13) 

for presently available melt-grown GaAs. 

A first derivative surface photovoltage spectrum characteristic of 
crystal 1 is shown in figure 11. The well resolved subbandgap Si-peak permitted 
the identification of Si as the dominating impurity. The density of disloca- 
tions was determined from etch-pit densities (figure 12) and was confirmed 
by SEM-cathodoluminescence scanning. The cathodoluminescence image of a 
sample prepared from crystal 1 is shown in figure 13. The dark spots in 
this image correspond to dislocation regions characterized by enhanced non- 
radlative recombination. The room temperature cathodoluminescence spectrum 
of the same sample (dashed and dotted lines correspond to different points 
on the sample) is shown in figure 14. Both the position of the peak (band-to- 
band transitions) and its half-width are consistent with literature data on 
cathodoluminescence characteristics of GaAs of similar electron concentration. 

The low electron concentration of crystal 2 permitted the determination 
of the deep level characteristics utilizing deep level transient capacitance 
spectroscopy. A typical DLTS spectrum is shown in figure 15, ~*Jidicating the 
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Fig. 13 Dislocation pattern of nalt-grown crystal 1 

as rcvaalad by SEH-cathodolumintscanca. Mots 
that Figs. 12 and 13 corraspond to neighbor ing 
aagaants of crystal 1. 
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Fig. 14 Room temperature cathodolumlnescence spectra 

of crystal 1 measured at two points separated 
by 25 yn (dotted and dashed lines, respectively). 
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Fig. 15 Deep Level Transient Capacitance spectrum of 

melt--grovm GaAs crystal 2 showing oxygen-related 
level at 0.81 eV below conduction band. 
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presence of an electron tap 0.81 eV below the conduction band. This trap 

coincides with a level commonly interpreted as related to oxygen. The con- 

16 *3 

centration of these traps is determined to be about 10 cm , i.e.» typical 
for melt-grown GaAs. The characterization effort outlined above constitutes 
an assessment of the electronic properties of Che first GaAs single crystals 
grown with our newly constructed melt-growth apparatus. These properties 
represent an initial reference point for our future research on growth-property 
relationships. 

CHARACTERIZATION 

Deep-Level Transient Spectroscopy 

Deep-level transient spectroscopy (DLTS) is presently used as a standard 
characterization technique of wide bandgap semiconductors. This technique 
is capable of displaying the spectrum of traps in a crystal in the form of 
peaks (emerging from a flat baseline) as a function of temperature. From the 
position and the height of the peaks one determines the energy position of the 
levels, the capture cross-section, and concentration. Since the DLTS method 
is sensitive, rapid, and easy to analyze, it is extremely useful as one of the 
basic techniques for establishing the relationships between growth conditions 
and the properties of GaAs. The block diagram of our recently constructed 
DLTS system is shown in figure 16. The basic capacitance measurements are 
performed with a modified 1 MHz commercial capacitance meter. The rate window 
is implemented by a double boxcar averager. In this approach, the emission 
rate, 62 * corresponding to the maximum height of the DLTS peak, is precisely 
defined by settings of gates in the boxcar averager. The values of e^ at dif- 
ferent gate settings and the temperatures corresponding to the peak maxima 
allow determination of the thermal activation energy of traps. The measurements 
for testing the system were performed on MS structures (Schottky barrier) 
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Fig. 16 Block diagram of our Deep Level Transient 
Spectroscopy system. 
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Fig. 17 Typical Deep Level Transient Capacitance spectrum 

of commercially available undoped (n-type, 6.10^^cm“^) 
melt-grown GaAs. DLTS peaks correspond to electron 
traps located at 0.25» 0.35» 0.6 and 0.81 eV below 
conduction band. Note that in our melt-grown GaAs 
(crystal 2, Fig. 15) only an oxygen related peak 
(0.81 eV) was observed. 
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prepared on commercially available melt-grown GaAs (n>*6xl0 cm )• A 

typical transient capacitance spectrum Is shown In figure 17. This spectrum 

reveals four peaks which can be considered as characteristic of the melt- 
( 21 ) 

grown GaAs. 

The*"DLTS technique was also employed for the Investigation of MOS-GaAs 
structures. Preliminary results have shown that measurements of transient 
capacitance spectra of MOS structures vs. filling pulse magnitude provide a 
convenient means for Identification of Interface states and bulk levels. 

Thus, as shown In figure 18, besides the DLTS peak corresponding to 0.8 eV 
bulk level (ok served also In MS structures) an additional peak Is observed 
at high filling pulse voltages. This additional peak can be accounted for 
on the basis of an interface state with a discrete energy at 0.65 eV below 
the conduction band edge. Further Investigation of Interface states is in 
progress. 

Derivative Photocapacitance Spectroscopy 
The capacitance measuring system discussed above has also been utilized 
in conjunction with wavelength modulation of the Incident light as a new 
convenient and simple approach to the determination of photoionization charac- 
teristics of deep levels. The advantages of wavelength modulation for 

studying deep levels have recently been demonstrated in our derivative surface 

(8 9) 

photovoltage spectroscopy as applied to GaAs MOS structures. * The newly 
developed derivative wavelength modulated photocapacitance spectroscopy is 
applicable to MOS as well as to structures simpler than MOS, e.g., metal- 
semiconductor Schottky barriers and p-n junctions. In the course of these 
studies the technique was successfully employed for the determination of 
deep levels in melt-grown GaAs and the Cr level in high resistivity GaAs. 
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Fig. 18 DLTC spectra of MOS structure (prepared on epitaxially 

grown GaAs by anodic oxidation process) showing interface 
states and bulk levels. Note the dominating character 
of DLTS peak due to the interface states for high voltages 
of the filling pulse. 
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SEM-Cathodolumlnescence 

A cathodo luminescence measuring system has been assembled to perform 
experiments in dispersive and non-dispersive modes at temperatures ranging 
from 80-700 K. The system consists of an optical microscope (incorporated 
into our 'scanning electron microscope apparatus), of a monochromator and a 
photomultiplier which serves as a detector of light emitted by the sample 
irradiated by the electron beam. 

The image shown In figure 13 is a non-dispersive mode representation of 
cathodoluminescence of Si-doped melt-grown GaAs (crystal 1) , discussed 
earlier. Emission spectra (dispersive mode) of the same sample were 
given in figure lA. The non-dispersive mode studies provide an instantaneous 
means for qualitative imaging of recombination Inhomogeneitles. Quantitative 
dispersive mode analysis enables the determination of carrier concentration 
microprof lies and/or microprofiling of radiative recombination centers. In the 
course of the present studies the relationship between the half-width of 
the emission peak and the carrier concentration was used to develop a simple 
procedure for quantitative electron concentration microprofiling of n-type 
GaAs. Typical results obtained with a commercially available melt-grown Te- 
doped GaAs are given in figure 19. It is seen that the crystal exhibits 
large carrier concentration luhomogeneities (n varies by a factor of 2) in 
agreement with our results obtained with IR absorption scanning. However, it 
Is also seen that cathodoluminescence scanning allows the identification of 
Inhomogeneitles with extremely high spatial resolution of about 1 Um. The 
advantage of cathodoluminescence over optical excitation techniques Is In 
this respect unquestionable. 

GROWTH-PROPERTY RELATIONSHIPS 

Electronic Properties of GaAs on a Microscale 


Utilizing our tunable laser IR absorption scanning technique we 
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achieved, for the first time, the. simultaneous microprof lllng of semiconductor 
free carrier, donor and acceptor concentrations. A number of techniques 
utilizing spreading resistance, Schottky barriers, EBIC, electroreflectance 
and IR absorption or topography have previously been developed for profiling 
the carrier concentration in semiconductors. However, these techniques 
essentially fall to distinguish between changes in Impurity concentration 
and in compensation ratio. The present method is an application of our 
recent quantitative studies of the free carrier absorption in GaAs and InP 
as a means of determining the compensation ratio. The principle of the method 
is understood by considering the absolute value and the wavelength dependence 
of the free carrier absorption coefficient are experimentally measured quanti- 
ties which contain only two unknown parameters, namely the electron concentra- 
tion n ■ and the compensation ratio 6 - (N^ and are the con- 

centrations of ionized donors and acceptors, respectively.) 

Microprofiles of free carrier, donor and acceptor concentrations obtained 
with melt-grown (horizontal Bridgman) Ge- and Si-doped GaAs are presented in 
figures 20 and 21, respectively. The Inhomogeneltles seen in these figures, 
although of an amplitude larger than average, are typical of all crystals 
obtained from several sources. Actually, as a general rule, all presently 
available melt-grown GaAs is highly compensated and exhibits pronounced varia- 
tions in free carrier concentration. 

A striking feature of the results of figure 20 is that the electron con- 
centration exhibits large fluctuations of a magnitude comparable to its average 
value, while the totc^. concentration of Ionized impurities remains ess€uitially 
constant. It is also seen in figure 20 that the spatial variations of donor 
and acceptor concentration are out of phase, l.e., maxima in coincide with 

minima of N, . 

A 


0 0.4 0.8 1.2 

Distance (mm) 


Fig. 20 Electron concentration and ionized impurity 
microprofilea of Ge-doped oelt->grown GaAa. 
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The variations of donor and acceptor concentration shown in figure 20 
cannot be explained on the basis of standard impurity segregation behavior 
controlled by the crystal growth velocity or the diffusion boundary layer 
thickness y and they must be associated with factors affecting the relative 
concentration of Ga and As vacancies or stoichiometry of GaAs. 

Similarly, the free carrier concentration profiles in Si-doped GaAs 
cannot be accounted for by standard segregation behavior. Thus, in figure 21 
the donor and acceptor concentration changes are not nearly so pronounced. 

The results of figure 21 were obtained from the same Si-doped crystal as those 
of figure Aa, but a different segment. In this case, the fluctuations of the 
free carrier, the donor and the total ionized impurity concentrations 

are of similar nature and magnitude, but the ionized acceptor concentration 
undergoes relatively small variations. 

It is apparent from figures 20 and 21 that fluctuations in the Individual 
profiles of carrier, donor, acceptor and total ionized impurity concentration 
can readily be mistaken as related to standard impurity segregation behavior 
as they resemble well established impurity segregation variations, particularly 
in elemental semiconductors. Thus, all of these microprofiles must be simul- 
taneously considered in assessing the origin of inhomogeneities In compound 
semiconduc tors . 

Such simultaneous analysis performed on the basis of IR scanning absorption 
■:>howed clearly that the electronic characteristics of GaAs on a microscale are 
controlled by amphoteric doping and deviations from stoichiometry rather than 
by inqiurity segregation. This conclusion is not necessarily limited to GaAs, 
but most likely it is valid for other semiconductor compounds as well. 

Our preliminary catbodolumlnescence study of InP has demonstrated striking 
differences between microscale variations of free carrier and impurity 
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concentrations. Typical results obtainad with a p-type, 10^^ca~'\ saaple of 
melt-gro%m InP are shown in figure 22, and they are labeled as I, II, and III. 

In each case t%ra spectral dependences are given as measured for the neighboring 
points (separated by about 25 yn). In figure 22 I, the band-to-band lunines* 
cence peak' (B~B) does not vary, while a significant increase in a lower energy 
peak involving inpurity level transitions (B-L) is observed. In case II, the 
band-to-band peak increases while the peak involving the impurity level de- 
creases noticeably. In case III, the band-to-band transition is enhanced, 
while the Impurity level transitions remain unchanged. Thus, it is evident 
that the spatial variations of band-tc- l and luminescence and of impurity level 
luminescence in InP are not interrelated. The sensitivity of the band-to-band 
transition to concentration of free carriers, and of the in^urity luminescence 
to impurity concentration leads to the conclusion that in InP, like in GaAs, 
the concentration of free carriers and the concentration of impurities does 
not follow the standard pattern established by crystal growth dependent segrega- 
tion behavior, but rather it is controlled by growth patterns affecting amphoteric 
doping and stoichiometry, which cannot be as yet unambiguously defined. 

Compensation of InP 

Ue have recently initiated preliminary characterization studies on InP 
which has been receiving Increased attention as a material competitive to 
GaAs in a number of applications. Thus, we have extended our theoretical cal- 
culation of electron mobility and free carrier absorption of GaAs to the case 
of InP, and we have performed rigorous experimental study of DC and high fre- 
quency AC transport phenomena. Detailed results of this study are given in a 

(14) 

forthcoming publication. ' ' 

It is, however, of importance to underscore that according to our own 
results and the analysis of literature data, melt-grown InF exhibits much 
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lower compensation than that encountered In melt-grown GaAs. This difference 
is evident in figure 23 where the normal range of N. , /n is plotted as a func- 
tion of electron concentration for commercially available high-quality GaAs 
and for InP. In the low electron-concenti&tjcn region the melt-grown GaAs 
is highly 'bompcusated, with a total concentration of ionized impurities exceeding 
by an order of magnitude the free carrier concentration. InP, on the other hand, 
even at low concei cr :itlons can be obtained by melt growth with a concentration 
of ionized Impurities of the same order of magnitude as the electron concentration. 
Tills marked difference is probably related. In part, to the lower melting point 
of InP (1062®C for InP and 1238®C for GaAs) and thus to toe lower concentration 
of native point defects; point defects in compound semiconductors are known to 
participate in the formation of compensation centers. Furthermore, the covalent 
radii of In and P differ by about 40%, whereas those ofGa and As are about the 
same; thus, a dopant impurity is more likely to occupy preferentially In or 
P sites in InP than Ga or As sites in GaAs (amphoteric or compensating doping). 
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APPENDIX 

PUBLICATIONS 

Reprints and preprints of papers which appeared in the literature or were 
submitted for publication since our last annual report are appended. They 
provide a more detailed account of some of the work discussed in the text of 
the present report. 
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ABSTRACT 

A comparative study was carried out on the defect structure and electronic 
properties of GaAs layers grown under similar conditions by electroepitaxy 
(electromigration controlled^) and by thermal LPE. It was found that the den- 
sity of certain microdefects» commonly present in GaAs layers, is significantly 
smaller (about an order of magnituae) in electroepitaxially ^Town than in 
thermally grown layers. It was also found that electroepitaxial layers exhibit 
greater carrier mobility and diffusion length than the thermally grown layers; 
in addition, electroepitaxiallv grown p-n junctions exhibited lower satura- 
tion current than the thermally grown junctions. The differences in the char- 
acteristics of two types of layers are attributed to corresponding differences 
in the nature of the driving forces lor growth. 


Growth of semiconductor compounds from solution 
can be initiated and sxistained by passing electric cur- 
rent through the growth interface while the over-all 
temperature of the system is kept constant. Electro- 
epitaxial growth has been successfully achieved with 
InSb (1), GaAs (2, 3), GaAlAs (4, 5), InP (6), 
HgCdTe (7), and garnet layers (8). The advantages of 
electroepitaxy in obtaining controlled doping (2, 9), 
improved electronic structures (9, 5), and in studying 
growth and segregation (10) have already been re- 
ported. A theoretical model of ^wth kinetics has 
been recently formulated (11) which defines the con- 
tribution of the solute electromigration and that of the 
Peltier effect (at the solid-solution interface) to the 
over-all growth process. According to the m^el, the 
contribution of electromigration is dominant in the ab- 
sence of convection in the solution, whereas the con- 
tribution of the Peltier effect can be dominant in the 
presence of convection. This model was found to be in 
excellent agreement with experiment data on the ki- 
netics of the electroepitaxial growth of GaAs from 
Ga-As solutions. 

Under electromigration control it has been shown 
that growth takes place under nearly equilibrium 
(isothermal) conditions (11) Consistent with this re- 
sult it has been found that the surface morpholo^ of 
electroepitaxial layers (attained under electromigra- 
tion control) do not exhibit terracing commonly ob- 
served on layers grown by standard (thermal) LPE 
(12). In the light of these results, the present study 
was undertaken and aimed at the comparison of the 
defect structure and electrical properties of GaAs 
layers grown by thermal cooling and by electromigra- 
tion controlled electroepitaxy under similar over-all 
conditions. 

Experimentoi 

Electroepitaxial growth of GaAs was performed in a 
standard LPE apparatus modified to permit controlled 
current flow through the solution-substrate interface 
(2, 13). Ga-As solutions doped with Ge (2%) were 
used; their height in the growth cell was from 10 to 14 
mm. Growth was carried out on a 0.5 cm< area of 
(100) oriented 300 ton thick GaAs substrates which 
were either Cr-doped (semi-insulating in room tem- 
perature; or Si-doped (n = 2 x 10** cm®). Electrical 
contact to the substrate was made with a 150 
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thick layer of Ga (13). To insure electromigration con- 
trolled grow^ convective interference was eliminated 
as described elMwhere (13). 

Prior to growth, the Ga-As solution was thermally 
equilibrated over a dummy substrate at 900*C. After 
bringing the solution into contact with the substrate, 
the temperature of the system was increased by 8*C 
to dissolve the substrate surface layer which might 
have been depleted in As during the thermal equili- 
bration period (14). 

Epitaxial layers ranging in thickness from 20 to 40 
;«m were grown by employing current densities of 2, 5, 
10, 20, 40, and 50 A/cm® with the substrate having a 
positive polarity or by decreasing the temperature 
10*C at a rate ranging from 0.1* to 2*C/min; this cool- 
ing range was chosen to obtain the same microscopic 
growth rates as those obtained by the employed cur- 
rent densities. 

The microscopic growth rate was detennlned by 
interface demarcation (15); current pulses (80 A/cm®) 
of 0.5 sec duration and at intervals of a few minutes to 
one hour were passed through the growth interface 
during thermal cooling or they were superimposed to 
the base current during electroepitaxy.® The micro- 
scopic growth rate was found to be essentially constant 
during the first hour of growth; only layers grown 
under constant growth rate were used in the present 
analysis. 

After growth the electric contact to the base of the 
substrate was removed by lapping. To determine the 
growth rate, the substrate and epilayer were cleaved 
and etched in AB etchant for about 20 sec; interference 
contrast microscopy was used to measure the spacing 
between interface demarcation lines. Samples for 
carrying out Hall-effect measurements and for p-n 
junction characterization were obtained by cleaving 
the substrate and the epilayer into 1x1 mm squares; 
electrical contacts were made with In or Sn by heating 
at 400*C for 10 min in a Ht atmosphere; layers grown 
on semi-insulating substrates were used for the Hall- 
effect measurements. Saturation currents were deter- 
mined from the 1-V characteristics of the diodes 
(p-type epilayers on n-type substrates). The diffusion 
length was determined from the dependence of the 
electron-beam-induced current (in an SEM) on the 
distance of the p-n junction employing a 30 kV ac- 
celerating voltage (16). 

•By chsnglns tbtlr tfuntlon and frsqttMief tt was aatabUahad 
that currant pulaac do not aXact tha ovarall growth precaas or 
the propertlaa of tha matarlal grows. 
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Crystallographic defects present in the epilayers 
were revealed by 5 min etching in nonstirred AB 
etchant at 10*-15*C. Because the presence of terraces 
on thermally grown surfaces interfered with the ob- 
servation of etch pits, these layers were lapped with a 
5 Mm particle diam abrasive and mechanical-chemical 
polished prior to etching. Electroepitaxially grown 
layers were etched without any pretreatment. To 
examine defect distribution in the growth direction, 
the epilayers were 1* angle-lapped, mechanical-chem- 
ical polished with (commercial) bleach solution for 40 
to 80 sec; they were then etched in AB etchant for 5 
min at 10*-15*C. This process was repeated a few 
times to obtain the statistical defect distribution 
through the epilayer. Determination of etch pit densi- 
ties on as-grown 'epitaxial layers and on angle-lapped 
surfaces showed that the lapping process did not intro- 
duce defects in the material; in addition the location of 
the substrate-epilayer interface could be determined 
with an accuracy of from 1 to 2 mHI by etching angle- 
lapped and cleaved surfaces. 

Microscopic defects , — Four types of crystallographic 
defects were observed in electroepitaxially and in 
thermally grown layers similar to those previously 
reported in the literMture. Typical examples of these 
defects in a thermally grown layer with an average 
velocity of 1 Mm/min (TC/min) are presented in Fig. 
1. The conical pits shown in Fig. la have been attrib- 
uted to dislocations which intersect the (001) surfaces 
at 90* (17). The flat-bottom rounded pits in Fig. lb are 
similar in shaoe to the Mnrer nit« observed 


in GaP (18); their origin is not clear; they are 
probably associated with agglomerates of defects or 
residual impurities (18). The shallow line defects in 
Fig. Ic are attributed to dislocation segments lying 
along the (100) surfaces (17). The boat-shaped defects 
(Fig. Id), oriented along <110> directions, are related 
to the presence of Ge, since they are not present in the 
substrates and they were not observed in epilayers 
which were Sn- or Te-doped; they could be attributed 
to either dislocation loops or to stacking faults (10, 
20). Apparently, the formation of these pits is due to 
the segregation of Ge atoms at crystallographic defects. 

The above defects were found to be homogeneously 
distributed in the planes of growth (excepting at the 
edges). The distribution of the conical defects in all 
epilayers was the same as in the substrates in the 
initial stages of growth; their density decreased by 
more than one order of magnitude after a few microns 
of growth. This flnding indicates that the propagation 
of dislocation from the substrates to the epitaxial 
layers is signiflcantly suppressed during LPE growth. 
Suppression of dislocation propagation (by a factor of 
about four), although not as pronounced, has been 
reported during LPE growth on low dislocation den- 
sity substrates (21). 

The distribution of the flat-bottom defects in the 
growth direction is shown in Fig. 2. The density of 
these pits is highest near the substrate-epilayer inter- 
face and decreases below their density in the sub- 
strates as growth proceeds. The increase of the flat- 
bottom pit density at the interface might be related to 
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Mgreftlion of defects at the surface of the substrate 
during the heating cycle. For example, As vacancies on 
the substrate due to As sublimation during heating, 
prior to growth, could enhance this effect (14). The 
density of the conical and the flat-bottom pits in the 
epilayers does not depend on the growth velocity and 
is the same in thermally or electroepitaxially grown 
layers. 

The density of the shallow line defects is shown in 
Fig. 3; it is highest near the substrate-epilayer inter- 
face; these microdefects were not observed in the sub- 
strates or in the epilayers beyond 1 of growth. 
Their density decreases with increasing current den- 
sity; in layers grown at 50 A/cm* current their density 
was found to bt about one order of magnitude smaller 
than in the epilayers grown thermally at comparable 
growth velocity. 

The largest difference between thermal LPE and 
electroepitaxy was found with regard to the boat- 
shaped pits. Their distribution is presented in Fig. 4; 
the defect density in thermally grown layers with the 
maximum and minimum velocity employed in this 
study (0.2 and 5 Mm/min) is shown together with 
results obtained for layers grown at 10, 20, and 50 
A/cm*. The boat-shaped defect density is significantly 
smaller in electroepitaxial layers thim in thermally 
obtained layers (grown with comparable growth 
rates) and is decreasing with increasing current den- 
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sity. Near the interface the density of the boat-shaped 
pits in the thermally grown layers is about two orders 
of magnitude greater than in the layers grown at 50 
A/cm* current. With continued growth the density 
of these defects decreases significantly in the thermally 
grown layers towards the values found in the electro- 
epitaxial layers. 

Electrical properties.^The carrier concentration, 
mobility, and diffusion length in the layers grown by 
electroepitaxy and under the same conditions by 
thermal LPE are presented in Table 1; each pair of 
layers grown at sWlar rates was obtai^ in a single 
run from the same solution. The variations of the car- 
rier concentration in the layers make a direct com- 
parison of the mobility in different layers rather diffi- 
cult; for this reason the values of the mobilities were 
normalized to a concentration of 2 x 10** cm** and are 
shown in Table L Normalization was obtained using 
the experimental relationship between mobility and 
hole concentration (9). It is seen that the normalized 
mobility is slightly but consistently higher in the 
electroepitaxial layers than in those grown with the 
same velocity by thermal LPE. The ratio of mobility 
in the electroepitaxial layers and in the corresponding 
thermally grown layers appears to increase with in- 
creasing growth velocity. 

The diffusion length in all layers grown by electro- 
epitaxy is higher than in the corresponding thermally 
grown layers. In both types of layers the diffusion 
length increases with increuing growth velocity. 
With increasing current density, the ratio of the diffu- 
sion length in electroepitaxially and thermally grown 
layers is increasing. It is apparent that the observed 
difference in diffusion length for the two types of lay- 
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tn cannot be attributed to chang eo in the carrier dif- 
fusion due to minority carrier mobility variations 
(22); it must oe due to changes in the lifetime. The 
increase of lifetime in the electroepitaxially grown 
layers compared with that in the thermally grown 
layers indicates that electroepitaxy provides material 
with a smaller concentration of recombination centers. 

The dependence of the saturation current of the p-n 
Junction on the growth velocity and/or current density 
for the layers grown by electroepitaxy and by thermal 
cooling is presented in Fig. 5. With increasing cur- 
rent density (up to 50 A/cm>) the saturation cur- 
rent of the grown junctions decreases by about one 
order of magnitude, while for the same changes of 
the growth velocity in thermal LPE the saturation 
current remains practically constant. The value of 
the saturation current increases with increasing con- 
centration of electrically active defects (in the junc- 
tion plane) leading to leakage current (21). The 
lower saturation current in the electroepitaxial junc- 
tions and its decrease with increasing current density, 
indicates that the density of electrically active defects 
in the electroepitaxially grown layers is lower than 
in the thermally grown layers, and it decreases with 
current density. This decrease in the density of elec- 
trically active defects in the junction plane contribut- 
ing to the leakage current, could be associated with 
the observed lower density of boat-shaped defects in 
the initial stages of growth in the electroepitaxial lay- 
ers than in the thermally grown layers. 

Summory ond Conclusions 

The micro-defect structure was studied in the lay- 
ers grown by electroepitaxy and by thermal LPE 
under the same conditions, i.t., same solutions and 
growth velocities. It was fbund that the density of de- 
fects related to dislocations lying in the plane of 
growth and the density of boat-shaped Ge-related de- 
fects, especially in the initial stages of growth, is lower 
in the electroepitaxially grown layers than in the lay- 
ers thermaUy grown under comparable conditions. It 
wu also established that their density decreases in 
both types of layers with increasing current density, or 
growth velocity. Higher values of carrier mobility and 
of minority carrier diffusion length and lower values 
of saturation current of p-n junctions were found in 
the electroepitaxial layers as compared to those in 


thermally grown laytri; thcM values were fbund to 
Improve with increasing ctirrent density. In view of 
this finding it appears likely that the obaerved defects 
are electrically active and that their decreased den- 
sity in electroepitaxially grown layers leads to im- 
proved electrical properties. 

The pronounced differences in the defect structure 
observed in the layers grown by electroepitaxy and 
those grown by thermal LPE cannot be attributed to 
specific growth conditions. Rather, they must be as- 
sociated with fundamental differences in the growth 
mechanisms involved in the two growth methods. 
Since quantitative models of nucleation and of defect 
fbrmation have not been established for these meth- 
ods, the present results will be considered in the 
light of some geneial features characterising the 
methods. 

In electroepitaxy (electromigration-controlled) eu- 
persaturation takes place in the immediate vicinity of 
the growth interface and growth is controlled by the 
rate of transfer of solute, under the electric field, to 
the interface. Under these conditions, growth takes 
place under nearly equilibrium (isothermal) condi- 
tions (11). On the other hand, in thermal LPE a super- 
saturation takes place in the entire solution. Thus, the 
As distribution in the solution and the interfkdal phe- 
nomena are different in electroepitaxy and in thermal 
LPE. Growth under a constant and uniform fiux of 
solute driven to the growth interface by an electric 
field must be less susceptible to defect fbrmation than 
growth under thermal LPE conditions. 
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ABSTRACT 

£xptrimenUl rttulU art presented showing that outdiffusion of recom- 
bination centers from the GaAs substrate into the epitaxial layer takes place 
during growth. Such outdiffusion decreases the carrier lifetime in the epi- 
taxial layer to much lower values than the radiative recombination limit. 
Furthermore, it introduces a lifetime gradient acrou the epitaxial laver which 
depends critically on the growth velocitv and thermal treatment High rates 
of growth (such as those attainable in eiectroepitaxy) and high cooling rates 
can minimize the adverse effects of normally available substrates on the epi- 
taxial layers; however, good quality substrates are essential for the con- 
sistent growth of device quality layers. 


It has been established that the performance of 
many compound semiconductor devices is adversely 
affected by nonradiative recombination centers. Their 
presence in the device active region leads to a decrease 
in the device efficiency [as in the case of solar cells 
(1)]. Furthermore, the recombination-enhanced diffu- 
sion of such centers (2) causes a slow degradation of 
the device performance [u in the case of semicon- 
ductor lasers (3)]. 

The nonradiative recombination centers in GaAs 
are usually associated with residual impurities such as 
oxygen (4), transition metals (5), and nonstoichio- 
m^c point defects (6). The origin of nonstoichio- 
metric defects is not fully understood, although in 
OaAs it has been attributed to the solidus of the phase 

* Itoctrocbemlcsl Society Active Meaher. 

*Prete*.t eourc«e: kCA keaeercli L.«ooretorlet, DavlS Semoff 
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diagram (7) and/or to As depletion of the substrate 
during the heating cycle prior to epitaxial erowth (^ ). 
Growth at low temperatures and backmeltinK uf the 
substrate prior to growth (9) tend to reduce the con- 
centration of nonstoichiometric defecU and improve 
the quality of the grown epitaxial layer. 

Poor quality GaAs substrate material is a source of 
defecU which can diffuse into the epiUxial layer during 
the growth process. The available melt-grown GaAs 
substrates are in most instances highly compensated 
(10). exhibit large carrier concentration inhomo- 
geneities (11). and contain high concentration of non- 
radiative recombination centers (12. 13). 

The present work is concerned with the outdiffusion 
of recombination centers from GaAs substrates into 
GaAs LPE layers. Minority carrier diffusion length 
and lifetime profiles obuined with EBlC-mode (SEM) 
measurements are employed for the study of such out- 
diffusion. Growth conditions required to minimize 
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tht outdllTuiion of rtcombinalion centtrt art pointed 
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Experimtnial 

Epitoxiol prowfh and lampU prtparotion*— LPE 
growth of OaAs was carried out employing clauical 
thermal cooling techniques and electromigration-con- 
trolled electroepitaxy (14, 15). Cd-doped (p & 2 x 
cm**) 900 ^m thick (100) substrates were used. 
Epitaxial growth was performed in a two-well graphite 
boat In one of the wells the substrates were ^ck- 
melted (about 25 ^) with an under-saturated solu- 
tion to remove the As depleted layer formed during 
the heating cycle. After backmelting the substrate 
was positioned in the second well containing 2.5g of 
an unioped Ga-^s solution with a GaAs source on 
top. In the electroepitaxy experiments layers ranging 
in thickneu from 10 to 100 were grown at 900* C 
on an area of 0.5 cm* at a rate of about 6 by 

passing an electric current of 60 A/cm* for a period of 
2-20 min. Growth was terminated by turning the 
current off. In the thermal growth experiments the 
growth was performed from equilibrated solutions 
by cooling from 910* to 900*C at a rate of l*C/min. 

In both procedures, after growth was completed, the 
system was kept at the growth temperature for time 

E »riods ranging from 1 min to 3 hr with the grown 
yer in contact with the solution. The experiment 
was terminated by quenchirj the system to room 
temperature with an initial cooling rate (from 900* to 
700*0 of about 70*C/min. The solution was left over 
the epitaxial layer or wiped by the moving slider 
prior to cooling. 

After cooling, if the exceu Ga-GaAs mixture was 
not wiped, it was removed from the epitaxial layer 
surface by boiling in HCl; the epilayer was cleaved 
and etched in AB etchant for about 30 sec. Interference 
contrast microKopy was employed to determine the 
thickness of the epitaxial layers. Whenever necessary, 
the part of the epitaxial layer grown during cooling 
to room temperature (about 5-10 Mm in thickness) was 
etched away in 5 parts H2SO4 4- 1 HtO ^ 1 Hi 02 . 
Ohmic contacts were soldered on the substrate and the 
epilayer employing Sn and In in an Hi atmosphexe. 

The epitaxial layers were n-type with a carrier con- 
centration of about 5 X 10>* cm**. Schottky barriers 
required for EBIC-mode measurements were made by 
evaporating aluminum or gold on the surfaces of the 
epitaxial layer. 

Minority carrier diffusion length and It/etime.^Mea- 
surements of the minority carrier diffusion length were 
performed utilizing electron beam excitation as shown 
schematically in Fig. 1. The diffusion length was ob- 
tained from the dependence of the electron beam-in- 
duced current (EBIC) on the distance between the 
generation position and the collecting Schottky barrier 
or the p-n junction (16). Reabsorption of photons re- 
sulting from radiative recombination in general intro- 
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duces uncertainties in the determination of the minor- 
ity carrier diffusion length. However, in the pre a tn t 
case, such uncertainties are of no consequence, since 
the lifetime of the minority carriers wu signiflcantly 
lower (up to two orders of magnitude) than the radia- 
tive recombination limit A 35 keV electron beam en- 
ergy was used in all experiments to minimize the effect 
of surface recombination on the measured diffusion 
length (17). 

The configuration shown in Fig. 1 made possible the 
determination of the diffusion length, Lp» of the minor- 
ity carriers in the epitaxial layer near the epilayer 
surface (with the electron beam at position 1 and 
vith current circuit 1), and near the substrate (with 
the electron beam in position 2 and with current cir- 
cuit 11); the minority carrier diffusion length wu also 
determined in the p-type substratu (with the elec- 
tron beam in position 3 and with current circuit 11). 
Valui s of lifetime, t, were estimated from the stan- 

( fcT \*'* 

— Mt J , where ^ is the 

minority carrier mobility and kT/e equals 0.026V at 
300*K. In p-type GaAs (substratu in the present case) 
the electron mobility can be noticeably smaller than 
the mobility of electrons in n-type material of similar 
free carrier concentration (19). Accordingly, in evalu- 
ating the electron lifetime in the substratu recently 
calculated theoretical valuu of electron mobility in 
p-type GaAs were used (18). In thue calculations the 
contribution of electron scattering by heavy holu and 
the difference in scruning energiu between holu and 
electrons have been taken into account. In the case 
01 holes their mobility u majority or minority car- 
riers is expected to be similar at room temperature and 
for free carrier concentrations below 10>* cm** (since 
screening effects and carrier-carrier scattering are not 
significant). Accordingly, in evaluating the hole life- 
time in epitaxial layers, the hole mobility value of 
Mg 300 cm*/V-sec wu used, i.e., the hole mobility 
in p-type GaAs wit> hole concentration of 5 x i0>* 
cm** (19), which is similar to the electron concen- 
tration in the epitaxial layers. 

Ruults ond DiKuuion 

Substrate.^The minority carrier lifetime in the 
substrate material was found to be about 10*^* sec, 
i.t., two orders of magnitude smaller than the lifetime 
value expected from band-to-band radiative recom- 
bination (20). It wu also determined that the sub- 
strate material is highly inhomogeneous with local 
lif;iime fluctuations excuding a factor of 4. These 
findings clearly indicate that the lifetime in the sub- 
strates is controlled entirely by nonradiative recom- 
bination and thus the substrate must contain a high, 
nonuniform, concentration of recombination centers. 

OutdiSusion oj re combination centers from the fub- 
ftrate.— The lifetime in the epitaxial layers wu maa- 
surtd u a function of position on a given plane parallel 
to the surface, u a function of distance from the sub- 
strate-epitaxial layer interface, and u a function of 
time the layer wu kept at the growth temperature. 

Typical results of the lifetime meuurements u a 
function of position on a plane parallel to the surface 
of the epitaxial layer are given in Fig. 2. The upper 
curve corruponds to positions near the surface of the 
epitaxial layer. The layer wu 40 Mm thick; it wu 
grown electroepitaxially and wu kept 20 min at the 
growth temperature (from the beginning of growth to 
the beginning of cooling). It is evident from Fig. 2 that 
lifetime inhomogeneities of a similar nature are pru- 
ent near the substrate-epitaxial interface and near the 
surface of the layer; however, near the surface the 
lifetime values are signiflcantly hijher. It should be 
noted that even the highut values of lifetime are well 
below the valuu (of the order of 10** sec) expected 
from radiative band-to-band recombination (20). 
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The Above behtvior of the minority cerrier lifetime 
can be understood if one considers that the substrate 
represents a source of recombination centers which dif- 
fuse into the epitaxial layer during the growth process 
and thermal treatment. On this basis, the difference 
between the lifetime near the epitaxial layer surface 
and the lifetime near the epitaxial layer-substrate 
interface should decrease with increasing exposure 
time of the system to the growth temperature. Simi- 
larly, the lifetime in the epitaxial layer should de- 
crease with increasing exposure time to the growth 
temperature. As is seen from Fig. 3, such behavior is in- 
deed found in electroepitaxially and in thermally 
grown layers. The lifetime in both electroepitaxially 
and thermally grown layers behaves similarly upon 
exposure to high temperature after growth, as it is 
dominated by recombination at outdiffused centers. 

The fact that recombination centers outdiffuse from 
the substrate into the epitaxial layers is further sup- 
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ported by the results shown in Fig. 4. Here the lifetime 
measured near the epitaxial layer surface is plotted u 
a function of thickness of the epitaxial layers, all of 
which have been exposed to a high temperature 
(growth temperature) for the same period of time 
(30 min). The layers were grown by electroepitaxy 
and at the same growth rate of about 0 ^im/min. It is 
seen tnat the liietime increases by about two orders of 
magnitude (near the surface) as the distance from 
the substrate (thickness of the layer) increases from 
0 toOO^m. 

If the simplified assumption is made that the life- 
time is inversely proportional to the concentration of 
the recombination centers, the data of Fig. I and 4 
can be used to estimate the diffusion constant of the 
recombination centers. Thus, treating the substrate- 
epitaxial layer interface as a limited diffusion source 
(21) the lifetime becomes 



where x is the distance from the substrate-epitaxial 
layer interface, t is the time, and D is the diffusion con- 
stant of recombination centers. From this expression 
the solid line in Fig. 3 is obtained* by taking x = con- 
stant and the solid line in Fig. 4 by taking t = constant. 
From these plots a value for the diffusion constant of 
the recombination centers of approximately 5 x 10~* 
cmVsec is obtained (900*0. The same value of 9 x 
10 cmVsec has been reported for point defect diffu- 
sion in GaAs (gallium vacancies) at 1000*C (22). 

It is important to note that the observed lifetime 
behavior which is consistent with the above diffusion 
model, cannot be explained on the basis of impurity 
segregation effects. A recent experimental and theo- 
retical analysis of impurity segregation in OaAs dur- 
ing electroepitaxy (23) has shown that the maximum 
changes of the dopant segregation coefficients do not 
exce^ 40^, ie., they are two orders of magnitude 
smaller than the changes required to account for the 
presently reported lifetime behavior. 
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Growth vtlocity and tha affacU of outdiffuiion,^Tha 
diffusion constant of impuriUtt and point defects de- 
creases exponentially with decreasing temperature, 
and thus lowering the growth temperature should 
drastically reduce the effects of outdiffusion from the 
substrate. However, there are limitations to the lower- 
ing of the growth temperature as the attainable 
growth rate decreases signiAcantly and single crystal 
growth becomes problematic. 

On the other hand, it is evident that the effects of 
outdiffusion can be reduced if the growth velocity, R, 
is much greater than the velocity of the diffusion front 
propagation vd- In a diffusion process a constant con- 
centration profllf can be approximated as x^/iDt = 
const.; consequently, vo = dx/df = const. \/B7i Thus, 
the condition Jl >> vo can be expressed as 

/ D V'* 

*»(t) 

or 

n»j (ibj 

Equation [2a] represents the case where growth is 
performed for certain time t and Eq. [2b] represents 
the case where the growth of a layer of a thickneu 
d is required. 

In the present case the outdiffusion constant of re- 
combination centers was found to be 5 x 10** cmVsec. 
Thus, according to Eq. [2b], for a layer 40 ^m thick the 
growth rate required to rHuce signiflcantly the effect 
of outdiffusion is R » 0.73 nm/min. In the present 
experiments the electroepitaxial growth velocity was 
about 6 Mm/min, which is high enough to satisfy the 
above condition for reducing signiAcantly the effects 
of outdiffusion of recombination centers. Accordingly, 
it is not surprising that, without a postgrowth heat- 
treatment epitaxial layers 40 >m thick were consist- 
ently grown with a lifetime of about 10*^ sec. 

The growth velocity in thermal LPE is limited by 
the slow solute transport to the growth interface and 
by the fact that the relatively small supercooling must 
be used to prevent spurious nucleation in the solution 
(24). These limitations are overcome in electroepitaxy 
(15). 

Sainmory ond Conclusions 

It was found that outdiffusion of recombirution cen- 
ters from a substrate to the epitaxial layer takes place 
during liquid phase epitaxial grovnh. iS-om the study 
of the time dependence of the diffusion process a dif- 
fusion constant of S x 10** cmVsec at 900*C was ob- 
tained for the outdiffusion of recombination centers. 
This value is the same as that reported for the diffusion 
constant of Ga vacancies in CaAs at lOOO'C. It should 
be noted that outdiffusion of point defects from the 
substrate has recently been suggested (25) on the 
basis of measurements on the distribution of residual 
deep levels in LPE CaAs crystals employing photo- 
capacitance in conjunction with step-etching. 

According to the present results, the recombination 
of minority carriers at the outdiffused defects con- 
stitutes the limiting factor for the minority carrier 
lifetime in the epitaxial layers. It has also been shown 
that the substrate-epitaxial layer interaction during 
the growth process can be reduced by increasing the 
growth velocity. Thus utilizing high growth velocities 
attainable in electroepitaxy, an improvement of as 
much as two orders of magnitude in the minority car- 
rier lifetime was observed. The present results could 
account for the improved characteristics of GaAs 
lasers grown under high growth velocities (26). 

Thus, the recombination characteristics of excited 
earriers in epitaxial layers are signiAcantly inAuenced 
by compositional and structural defects in the sub- 


strate. Since it is generally accepted that the quality of 
melt-grown GaAs substrates is poor (in terms of do* 
fects, residual impurities, and compositional inhomo- 
geneities) (12) improved GaAs substrates are oMential 
for the consistent growth of device quality GaAs LPE 
layers. 
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Liquid-phate electroepitaxy: Dopant aegregatlon 

J. LagowsW, L JMlrzebiki,** and H. C. Qatot 

DepQrtmentcfhiattriQhSckncfQndEngirtetring, Mauoehuiem ifatituit cfTfchi^oloiy, Combridgi, MQsmcbuMm02l39 
(Received 21 Miy 1979; accepted for publication 20 August 1979) 

A theoretical model it presented which accounts for the dopant segregation in liquid-phase 
electroepitaxy in terms of dopant transport in the liquid phase (by electromigration and 
diffusion), the growth velocity, and the Peltier effect at the substrate-solution interface. The 
contribution of dopant electromigration to the magnitude of the effective segregation coeffictent 
is dominant in the absence of convection; the contribution of the Peltier effect becomes 
significant only in the presence of pronounced convection. Quantiutive expressions which relate 
the segregation coefficient to the growth parameten also permit the determination of the 
diffusion constant and electromigration m^lity of the dopant in the liquid phase. The model 
was found to be in good agreement with the measured segregation characteristics of Sn in the 
electroepitaxial growth of OaAs from Oa- As solutions. For Sn in Oa- As solution at 900 *C the 
diffusion constant was found to be 4x 10* * cmVs and the electromigration velocity (toward the 
substrate with a positive polarity) 2x10*’ cm/s at a current density of 10 A/cm’. 


PACS numbers: 6S.SS. -f b, 6l.70.Tm 
INTRODUCTION 

It has been re^ .;; nixed for many years that the dopant 
distribution in semiconductors could be controlled by pas- 
sage of electric current through the solid-melt interface dur- 
ing crystal growth. Peltier cooling (or heating) in the vi- 
cinity of the interface ' and/or electromigration in the 
liquid phase have been considered as contributing factors 

to changes in the crystal growth velocity and/or the impuri- 
ty segregation. In the case of the growth of bulk single crys- 
tab from the melt the proposed approaches utilizing electric 
cunent have been proven of limited applicability because of 
convection instabilities in the melt and other interfering pro- 
cesses caused by Joule heating. 

In the liquid-phase-epitaxy (x PE] configuration, where 
the dimensions of the substrate and the solution are relative- 
ly small, the effecu of Joule heating resulting from dc can be 
minimized. Actually, it has been demonstrated that electro- 
epitaxy (LPE growth by passing dc through the substrate- 
solution interface) can provide control of the growth veloc- 
ity as well as of the dopant segregation and defect 
structure in the grown layers. In standard, thermal LPE, 
such control is still a complex problem. 

Recently a quantitative model has been developed for 
the kinetics of electroepitaxial growth. ' The model defines 
the contribution of the Peltier effect (Peltier cooling or heat- 
ing at the growth interface) and that of solute electromigra- 
tion to the overall growth proceu; it also provides expres- 
sions which relate directly the growth velocity to growth 
parameters. This model was found to be in excellent agree- 
ment with extensive experimental dau cm the electroepiux- 
ial growth of GaAs. ' 

In the present study of a theoretical mcxlel is presented 
for segregation of dopants in electroepitaxial growth. This 
model CO nbines the characteristics of growth kinetics and 
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segregation phenomena involving transport of dopants in 
solution (electromigration atid diffusion) and Peltier cooling 
(or heating) at the interface. The segregation model is ap- 
plied to the analysis of the experimental results obtained on 
the dopant segregation in electroepitaxial growth of GaAs 
from Sn-doped Ga-As solutions. 

ELECTROEPrr AXIAL GROWTH PROCESS 

In thermal LPE, growth is initiated by lowering the 
temperature of an equilibrated substrate-solution system by 
jr(i.e., Tq-^Tq - d^)andthusbringingaboutsupersatur- 
atiof * which leads to solidification taking place preferentially 
(although not always exclusively) on the substrate. The de- 
pletion (or accumulation) of the solute at the growth inter- 
face creates a concentration gradient in the solution and thus 
solute transport for the growth proceu is provided by diffu- 
sion. External control of diffusion is hardly possible and 
consequently control of the overall growth proceu is limited 
to the control of d T u a function of time. 

In electroepitaxy, the substrate-solution system is kept 
at the equilibnum temperature Tq , while an electric current 
is passed through the solution-substrate interface. ' Since the 
substrate (e.g., GaAs) and the solution (e.g., Ga-As solution) 
are dissimilar conductors, the electric current (in the appro- 
priate direction) leads to Peltier cooling at their interface 
which in turn leads to superuturation and growth u in sun- 
dard LPE. At the ume time, current flow (in the appropriate 
d ;ection) resulu in electromigration of the solute toward 
the substrate due to electron-momentum exchange and/or 
elect rosutic field forces. Such electromigration leads also 
to superuturation at the substrate-solution interface and 
thus to crysul growth. Thus, Peltier cooling as well as solute 
electromigration can contribute to the overall electroepiux- 
ial growth velocity v. It has been shown ' by taking into con- 
sideration solute transport through diffusion and electromi- 
gration that 

v^VfMEAi)^v,, ( 1 ) 
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where Vf is the contribution to the overall velocity by diffu- 
sion transport and Vg is the contribution by electromigra- 
tion; these velocities are expressed as follows; 


Ur = 


dC\ 


l>E 


C.-C, dT\, 
C, 


C.-C, 



(2a) 

(2b) 


where A is the temperature decrease due to Peltier cool- 
ing, C, and C| are the concentration of the solute in the solid 
and at the substrate-solution interface, respectively; D is the 
diffusivity of the solute in the solution; is the mobility of 
the solute under theHdectric field £; dC /dT ||. is the slope of 
the liquidus line. The function in Eq. (I) accounts for 
solute transport by processes other than diffusion (e.g., con- 
vective flow) and in metallic solutions of primary interest 
here = 1 in the absence of convection, and 

= {nDt ) ‘''Vi in the presence of convection, where i is the 
thickness of the diffusion boundary layer. 

In the absence of convection (4 = 1) the Peltier contri- 
bution, Vr, decreases with time, whereas is independent 
of time [^.(2b)]; thus after the initial transient period the 
overall growth velocity is practically constant, and it is con- 
trolled by the solute electromigration. Since Peltier cooling 
is usually small (of the order of 1 *C), C, in Eq. (2b) can be 
uken to be equal to Q , and the electromigration growth 
velocity becomes 


^ . ( 3 ) 

where Cq is the equilibrium solute concentration at the con- 
stant growth temperature Tq, 

In the presence of convection,/* is inversely propor- 
tional to 6 which in turn decreases with increasing convec- 
tive flow. Thus, in the presence of convection the contribu- 
tion of Peltier cooling to the overall growth velocity can be 
significant. 

The agreement of the above model with experimental 
results has been discussed in detail elsewhere. * 


SEGREGATION MODEL 

In considering dopant segregation in electroepitaxy it is 
assumed that the dopant concentration in the solution is 
small enough so that it has no appreciable effect on the phase 
diagram of the host system as well as on the electromigration 
and diffusion of the solute; i.e., it is assumed that the impuri- 
ty has no effect on the crystal growth velocity. 

Two representative cases will be analyzed; (1) convec- 
tion in the solution is \irtually absent; in this case the growth 
kinetics is controlled by electromigration and the growth 
velocity, i; = V£, is independent of time; (2) convection is 
present in the solution; in this case the role of the Peltier 
effect becomes significant, and after a transient period of 
ic^S VD growth proceeds essentially under steady-state 
conditions. Both of these cases are readily attainable experi- 
mentally as discussed elsewhere. 

In developing the model of dopant segregation in elec- 


troepitaxy, a similar approach is used to that employed in the 
theoretical treatment of growth kinetics. Thus, segregation 
will be related to mass-transport considerations involving 
electromigration and diffusion of impurities in the solution 
and the advancement of the growth interface. Diffusion of 
impurities within the solid and also interface phenomena 
(e g., adsorbed layers on the substrate) will not be taken into 
consideration. 

The incorporation of impurities into the solid will be 
described by introducing an interface segregation coe."^cient 
^0 = C„/C, . 0 where C„ and CJ, . o impurity 
concentration in the solid and in the solution at the interface, 
respectively. 

The pertinent equations are the transport equation and 
the relationship ensuring conservation of fluxes at the 
interface. 

The transport equation is expressed as 


^C, 




Bx 


BC, 
Bt ' 


(4a) 


where is the diffusion coefficient of the impurity in the 
solution, C, is the impurity concentration in the solution, x is 
the distance from the interface, v is the growth velocity, /i, is 
the mobility of the impurity (under an electric field) in the 
solution, and / is the time. The sign of the electromigration 
velocity /i,£ depends on the direction of current flow; in the 
present case/i, E is taken as positive when impurities electro- 
migrate toward the interface. 

Conservation of fluxes at the interface is expressed as 




(4b) 


The boundary conditions describing impurity concentration 
are: 


C.=C,, forr = 0 (4c) 

C, == for X = 00 (absence of convection) (4d) 
or 

C, = for (presence of convection). (4e) 

The initial temperature of the system is Tq\ upon cur- 
rent flow the temperature changes, due to Peltier effect, to 
r, = To 4 . Accordingly, upon current flow the inter- 
face segregation coefficient k changes from ^o) 

k '(T, ) = *0 -I- /dT I cases defined 

above (absence and presence of convection in the solution) 
the solution of Eq. (4) and, thus the effective segregation 
coefficient k^ = can be obtained in an analytical 

form. The Laplace transform procedure of solving Eq. (4) in 
the absence of convection is given in the Appendix. In the 
presence of convection, the steady state solution obtained in 
Ref. 5 will be adopted here. 


AbMoc# of convtction oloctromlgratlon-controllod 
growth 

According to the derivation given in the Appendix, the 
effective segregation coefficient is given by 
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P-a 


-erfc[-^(AO'''] 


2(0+^) 

+ — “exp[D,r(a*-^’)]erfc[a(X),0''*l|. (5) 

-♦“P I 




ID. 


It is of interest to note that if the impurity transport is 
not affected by the electric field (fi, = 0) and for =0, 
Eq. (S) reduces to 




Xexpl 


Equation (6) is identical to the standard expression for diffu- 
sion-controlled segregation under conditions of constant 
growth velocity. 

The complex form of Eq. (S) reduces to much simpler 
forms in the cases of high or low electric fields. For a high 
electric field, and under the assumption that^,£r 
>2(D,0 (5) reduces to 



Thus, the effective segregration coefficient increases expon- 
entially with time to a value = A o + Such be- 

havior is expected only in nonmetallic solutions (e.g., elec- 
troepitaxial growth of garnets '*) where the electric field can 
be of the order of V/cm. 

In metallic-type solutions (such as those employed in 
epitaxial growth of III-V compounds) the values of|i,£ lie 
between 10 “ ’ and 10 “ ^ cm/s; accordingly, the low-field 
condition,/i,£r is applicable. The approximate so- 

lution for low electric fields can be obtained from Eq. (5) by 
expanding the erfc and the exponential functions around 
zero [for small z, erfc I - (2/tr '^^)z] and by neglecting 
the products proportional to £" with n > 1 . Since A is 
proportional to the electric current density and thus to the 
electric field £, the linearization of Eq. (5) leads to the 
expression 



( 8 ) 


or 

*rfr/*o ^\^A^ ^A^'^A^ (9) 

where J and J ^are the contributions of the electromi- 
gration of impurities of the growth velocity and of the Peltier 
effect, respectively. 



Diatonca (/itn) 



FIG. 1. Time dependence of the effective segregition coefficient 
cilculated from Eqs. (9) and (10) for electroepitaxial growth of GaAs from 
Ga-rich solution at 900 *C under conditions of negligible Peltier cooling at 
the interface (J ^0). Dopant is electromigrating toward the substrate 
0) or away from the substrate (/i,£ < 0) with a velocity 
|^,£ I s 10 ' cm/s. Other parameters used in calculations: £>, s 10 
cm Vs, electromigration velocity of As; fi£ ^ 1 .3 X 10 ' cm/s, distance 
from the original growth interface is shown on the upper scale. 


It is apparent that if Acq > I* the effect of growth velocity 
can be more significant than the effect of electromigration 
since the ratio of \A ^ /A *'| becomes very small (the depletion 
of the impurity at the interface due to segregation becomes 
large compared to the amount of impurity transported to the 
interface by electromigration). Similarly, if < l.Ac^ does 
not necessarily increase with growth velocity, as it does in 
thermal LPE, because the sign of the electromigration term, 
A ^ is independent of the value of k^ but depends on the 
electromigration direction (i.e., it is positive if impurities mi- 
grate toward the growth interface and negative if they mi- 
grate away from the interface). Finally, the sign of the Peltier 
term, A ^ is determined by the direction of currrent flow (for 
a given conductivity type of substrate) and by the tempera- 
ture dependence of the segregation coefficient. 

The time dependence of 1 J ^ -I- J *’ calculated from 
Eqs. (9) and (10) is shown in Fig. 1 for electroepitaxial 
growth of GaAs from a Ga-rich solution at 900 *C and for a 
dopant electromigrating with a velocity |/i£ | = 10 ~ ^ cm/s. 
It is seen that in both instances, i.e., 1 ff^d Acq > 1, the 

value of 1 -f- ^ ^ -1- J is greater than unity (and thus k^ is 
greater than Acq) when impurities electromigrate toward the 
interface (^, £ > 0), and it is smaller than 1 (k^ smaller than 
Ato ) when impurities electromigrate away from the interface. 
Such behavior can be generally expected in electroepitaxy 
unless the electromigration velocity is lower than 
\v{ko — 1)1 and/or when the thermal contribution^''^ is sig- 
nificant. In most actual cases Acq is less than unity and also 
is expected to be much smaller than /i,£. (Note that v is 
typically one order of magnitude smaller than the electromi- 
gration velocity of the solute [see Eq. (3)]. Thus, the growth 


366 


J. AppI Phyi . Vol 51. No 1. January 1900 


Lagowaki. Jastrzebski, and Gatos 


366 



velocity contribution to dopant segregation in electroepitaxy 
plays a secondary role, while dopant electromigration con- 
stitutes the major contributint factor. 

For comparison with experiment it is convenient to in- 
troduce the average change of the segregation coefficient, 
Ak for a growth time, t defined as 

4*5 = - 00 

t Jl 

which according to Eqs. (8)-(IO) is 


4*5 

*0 


I 


*0 BT 


4F. 


(Ila) 


Pretenct of convtction 


The expression for steady-sute segregation when con- 
vection is present in the solution can be obtained by solving 
Eq. (4) with a boundary condition of Eq. (4e) for time 
t>6 VD, . By following the procedure given in Ref. 5 it can 
be readily shown that 


*f» 


— _j| . 




( 12 ) 


when Peltier cooling is neglected and Eq. (12) is 

the same as that originally derived in Ref. (S). In addition, if 
the impurities exhibit no mobility in an electric field (/i, = 0) 
Eq. ( 1 2) converts into the standard Burton-Prim-Slichter ex- 
pression ’’ for the effective segregation coefficient. 

For a low electric field approximation, valid in metallic 
solutions as those presently being considered, Eq. (12) 
becomes; 


= *o(l + 




It is seen that in Eq. ( 1 3) the impurity segregation coefficient 
can be expressed in a form similar to that of Eq. (9), appro- 
priately modified to include the thickness of boundary layer 
thickness, S. 

The calculated dependence of k^/k^ on the boundary 
layer thickness 6 is shown in Fig. 2, for the same parameters 
as used in Fig. 1 . It is seen that for larger 6 of the direction of 
change in segregation coefficient is determined primarily by 
the sign of electromigration velocity /i. £. However, the mag- 
nitude of d ^ -I- J decreases with decreasing 6 and thus for 
6—^ (strong convection) the Peltier term, A which is in- 
depndent of 8, becomes dominant. 



FIG. 2. Dependence of k,. /A„ on the boundary layer thickneu, d. calculat- 
ed from Eq (1 3) for electroepitaxial growth of GaAs from Ga-rich solution 
at 900 was taken as - 3*C; (l/k..KdA../dr) « 2x 10 ^ Other 
parameters are as in Fig. 1 . 


locity, and change in interface temperature; these contribu- 
tions depend linearly on the electric current, (b) The growth 
velocity component reflects a dependence of the segregation 
coefficient on growth velocity similar to that encountered in 
a standard segregation behavior, (c) Impurity electromigra- 
tion leads to an increase or decrease of the segregation coeffi- 
cient depending on the direction of electromigration rather 
than the sign of (/cq — 1). (t/ ) In the absence of convection in 
the solution, the contributions of electromigration and 
growth velocity to k^ depend on time as t the Peltier 
contribution is independent of time, (e) The presence of con- 
vection leads to steady-state contributions from electromi- 
gration and growth velocity which are proportional to the 
boundary layer thickness, 6; the Peltier contribution is inde- 
pendent of 6; accordingly, with increasing convection (de- 
creasing 6 ), the relative contribution of Peltier effect is 
enhanced. 

In the present treatment of impurity segregation in elec- 
troepitaxy from equilibrated solutions two types of tran- 
sients were not considered; (a) the relatively fast transients 
associated with the establishment of a new interface tem- 
perature under the influence of the Peltier effect; typical time 
constants of these transients are of the order of 1 s; (b) the 
relatively slower transients (with time constants of the order 
of 1 min) during which the growth velocity relaxes to a value 
(see Ref. 8). These cases require a numerical approach 
and they will not be pursued at present. 


Oantral ftaturat of Impurity Mgrogatlon In 
•loctroapltaxy 

The essential features of impurity segregation in elec- 
troepitaxy from metallic solutions can be summarixed as fol- 
lows: (a) The segregation behavior can be analyzed in terms 
of the additive contributions of electromigration, growth ve- 


EXPERIMENTAL PROCEDURE 

Electroepitaxial growth experiments were carried out 
from Sn-doped Ga-As solutions in a standard LPE appara- 
tus (horizontal graphite boat) modified to permit passage of 
electric current through the solution-substrate interface as 
described elsewhere. 
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The effect of electric current on changes in the impurity 
segregation coefficient was studied in the absence and in the 
presence of convection in the solution. Thus, growth was 
performed in a configuration in which convectional flow in 
the solution was either essentially eliminated as described in 
Ref. 16 or was intentionally introduced by superimposing a 
2.S *C/cm horizontal temperature gradient across the solu- 
tion. The magnitude of convection was further varied by 
varying the solution height (Grashoff number). The value of 
the solute boundary layer thickness was obtained from an 
analysis of the growth velocity using the procedure de- 
scribed in Ref. 8.^ * 

Semi-insulating (at room temperature) Cr-doped sub- 
strates were used; their thickness ranged from 2S0 to \500 
/im. The desired value of the interfacial temperature de- 
crease (due to Peltier cooling), for a given current, was ob- 
tained by selecting the proper thickness for the substrate. As 
reported in Ref. 8 varying the thickness of Cr-doped sub- 
strates from 250 /im to 1.5 mm, J is varied from 0.5 to 
3.0 *C for the same current density of 10 A/cm ^ In order to 
remove the As-depleted substrate surface formed during the 
heating cycle, the substrates were back-etched prior to 
growth by a 3 *C undersaturated solution. 

The electronic characterization of the grown layers in- 
cluded the determination of carrier concentration on ma- 
croscale and microscale and the compensation ratio, 

6 = /Np , where and Np are the concentration of the 

ionized acceptors and donors, respectively. The compensa- 
tion ratio in LPE layers of GaAs ranges from 0. i to 0.8. Sn 
introduces shallow donors and compensating deep acceptors 
with 6 typically - 0.3 in layers grown at 700 *C. Accord- 
ingly, for a reliable analysis of impurity segregation, knowl- 
edge of the d values is required. In the present study, 0 was 
determined by recently developed theoretical and experi- 
mental methods based on electron mobility and free carrier 
absorption. When the electron concentration, n, and the 
compensation ratio, 0, are known, one can determine the 
concentration of ionized donors, Np, and acceptors, . 
Since 

e= — and n^No-N^, (I4a) 

No 

it is apparent that 

Np = — - — and = n — (14b) 

In nondegenerate material the value of Np, given by expres- 
sion (14b), can be used as the total concentration of donor 
impurities. In degenerate material (n > 6x 10 cm ‘ ^) a 
correction must be introduced accounting for the occupa- 
tional statistics of impurity levels. 

Charge carrier concentration profiles in the electroepi- 
taxial layers perpendicular to the growth direction were de- 
termined by ir absorption scanning which provides a 
spatial resolution of about 20 /im. 

For the analysis of the Peltier effect contribution to the 
segregation coefficient, knowledge of the temperature de- 
pendence of Au is required. The values of for Sn were 
estimated from measurements on epitaxial layers grown 


thermally (no electromigration involved) by ramp cooling 
from a temperature of Tp -f 2.S *C to a temperature of 
To — 2.5 *C The cooling rate was varied betwen 0. 1 and 
3 *C/min; this rate was found to have no effect on the segre- 
gation coefficient during thermal growth. The values of kp 
obtained for different temperatures were used to determine 
(1/^0 /dT) which at 900 *C was estimated to be 
+ 0.01 *C “ * (at 900 *C kp was aout 5x10“ % 

In the analysis of the experimental data the average val- 
ue of the effective segregation coefficient will be used for 
donor-type impurities: 


^ 1 _ 

*0 *0 Sp(J^Q) 


(15) 


where Np{J) and Np(J = 0) are the average donor concen- 
tration in an electroepitaxially grown layer and in a thermal- 
ly grown layer, respectively. (7 denotes the current density). 


RESULTS AND DISCUSSION 
Abtanco of convoction 

It was found that the carrier concentration in the elec- 
troepitaxial layers, grown from solutions in which convec- 
tion was essentially eliminated, was very uniform in the di- 
rection normal to the growth direction. Typical electron 
concentration profiles obtained by ir scanning are shown in 
Fig. 3. The layers were grown from 3.5% Sn-doped solution 
at 950 *C with current densities of 1 5 and 50 A/cm ^ applied 
for a period of 30 min. Except for the layer edges (not shown 
in Fig. 3) the electron concentration in the growth plane does 
not change by more than 1%. 

The effect of current density on the change of the aver- 
age segregation coefficient is shown in Fig. 4. Growth was 
carried out for 20 min from 2.5% Sn-doped solutions at 
900 *C. Two substrate thicknesses were used, 300 and 1000 
/im; the corresponding values of were — 2.5 and 
— 8 *C, respectively, for 40 A/cm ^ It is seen that the segre- 
gation coefficient increases linearly with current density. 
For the thicker substrate, the values of Ak ^/Aco are slightly 
lower due to the higher contribution from Peltier cooling 
[last term in Eq. (11)]. However, the difference between the 



FIG. 3. Electron concentration proSle pcrpeni!icuUr to the growth direi 
tion for Sn-doped elect roepitaxial layeri (see test) 
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FIG. 4. Relative change of the effective segregation coefficient Ak 
versus current density for Sn-doped layers grown electroepitaxially at 
900 *C in the absence of convection on Cr-doped substrates 300/im thick ({) 
and 1000 /im thick (^). respectively (corresponding values of were 
- 2.S and — 8 *C). Growth was earned out for 20 min. 


results obtained with thick and thin substrates (in spite of the 
difference in J ) is rather minor, which indicates that seg- 
regation is controlled by impurity electromigration and the 
growth velocity. 

The electroepitaxial growth velocity contribution to the 
results of Fig. 4 estimated from expression (1 1) does not 
exceed 15%. It is thus apparent that the observed segrega- 
tion coefficient changes induced by the electric current are 
due to electromigration of Sn in the solution. The increase of 
segregation coefficient indicates that electromigration of Sn 
is directed toward the positive polarity of the substrate, 
which is consistent with reported electrotransport direction 
of Sn in liquid gallium. Using expression (1 1) a value of 
l/D'/^yi,E -b i<l — *0 ]:^3.8x 10“ ^ s“ is obtained 
from the slope of the line of Fig. 4 for a current density of 10 
A/cm ^ 

The experimental results of the dependence of the aver- 
age segregation coefficient on time are given in Fig. 5. 
Growth was performed with a current density of 30 A/cm ’ 



FIG S. Growth lime dependence of the relative change in Kgreption coef* 
flcient k ^ /k„ for layen electroepitaaially grown with current density 30 
A/cm in the absence of convection The substrate thickness was 3(X)/im. 
other parameters the same as in Fig 4 Dotted line is calculated from Eq 
(1 l)(see test) 


and with substrates of 300/im thickness. Other experimental 
conditions were identical to those employed in the experi- 
ments of Fig. 4. Using the value of 
\/D]'^[^,E -h t<l - iko j2s3.8X 10“ » s-“ (as deter- 
min^ above) the theoretical dependence of AkH/ko on 
time was calculated from expression (1 1). The calculated 
curve is shown in Fig. S. It is seen that the theoretical behav- 
ior is in good agreement with experiment. 

On the basis of the above results it is apparent that the 
dependence of the effective distribution coefficient on cur- 
rent density for constant growth time as well as the depen- 
dence of the effective distribution coefficient on time for con- 
stant current density are quantitatively accounted for by the 
present model of segregation in electroepitaxy. It is also ap- 
parent that the electric current effect on Sn segregation in the 
GaAs : Ga-As solution system is due almost entirely to Sn 
electromigration. 

Prattnet of convoction 

The electroepitaxial layers grown in the presence of 
convection were found to exhibit a pronounced variation in 
electron concentration in the direction perpendicular to the 
growth direction. Typical results obtained with ir scanning 
are shown in Fig. 6. The layers were grown at 950 *C with a 
current density of 30 A/cm ^ from 3.5% Sn-doped solutions; 
the solution heights were 0.6 and 1.2 cm, and the substrates 
were 300 /im thick. 

For relatively small convectional flow (0.6 cm solution 
height) the variation in impurity concentration across the 
layer is about 10% and increases to 20% with increasing 
convectional flow (1.2 cm solution height). These changes 
reflect an increase in the variation of the solute boundary 
layer thickness along the substrate with increasing convec- 
tional flow in the solution. It has been reported earlier that 



FIG 6 Electron concentration profile of electroepitaAial layen grown in 
the presence of convection The layen were grown a 950 XT with a current 
density of 30 A/cm ^ Cr-doped substrates were used with a thickness of 300 

fim 
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FIG 7 Changr in Sn trgrefation coefficient JA.n /ik., as a function of cur- 
rent density in the presence of convection (see leal). 


vtriations in solute-boundary layer thickness can also lead to 
variations in the thickness of the grown layers. 

The experimental results of the steady-state segregation 
coefficient as a function of current density are shoWn in Fig. 
7 for two values of solution height, 0.6 and 1 .2 cm; the corre- 
sponding values of 5. as estimted from the growth velocity, 
are 0.1 1 and 0.05 cm. respectively. Growth was performed 
on l.S-mm-thick substrates at 900 *C from 2.59^ Sn-doped 
solutions. For comparison, the changes of the average segre- 
gation coefficient obtained under convection-free conditions 
is also shown in Fig. 7. 

Thus, consistent with the theoretical model, the magni- 
tude of segregation coefficient changes in electroepitaxy is 
noticeably reduced in the presence of convection. It should 
also be noted that with increasing convection (decreasing 6 ) 
the current-induced change of the segregation coefficient be- 
comes negative (Fig. 5). Such behavior is readily understood 
in the light of expression ( 1 3): The contribution from electro- 
migration and from growth velocity decreases with decreas- 
ing 6, to the extent that the thermal contribution J ' be- 
comes dominant. However. J ^ is negative (J < 0). 

From Eq. (13) and the experimental data of Fig. 7, 
^,E/D, can be obtained. Thus, considering that (1/Aco) 

(dk /dT)i^ 10 - M/ C. J /y-0.3 C/A cm';« = 0. 11 and 
0.05 cm. the value of/i,£ /D, \ cm is ob- 

tained from the slopes of the two curves of Fig. 7. 

This value, together with the parameters determined 
earlier, yields D, X 10 " ^ cm Vs and Ai.f | ,o A/cm* 
ss2.3 X 10 ~ ’ cm/s for Sn in Ga-As solution. Both values 
seem quite reasonable. 

The effect of increasing convection flow on the changes 
of the segregation coefficient is shown in Fig. 8. The layers 
were grown under constant current density (40 A/cm on 


300-/im-thick substrates (^ « — 3 *C) at 900 *C from 

2.3% Sn-doped solution with heights ranging from 0.3 to 1 .2 
cm, corresponding to a range of 6 from 2.0 to O.S mm. The 
dotted line represents the dependence calculated from the 
theoretical model [Eq. (13)] using parameters as determined 
above. The contribution of the Peltier effect ot the segrega- 
tion coefficient changes is small, compared with the contri- 
bution of electromigration. However, the effect of electromi- 
gration decreases (consistent with theory) with the 
decreasing thickness of the solute-boundary layer (increas- 
ing convection in the solution); thus, for sufficiently large 
solution heights a change in the sign of is observed (see 
also Fig. 9). 

The effect of Peltier cooling on the changes of the segre- 
gation coefficient in the presence of strong convection in the 
solution is shown in Fig. 9. Growth was performed at 900 *C 
with a current density of 40 A/cm ^ from 2.3% Sn-doped 
solution with a height of 1.2 cm (£t:^.03 cm) on substrates 
with thicknesses ranging from 200 m to 1 .3 mm, correspond- 
ing to a range of J values from — 2 to — 13 *C. As pre- 

dicted from the model, for small values ofAT^, the change of 
the segregation coefficient is controlled by electromigration, 
but with increasing substrate thickness {A T^) the contribu- 
tion of the Peltier effect becomes comparable to that of elec- 
tromigration; for thick substrates ( > 1 mm) the Peltier effect 
becomes dominant and leads to a change in the sign ofAk^. 

SUMMARY 

Dopant segregation in liquid-phase electroepitaxy was 
investigated employing GaAs epitaxial layers grown from 
Sn-doped Ga- As solutions over a broad range of experimen- 
tal conditions. The results were analyzed on the basis of a 
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Chong* in interfoc* ttmp. *Tp(*C) 



FIG. 9. Chinge in Sn segrefition coefficient versus thickness of 

subetnie lor change in interface temperature. A caused by Peltier effect) 

for layer growth in the presence of convection. 


theoretical model which relates quantitatively the dopant 
segregation coefficient and growth parameters such as tem- 
perature, convection flow in the solution, and current densi- 
ty or growth \elocity. In the theoretical analysis dopant 
transport in the liquid phase by diffusion as well as by elec- 
tromigration was taken into consideration. 

Consistent with the theoretical model it was found that, 
in the absence of convection in the liquid, the effective segre- 
gation coefficient ofSn increases with increasing current 
density and also increases with increasing time (for constant 
current density). This behavior is encountered when the 
electromigration term is positive, i.e., when electromigration 
of the dopant proceeds toward the growth interface. The 
effective distribution coefficient decreases with increasing 
current density if the eiectromigration term is negative. 

In the presence of convection, the effective segregation 
coefficient decreases with increasing magnitude of convec- 
tive flow (as assessed by the height of the solution) for a 
constant current density even below the value of the equilib- 
rium segregation coefficient. Since, in the presence of con- 
vection, the contribution of the Peltier effect to the growth 
process becomes significant, and since Peltier cooling in- 
creases with increasing substrate thickness, the value of the 
effective segregation coefficient depends on the substrate 
thickness for a constant current density and solution height. 

From the epitaxial data the diffusion coefficient and the 
eiectromigration velocity of Sn in Ga-As solution at 900 *C 
were determined on the basis of the theoretical model. The 
values obtained were D = 4 x 10 ’ * cm Vs and/i£ | 
^2x10”* cm/s. 
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APPENDIX: SEGREGATION OF IMPURITIES 
MIGRATING IN AN ELECTRIC FIELD 

In most instances electroepitaxial growth proceeds 
with a growth velocity, u. which is independent of time. Un- 
der such conditions, the impurity distribution in the liquid 
phase C. (x.r ) can be obtained in an analytical form. In the 
absence of convectional flow, the pertinent set of equations 
becomes: 


diffusion and eiectromigration 



II 


dC, 

dx 



(AI) 


where II = 

conservaton of fluxes at interface (x = 0) 


vC,, — uC, 


.-0 ‘ dx 


MmO 


boundary conditions: 

C. = C„ forr = 0 

and 

C, = C„ for X = 00 for any t. 

By introducing a quantity C * defined as 

Eq. (AI) reduces to a standard diffusion equation 

n = i£l 
' dx‘ dt 

and the equation of fluxes becomes 




where 


A = (Aco-I)t' + 


v-n,E 


and 


(A2) 

(A3) 

(A4) 

(A5) 

(A6) 

(A7) 


^0 — |« - O' 

Equation (A5) will be solved using the Laplace transform 
method. The Laplace transform of C * is 

C’= JT exp( — pr)C*<//. (A8) 

Multiplying both sides of Eq. (AS) by exp( — pi) and inte- 
grating over time one obtains 

D,^[ txp{-pi)C'dt 
dx^ Jo 


exp(-pr)C* 1° +pj^ exp(-pi)C'Jt. (A9) 


However, cxp( — pi )C* vanishes for I = oo and for I = 0 it 
equals C„ t%p[(u/lD^)x\, Thus (A9) leads to 

Similarly, the equation of fluxes (A6) becomes 

D,^^\ = /IC*| (All) 

dx 
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The solution of Eq. (A 10) it 

f *-yle*p[-(p//),)''‘*l + 


C.exp(«/2/)> 


p — «V4i), 

(AI2) 

where the constant A is obtained from the boundary condi* 
tion(All) 

C^{y-2R) I 


A = 


2i>? (p/D,-0^)[a + (p/D,y'^y 


(AI3) 


where a = R/D, and P = u/2D, . Thus 
C^(u-2K‘)exp[-(p/A)'^*l 
2Df(p/D, -^‘)(a + 0»/A)''*) 

C*exp(^x) 


(A14) 


Expansion of Eq. (A14) in partial fractions leads to 

1 exp[-(p/A)''^*] 


_ _ C^{M-2R){D,r^ ( 
ID] V 


2ia-¥p)P {p''^-PD]'^ 

1 exp[-(p/A)''*Jt] 


Wifi -a) 

1 exp[-(p//),)'^xl ^ 

exp(-(ux/2Z).)l 


+ C« 


(p-fi^D.) 


(AI5) 


For all components of Eq. (A 1 5) the inversions of Laplace 
transforms are readily obtainable.^' Using the tabulated 
forms of the inversions of corresponding Laplace transforms 
and the relation (A4) one obtains from Eq. (A IS) the solu- 
tion for the impurity distribution C,(jc,f ) as follows; 

(fi-a) np(-2xfi) 

2{a+fi) 


Xerfc 

i 




a + fi 


exp[x(a - /?) + D,t{a^ - 


The effective segregation coefficient and the interface 
segregation coefficient, Aq defined as 


S Hi. 

*wr = 77- *« = 


(AI7) 

C.(0.r) 

Thus Uiking Eq. (AI6) for x » 0 and multiplying both sides 
by k(,k^/C„ one obtains 

A- 


*o(l-ierfcl^(AO‘'*l 


fi-a 


2(a+fi) 

+ — ^«*i 

a^^■fi 


ttfcl -fi(P,t)''*] 


p[DMa* -fi*)] erfc(a(AO'^l) • (AI8) 


Equation (A18) describes the segregation behavior for iso- 
thermal electroepitaxial growth at equilibrium temperature 
Tq. If, due to the current flow, the interface temperature 
changes (Peltier effect) by the interface segregation co- 
efficient ko in Eq. (AI8) be replaced by A: o * ^ 
^dko/dr\rAT^. 
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Electroepitaxy of multicomponent syetemt: ternary and quarternary 
compounds 
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A theoretical model it pretented which accounts for the electroepitaxial growth kinetics and 
composition of multicomponent compounds in terms oT mau transport in the liquid and phase 
diagram relationships. The mass transport in the interface it dominated by electromigration in 
the absence of convection and by diffusion in the pretence of convection, llie composition of the 
solid it controlled by the Peltier effect at the growth interface and by the diffusion and mobility 
constantt^of the solute components and the growth velocity (current density). Thus, for a given 
solution composition, the composition of the solid can be varied by varying the current density. 

For a given current density the composition remains constant even in the case of rdativdy thick 
epitaxial layers. All aspects of the model were found to be in good agreement with the growth and 
composition characteristics of Oa. _ , AI, As layers. 

PACS numben; 6I.S0.a 8l.l5.Lm. 64.73. + g. 7l.2S.Tn 


L INTRODUCTION 

Electroepitaxy represents a novel approach to liquid 
phase epitaxy (LPE) of semiconductor compounds in which 
growth is achieved by passing an electric current through the 
substrate solution interface while the overall temperature of 
the system is maintained constant.' 

The feasibility of achieving precise control of the crystal 
growth process by simply controlling the current density has 
stimulated extensive experimental and theoretical studies of 
electroepitaxy.^** 

In recent theoretical investigations of binary systems a 
model of the growth kinetics’ and impurity segregation* has 
been developed. This model defines the contribution of the 
Peltier effect at the solid>solution interface and that of solute 
electromigration to the overall growth process. Thus, it has 
been shown that electromigration of solute species to the 
growth interface is the primary means of mass transport re- 
quired for supersaturation and continued growth; the contri- 
bution of electromigration to growth is dominant in the ab- 
sence of significant convection in the liquid. The presence of 
convection in the solution enhances the contribution of the 
Peltier effect (cooling or heating at the growth interface) 
which otherwise is of importance only in the initial stages of 
growth. The primary reason for this enhancement is the for- 
mation of a boundary layer (at the interface) with pro- 
nounced concentration gradients which increase mass trans- 
port by diffusion. 

The presence of electromigration and of the Peltier ef- 
fect render electroepitaxial growth sensitive to the direction 
and density of the electric current, to the electrical resistivity 
of the solution, to the conductivity type (h or p) of the sub- 
strate, and to the thickness of the substrate; as has been 
shown,’ these parameters provide a high flexibility in con- 
trolling and studying the growth process which is unattaina- 
ble in standard LPE. 

*^ennmcni iddroc InMituic oTPhyiia. PolUh Academy of Sciences. 

Wtnaw, Poland 


Electroepitaxial growth of multicomponent systems 
(i.e., Oa, _ , Al. As), with a remarkable stabilization of the 
composition in the direction of gn vth and improved prop- 
erties of related optoelectronic structures, has been report- 
ed’’*; however, these studies, by and large, are of an empiri- 
cal nature. Yet multicomponent systems play a key role in 
optoelectronic device structures, and their technology relies 
exclusively on epitaxial growth. Standard LPE provides 
only limit^ flexibility in controlling the microscopic growth 
velocity; thus, precise control of the composition of ternary 
and quartemary compound layers is still a problem. 

In the present study it is ^own that the same phenom- 
ena which dominate the electroepitaxial growth of binary 
compounds are also of fundamental importance in the 
growth of ternary and quartemary semiconductor com- 
pounds. Thus, a theoretical electroepitaxial growth model of 
multicomponent systems is developed based on electromi- 
gration in the solution and the Peltier effect at the growth 
interface. Quantitative relationships are derived for the de- 
pendence of the growth velocity and the composition of the 
layers on the current density and other growth parameters. 
The model is successfully used for the analysis of experimen- 
tal results on the electroepitaxial growth of Oa, Al^ As. 

II. MASS TRANSPORT MODEL 

A. ProcoMM and phMO diagram ratodonahlpa 

Electric current flow through the substrate (seed) in 
contact with a liquid solution of its components leads to a 
change d at the interface temperature due to the Peltier 
effect; AT^ u proportional to the current density and also 
depends on substrate characteristics such as the Peltier coef- 
ficient, thickness, and thermal conductance. In addition, the 
electric current flow leads to an electric field £ in the solu- 
tion, which is also proportional to the current density; this 
field causes electromigration of the species in the solutions 
with a velocity //£, where /i is the effective mobility. In gal- 
lium liquid solutions, nearly all solutes migrate under an 
electric field toward the anode. If the substrate has a posi- 
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SOLUTE 

CONCENTRATION C, (Al) 

nc. I . Schcmaiic rcproenution of the ptnincni pan of the phaac diagram 
of a (a) binary and a W ternary lyttcm and iranuiiont involved in eleciroe 
pitaiial growth. and x, arc the wlidui lines (mc teat). 


five polirity, electromigrition provides an effective means of 
solute transpon to the substrate leading to supersaturation 
and continued growth. It is thus apparent that electric cur- 
rent affects two key factors in epitaxial growth, i.e., the inter- 
face temperature and solute transport. 

The above two factors have been modeled’ quantita- 
tively in the electroepitaxial growth of binary compounds 
utilizing the widely accepted treatment of isothermal trans- 
port in LPE."-'^ Considering the binary phase diagram in 
Fig. 1(a), electroepitaxial growth is approximated as a tran- 
sition from point A„ (which is defined by the original equilib- 
rium temperature Tq and solute concentration C**) to point 
A, on the liquidus which is uniquely determined by the new 
interface temperature T, = To A- AT^ (AT^ is negative 
when T, < To) and defines a new solute concentration C ' at 
the interface. The solid composition of a binary compound is 
known, by definition, and thus the growth velocity can be 
obtained by solving the solute transport equations which in 
electroepitaxial growth are modified by including a current- 
dependent electromigration transport term. The basic as- 
sumption made here is that growth proceeds under thermal 
quasieqilibrium, i.e., the transition Aq -*A , follows the Uqui- 
dus line; the validity of this auumption in electroepitaxial 
growth has been experimentally verified.’ 

LPE of multicomponent systems (e.g., ternary and 
quartemary compounds) is more complex than that of bina- 
ry ones. Consider, for example, a ternary compound 
Ga, _ , Al. As in equilibrium with the (Ga-rich) Ga-AI-As 


solution. For a given temperature (e.g., Tq) the concentra- 
tion of the two solutes As and Al is not fixed (more degrees of 
freedom than in binary compounds), but can vary along the 
corresponding liquidus isotherm (To) u illustrated in Fig. 
1(b). The requirement that the solution be in equilibrium 
with a solid of a given composition Xo defines the equilibri- 
um solute concentration and Cj. and thus the starting 

point /<o- 

Ifsolution Ao (Fig. 1(b)] is thermally equilibrated at Tq 
and the interface temperature is decreased to T, , the growth 
point A , (in contrast to the case of a binary compound) is not 
uniquely defined but lies somewhere on the new liquidus 
isotherm T , , depending on the composition x, of the grow- 
ing solid. In order to determine the growth path Ao-*A, 
and, thus, AC,,ACj, and the new solid composition x, , an 
additional quantity must be introduced such as the differen- 
tial growth path function, ym dC] /AC, . It is apparent that 
the ratio ACj /AC, (at the growth interface) will depend on 
the rates at which the solutes are transported to the growth 
interface. In thermal LPE, these rates are determined by the 
diffusion proceu. In electroepitaxy they depend also on the 
electromigration properties of the two solutes. 

In the present treatment the differential growth path 
functions are obtained through a consistent treatment of the 
solute transport (diffusion and electromigration) to the in- 
terface and the phase diagram (liquidus and solidus iso- 
therms) relationships. As in the case of the binary system,’ 
here also a thermal quasiequilibrium approximation at the 
growth interface will be used. 

Since the Peltier effect-induced changes of the interface 
temperature are small” (of the order of I *C) the corre- 
sponding compositional changes of the solution are also 
small (i.e., |dC, |<C{*). Accordingly, the q^asiequilibrium 
approximation of the growth transition A' —A, can be ex- 
pressed in a differential form: 

dr=^dc, 4-^dc,-h...+ -^dc,., (I) 
or 

I 

where a, m dT/dC, lc, - cr .r, •** •*** liquidus slope compo- 
nents dC, is the change in concentration of a given solute 
component at the interface, n is the number of elements in 
the compound, and in the liquid phase the solvent concentra- 
tion is C and the solute concentrations are C, , where 

i B 1 n —I. It is apparent that the change of the solvent 

concentration it determined by the changes of the solute 
concentrations, i.e., — Xrj| dC,.Thus,informulat- 
ing mats transport in the solution only the solute compo- 
nents need to be taken into account, or only n —I elements. 


■. Sohita transport 

Assuming a dilute solution, the transport of each solute 
component can be treated independently. Thus, the trans- 
port equation for thermal LPE ” modified to include the 
electromigration term is expressed as 
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_ d»c, dr., dc, dc, 

(3) 

where i it the ditunce from the advancing growth inteifaoe, 
D, tt the diffusion contlant.v it the growth velocity, /I it the 
mobility, and r it the time. The sign oT/<,£dependt on the 
direction of current flow and it taken positive when the so- 
lute species electromigratet toward the growth interface (for 
convenience, the convention used here for the sign of /iff it 


opposite to that used in Ref S). 

The boundary conditions for solute concentration 
as follows: 

are 

at/>0 



C.-C? 

for all X, 

(4a) 

atr>0 



C.-C? 

for 2 - 00 (absence of convection). 

(4b) 

oratr>0 



c,-c? 

for i>6 (presence of convection). 

(4c) 


(S it the solute boundary layer thickness), and at r > 0 and 
r-0, 

C. - C;, (4d) 

(superscript I denotes values at the interface). 

Flux conservation at the interface requires that 

^C. 


+(!>+ /i,ff)c;, 

dz jfo 


(5) 


where C | is the concentration in the solid. 

Neglecting the term v(^C, /dx) in Eq. (3) as very small in 
comparison with the electromigration and the diffusion 
terms,’ the solution of Eq. (3) for conditions defined by Eqs. 
(4) is (see Ref. S) 


/D \'/J 

or i< • 0 \ irr / 


(6) 


where/ (ff,d,r ) is similar to the function /« (ff,d,r ) defined in 
Ref. 5. 

In the low electric field approximation 
[/i, E<2{D, // )'^’ ] which is valid for metallic solutions such 
as those presently considered, the function/ reduces to 

/ > I in the absence of convection (7a) 

/ ■ (irf),r)'^Vd in the presence of convection. (7b) 

From Eqs. (S) and (6) one obtains A C , , and substituting 
the resulting expression for AC, into Eq. (2), one obtains for 
the growth «'elocity 




obtained, and they will be discussed below in ooqjunction 
with compositional changes in the solid during 
electroepitaxy. 

C. (3rowth velocity 

It should be emphasised, at this point that the unique 
characteristics of electroepitaxy permit a significant simpli- 
cation of the crystal growth pr^lem, i.e., the growth veloc- 
ity can be treated independent of compositional changes in 
the solid. 

In contrast to thermal LPE, electroepitaxy allows the 
growth of thick layers (e.g., — 100 /im) without significant 
change in interface temperature. As pointed out earlier, 
AT^, associated with the Peltier effect at the growth inter- 
face, is typically of the order of I *C. The coresponding com- 
positional changes in the liquid and in the solid are also 
small, and AC] and AC] are much smaller than the equilib- 
rium values Cf and CJ. Thus, in Eq. (8), the concentration 
can be substituted by the equilibrium concentrations C f and 
Cf, which correspond to the initial temperature Tq. In con- 
trast, in thermal LPE, the overall temperature decrease dur- 
ing growth is often as high u 13-20 *C. Under such condi- 
tions the treatment of the growth velocity independent of 
compositional changa in the solid constitutes a very rough 
approximation. Nevertheleu. such an approximation re- 
duces the extent of numerical calculations and is often used 
in thermal LPE.” 

In the low electric field approximation, the electroepi- 
taxial growth velocity, in the absence of convection is 

47-0 




(«) 


a,(cr»-c?) y' 

If (A)’'* / 

(Av-» j • 

and in the presence of convection it is 

It is apparent that for n^l, Eqs. (9a) and (9b) become 
identical to those derived for the electroepitaxial growth ve- 
locity of binary systems’; if, in addition, /i, v 0, they, then 
reduce to the standard LPE expressions for diffusion<on- 
trolled growth of a binary compound. " ” It is of interest to 
note that for/i, ~ 0 or ff m 0 (no electromigration), and 
AT^ m AT, Eqs. (9a) and (9b) become 



a,(c;-c;)\-' 

V. 1 

f t' 

a,(C--C?)\ 

(«)•'> vf 

1 (AV'*/ ) 

-(rr)''*' 

l.f 

(A)*'* / 


(lOa) 


•(■t.sSXiNg;?)". 

In general since C,' and C* in Eq. (8) aic not known, a 
complete treatment of the problem must involve the growih 
path function ^ /AC, (for components k and / ) and 

the solidus relationships. Expressions for y,,, can be readily 


and 


They represent generalized expressions for standard thermal 
LPE growth velocity of multicomponent systems They can 
be used in conjunction with step-cooling thermal growth of 
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no. 2. Ekcmmiiiration ocMnbutwn to irowth velocity of 0«, . .AI.Ai 
if per unit currtnl demity celcuUled from E^. (9e), ukin| D^, " 
mix lO’ cm Viand - 1.1 X 10 *(cm/iKA/cm’) . 


ternary or quartemary compounds provided that d T is rea- 
sonably low (of the order of I *C). 

The first terms in Eqs. (9a) and (9b) define the Peltier 
contribution Vf to the growth velocity, while the second 
terms define the contribution from eiectromigration Vg. 
Both terms are proportional to the electric current density J 
(through d r, and E, respectively, as pointed out above). In 
the absence of convection, the Peltier contribution Vj- de- 
creases with time, while the eiectromigration term Vf is inde- 
pendent of time. Thus, the Peltier contribution is of signifi- 
cance during the initial suges of growth. When prolonged 
electroepitaxial growth is performed on thin substrates (of 
the order of (0. 1 mm), the electroepitaxial growth velocity is 
dominated by the eiectromigration term Vg . The presence of 
convection in the solution enhances the contribution from 
the Peltier effect, which is proportional to I/d. 



FIG 3 The Peltier efSoct ooniribut*jfi toekctroepiUAiel growth veloo«y of 
Ga, , Al. At ut (lae leit). calculated from Eq (9a). taking 
■ 5 X 10 ' * cm Vi and - 3D^. 


D.Oag.«AI.Ai 

In order to calculate the electroepitaxial growth vdoc- 
ity in ternary or quartemary compounds, knowledge of such 
parameters as the diflbsion constants D„ the mobilities^., 
and the resistivity of the solution is required. These param- 
eters are not known even for the most commonly used com- 
pounds. For this reason, the discuuion of Eqs. (^) and (9b) 
will be limited primarily to the role of the phase diagram, 
using Oa, _ « Al^ As as an example.'* 

The eiectromigration groivth velocity Vg, normalised 
with respect to the electric current density, is given in Fig. 2 
as a function of the solid composition x for different growth 
temperatures. In these calculations the diffusion constant of 
Al is taken equal to the diffusion constant of As in Oa-rkh 
solutions; the mobilities of As and Al are assigned relative 
values. It is seen that the eiectromigration grosrth velocity 
decreases with increasing Al concentration and decreases 
svith decreasing temperature. 

It is also seen in Fig. 2 that for the auumed Al mobility, 
equal to that of As. the growth velocity of Oa, _ , Al. As is 
smaller than in the case of pure OaAs. This fact is primarily 
due to the decrease of the As solubility in the liquid when Al 
is present; it clearly indicates that the Al cannot be treated as 
an impurity even when the concentration ratio Cg,/Cg, in 
the solution is of the order of 10 ~ 

The contribution of the Peltier effect to the growth ve- 
locity [Eq. (9a)]. normalised with respect to time and Peltier 
cooling (i.e., Vft *^Vd T, ). is plotted in Fig. 3 as a function of 
x. Here alto there are marked differences between binary 
and ternary compounds. 

It should be noted that the Peltier-induced change of 
the interface temperature AT, it determined primarily by 
the properties of the substrate (OaAs in the present cast). 
The contribution of the solution to A T, seldom exceeds 
10%. In the present case. A T,, corresponding to a given 
current density, substrate thickness, and electron (hole) con- 
centration. can be determined from the experimental data 
given in Ref. 3. 

An unusual feature of electroepitaxy is that growth can 
be carried out even in the presence of Peltier heating at the 



FIG 4 Toulcl«ctra(piiaualgiowtliv«iactlyorO«.,Al„,A«vcfi«M 
irowth lime for difermi valuci of ihc Pcltm bcalini (J f, > 0) U Ihc inter- 
face The panmetm umd in cnkiilsciont wm m 

P., bSxIO ’ cm Va, and /<..£//■ 11 x 10 *(cm/»NA/cm') '. 
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interface (jr^ pocitivc). This feature is illustrated in Fig. 4 
where the growth velocity of Oa, Al« As, calculated from 

Eq. (9b), is plotted u a function of time. Such behavior can 
be rcalited for Oai _ . Al. As by utilizing thinp-type OaAs 
sub*trates and "plus” polarity on the substrate. 

It should be noted that the calculations shown in Fip. 2 
and 3 can be very readily extended to the evaluation of the 
electroepitaxial growth velocity ofOO) _ , Al^, As in the pres- 
ence of convection, where a boundary layer o( thicknew i is 
present. The required modification are apparent from Eqs. 
( 9 a) and ( 9 b); this case will not be discussed further. 


A/ ~C;-(l-»-/j,£/u)C, ’ 


(II) 


wherep « )or I in the absence and in the presence of convec- 
tion, respectively. 

By denoting one specific solute component by C , , one 
obtains from Eq. (2) 

x'a.n,. ( 12 ) 

« • I 

Utilizing the common notation ^ , C. , for the sol- 

id phase of the compound (for ternary compound p or x 
equals I) and introducing differential segregation coeffi- 
cients * ,* - Sx/3C, I c . c„r„ and * / ■ ^/^C, | c - a.r„» 
the change of solid composition due to small changes of so- 
lute concentrations JC, becomes 

Ax - £' * f AC . , Ay - 'f k .”AC . . (13) 

1*1 * • I 

From this expression, the definition of Xa,. and Eq. (1 1), one 
obtains 


Ax 




(I4a) 



FIG 5 The microicopic vtlim of the aluimnufii concemraiion i in 
Ga. , Al. At. the differential growth path function y,, . and the growth 
vfkxit) cakubtad vmut dioance The initial temperature T, • IIO*C. 
the Peltier-induced change in interface temperature AT, ^ - 1 *C. and 
current dcntity 7-10 A/cm * Other parameterb are the tame at in Fig 4 



FIG 6 Am/ A T vcnui solid compoaition a of Oa, . . Al, At calcubtcd from 
the pment model for ebctrocpiuibl growth (solid lines) and thermal 
growth hy step cooling (broken lines). For electroepiuiial growth the Pel- 
tier effect contribution to the average growth velocity was taken as 
Vf - 0 25e, Other parameters are £>«, - - Va«« 

"SxlO ' cm Vs, and ^ A. II X 10 * (cm/sXA/cm^) 


and 


Ay 




(14b) 


where 7 ^, [Eq. (II)] is uken at the initial equilibrium point 
[A ® in Fif . (lb)], where T ■ To, C JarC J®, and C, For 
the reference solute component C, • it is convenient to select 
the element characterixed by the largest differential segrega- 
tion coefficient ft, (e.g., A1 in the case of Ga, ^ Al„ As). 

Equations (14a) and (14b), together with Eqs. (9a). 
( 9 b), and ( 10 ), describe the composition, growth path, and 
growth velocity in electroepitaxial growth of multicompon- 
ent semiconductors. 

It is of interest to note that in the low mobility l*mit 
yt, Eq. ( 10) reduces to the simpler form. 



which describes the differential growth path function in 
thermal LPE (in most instances convection is negligible and 
p » p. This growth path function is independent of the 
growth velocity. It is always positive and in most LPE cases 
its value is determined by the composition of the solid; usual- 
ly (D,/D^ is a weak function of the relative concentra- 
tions in the solution. It is thus evident that in diffusion-limit- 
ed thermal LPE, control of the composition during growth is 
primarily limited to the control of the cooling rate AT(t). 

In electroepitaxy, on the other hand, the differential 
growth path functions depend on the solute electromigra- 
tion and the growth velocity [Eq. (11)]. The ratio of the elec- 
tromigration velocity to growth velocity is typically about 
10. and thus, the solute concentrations C4 and/or Cf also 
contribute to the value of >^ 4 /. Another unique feature of 
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elcctroq>ilaxy it that the differential growth path functioni 
can be varied by varying the conductivity type (n or p) and 
the thickneu of the substrate. The overall electroepitaaial 
growth velocity v is the sum of the Peltier effect and electro- 
migration contributions; the magnitude and the sign of Pel- 
tier effect contribution can be varied by varying the substrate 
parameters. Utilizing Eq^ (9a), (I I), and (14), the solid com- 
position X was calculated as a function of distance from the 
original interface using parameten corresponding to 
Oa, . , Al« As growth from Ga-rich solutions, and the re- 
sults are shown in Fig. S, together with corresponding values 
of the microscopif, growth velocity v and the differential 
growth path function . It is seen that electroepitaaial 
growth leads to a remarkable stabilization of the solid com- 
position Thus, in elect roepitaxy, it becomes feasible to ob- 
tain thick layers with very small compositional gradients in 
the growth direction. It should be noted that this finding is in 
contrast to the originally proposed advantages of cunent- 
controlled growth, whereby the composition of the solid (do- 
pant concentration) could be varied (concentration gradi- 
ents) by varying the current density.'^ 

The composition stabilizing effect illustrated in Fig. S 
indicates that an electroepitaxial layer grown from a solution 
equilibrated at temperatures Tq with a solid composition x^, 
will exhibit essentially a constant composition 
jt « jtq -f J T^{Ax/A T ),^, where A is determined by the 
Peltier effect at the interface (it is a linear function of the 
electric current density and depends on the thickness, con- 
ductivity, and conductivity type of the substrate), whereas 
{Ax/AT)j^ can be obtained from Eqs. (9), (I l), and (14). 
The values of {Ax/AT)j^ corresponding to electroepitaxial 
growth of Ga, . . Al. As in the absence of convection are 
given in Fig. 6 as a function of solid composition for two 
values of the initial temperature Tq. For comparison 
(Ax/AT)f^ in thermal LPE, obtained from Eq. (15), is also 
shown 

It is seen that the value of Ax/ A T and, thus, the magni- 
tude of the current -induced char;ge in solid composition in- 
creases with increasing temperature and also increases with 



FIG 7 Currmf-mduerd changr in Aluminum conccntrAtion in 
Ga. . Al. A» for fleet rorpitAiiAl growth in the prcAmct of convAction in 
the solution. cAlculAfcd from the prnmt (tor teil) aa a function of the 
huundir) lA>rr ihit.knr%k S Diffusion coniunu And mobilitiCA »r cthc AAmc 
A% in Fig 6 
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FIO S. Fini-dchvAtivc photovoltAgc ptAki of hAndgnp trAntitiociA in 
Oa, ,A1.Ai In ihc upper portion o^ th» Sgurr the pCAkt Arc Ahown in An 
eipAndcd enerf> acaIc for eWctroepitAiiAl Uyen frown with different cur- 
rent denAity 


decreasing x. Thus, the optimum conditions for observing 
compositional changes as a function of current density can 
be d^ned as Tq 880 *C and solid composition 0. 1 -0.2. At 

To * MX) *C, and for solid compositions of 0.3-0.4, Ax/A T 
is 10 “ * C “ or less The low magnitude otAx/A T ex- 
plains why in previous electroepitaxy experiments per- 
formed with low current density (< 30 A/cm *), the current- 
induced changes in the composition of Ga, . . Al^ As could 
not be consistently resolved.*^* 

It is also seen in Fig. 6 that not only the magnitude, but 
also the sign of Ax/ A T depends on the initial growth condi- 
tions. Thus, for relatively k«w growth temperature. Ax/ A T 
can be positive in a region of low aluminum concentration 
and negative for higher values of x. This behavior clearly 
illustrates the significance of the phase diagram relation- 
ships in current-induced compositional changes of 
Ga, . , Al. As. 

In Fig. the current -induced compositional changes of 

Ga, . . Al. As are presented for electroepitaxial growth in 
the presence of convection leading to a boundary layer thick- 
ness 8. It is seen that the value of Ax decreases w'ith decreas- 
ing 8 and decreasing growth temperature T^. Thus, when 
growth is performed at relatively low temperatures and un- 
der pronounced convection in the liquid, the effect of the 
electric curro'.: on solid composition can be below the limit 
of experimenul detection. 

In summary, it is seen from the present model that in 
electroepitaxy, the solid composition x can be increased or 
decreased since the sign A can be changed without im- 
pairing the growth proceu. For certain solute composition 
and solute mobility ranges, the sign of Ax/ A T can also be 
reversed This flexibility of changing the growth path in elec- 
troepiuxy by selecting the growth parameters, e.g., the con- 
ductivity type of the substrate, n or p, (which are of no conse- 
quence in the thermal LPE) renders electroepitaxx 
particularly promising for studying the growth process 
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III. EXPERIMENTAL 

Electrocpitaxial growth experiments were carried out 
with Ga-rich Ga-AI- As solutions in a standard LPE appara- 
tus (horizontal graphite boat) modified to permit passage of 
electric current through the solution-substrate interface.'* 
The effect of current density on the growth velocity and on 
the solid composition was studied in a growth configuration 
in which convectional flow in the solution was essentially 
eliminated." Preliminary experiments were also carried out 
in a configuration where convection was intentionally intro- 
duced by superimposing a horizontal temperature gradient 
across the solution. The layers grown in such configuration 
exhibited noticeable compositional inhomogeneities along 
the growth plane and were found unsuitable for a meaning- 
ful analysis in terms of the theoretical model. 

The initial growth conditions, i.e., temperature Tq and 
composition of solution were selected so that (a) high 

current-induced relative changes or aluminum content in 
the solid, dx/xo, could be obuined, (b) knowledge of the 
phase diagram was assured, and (c) sufficient amounts of A1 
and As were present in the solution for the growth of layers 
about 20 /im thick without noticeable depletion of Al and As 
in the solution. Accordingly, experiments were performed at 
To = 870 ‘C. from a Ga-rich solution containing 2 g of Ga, 
1.47 mg of Al, and equilibrated with GaAs (C 
c^l.6x 10 ■ \ 10" ^). Current densities up to 80 

A/cm ^ were employed. In order to keep the total current at 
a level not exceeding 25 A (to avoid Joule heating in various 
parts of the growth system"), layers of only 0.5 cm in diame- 
ter were grown. 

For the verification of the electroepitaxy model, very 
precise measurements of compositional changes are re- 
quired. In Ga, _ , Al, As, for small values of x, the energy 
gap increases linearly with x and thus changes of energy gap 
provide a sensitive measure ofdx. In the present study, de- 
rivative photovoltage spectroscopy " was employed to deter- 
mine the energy gap transitions. For the analysis of the tran- 
sitions h the vicinity of the energy gap, the photovoltage 
generated on a MOS structure or on a Schottky barrier can 



FIG ^ Encrg> gapofGa, .Al, As versus solid composition jc □ indicates 
thepresem results. A is from Ref 20. o is from Ref 21; and I is from Ref 19. 
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be used. The advantage of the derivative (wavelength modu- 
lated) technique in comparison to a standard photovoltage 
threshold measurement (as used, e.g., in Ref. 20) is that it 
permits a direct and a far more accurate determination of the 
energy gap changes. These advantages are apparent from 
Fig. 8 in which typical first-derivative spectra obtained with 
Schottky barriers are given for Ga, _ , Al, As with different 
compositions. As shown in Fig. 8, the resolution of this ap- 
proach is sufficiently accurate for the determination of the 
current-induced changes of the energy gap (and thus of the 
composition) in electroepitaxially grown layers. 

In addition to derivative photovoltage spectroscopy, 
electron microprobe and Auger analysis, combined with ion 
milling, were employed for the compositional study of the 
electrovi axial layers. 

The rc. alts obtained on the energy gap as a function of 
composition (determined by electron microprobe analysis) 
are shown in Fig. 9, together with corresponding literature 
data.^ ^' Note that dE^/dx, determined from Fig. 9, is 
1 .2 ± 0. 1 eV. This value is used for the precise determination 
of current-induced compositional changes in Ga, .,A1, As 
from the small changes in the energy position of the first- 
derivative photovoltage peak. The electron microprobe itself 
could not be used for such measurements since a typical er- 
ror in microporbe measurements of x was about 0.01. 

IV. RESULTS AND DISCUSSION 
A. Growth vtlocity 

The experimental results on the electroepitaxial growth 
velocity of Ga, _ , Al, As as a function of electric current 
density is shown in Fig. 10; current density was the only 
parameter varied in these experiments. All experimental 
points correspond to the same initial composition of the so- 
lution, the sarie temperature (Tq = 870 *C) of the sys- 
tem, and the same substrate characteristics, i.e., (100) orien- 
tation, 500 ^m thick, n-type. Si doped with electron 
concentration n = 1 X 10 " cm “ \ at room temperature. 
Peltier effect-induced change in interface temperature per 
unit current density J T^/J was found to be about 
— 5 X 10 ^ *C cm VA. The results of Fig. 10 correspond to 
the average growth velocity , defined as = d //, where 
d is the layer thickness and / is the growth time. The vari- 
ation of growth thickness in the growth plane (typically 
± 10%) is responsible for experimental errors, which are 
marked with vertical bars. 

For a quantitative analysis of the results of Fig. 10 in 
terms of the theoretical growth velocity expressions, it 
should be noted that 

= — \ vdt = 2vr + Vc , (16) 

t Jo 

where Vj and are given by the first and second term, 
respectively, on the right-hand side of Eq. (9a). The calculat- 
ed average growth velocity as a function of current density is 
represented by the dashed line in Fig. 10. A very good agree- 
ment between theory and experiment is obtained by taking 
the diff'usion and mobility constants of Al about three times 
greater than those of As. The diffusion constant of Al in the 
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flG. 10. The average elect roepitaiial growth velocity of Ga, . , Al. As as a 
function of current density Points indicate expenmenlal and lines were 
calculated from E<,. (16) (see test). 


Ga-rich solution is not available in the literature. However, 
it haa been suggested, in conjunction with thermal growth of 
Ga, . Al, As, that it should be greater than the diffusion 
constant of As.^^ Other parameters, and 

used in the growth velocity calculations, are 1.1 X 10 
((cm/s)/(A/cm^)J * and 5 X 10 cm Vs, respectively; they 
are essentially the same as those previously reported for bi- 
nary Ga-As solutions.^ 

In the context of the compositional analysis discussed 
below, it should be noted that the value offi^^E/v obtained 
from the results of Fig. 10 is 10 ± 2. It is also of interest to 
note that the Peltier effect contribution to the average 
growth velocity (i.e., 2i;r) of Ga, _ , Alj, As was about 25%. 
This finding is consistent with the general characteristics of 
electroepitaxial growth kinetics in the absence of convection 
in the solution, whereby growth is expected to be dominated 
by solute electromigration. 

B. Composition 

Electroepitaxially grown Ga, . ^ Al^ As layers were 
fo. d to be highly homogeneous along the growth direction. 
Electron microprobe compositional profiles taken on sur- 
faces cleaved perpendicular to the growth plane, derivative 
photovoltage measurements (combined with chemical etch- 
ing). and Auger microprofiles combined with argon-ion 
milling showed that macroscopic or microscopic composi- 
tional changes across the layer thickness did not exceed 
djr/jr.^ :^3%. It should be noted that similar results have 
been reported in previous electroepitaxial experiments with 
Ga,..Al,As." 

This remarkable stabilization of composition in the 
growth direction is in excellent agreement with the present 
theoretical model of electroepitaxy (see Fig. 6 and related 
discussion in Sec. Ill ), and it is readily explained consider- 
ing that the major factors affecting the solid composition. 


Le., interface temperature and growth velocity, are essential- 
ly constant during electroepitaxial growth carried out under 
constant electric current density. The interface temperature 
is determined by the Peltier effect J while the growth 
velocity is controlled by the electromigration term Vg , which 
is independent of time. 

Experimental results on the composition of electroepi- 
taxial growth of Ga, . « Al. As as a function of current densi- 
ty are shown in Fig. 1 1. The experimental points correspond 
to the layers which have been analyzed with respect to 
growth velocity in Fig. 10. It is seen that the A1 content in the 
electroepitaxial layers decreases with increasing current 
density, in agreement with the present theoretical model. 
The energy gap for a given layer was determined with deriva- 
tive photovoltage spectroscopy using the value of 
dEg /dx = 1.2 ± 0. 1 eV (see Fig. 9). This procedure may in- 
volve some uncertainty in the determination of the absolute 
value of x\ however, the composition changes Jx as a func- 
tion of current density are determined with an error compa- 
rable to the error of dE^ /dx, i.e., better than 10%. This error 
is much smaller than the sum of other experimental errors 
such as those associated with the preparation of the solution, 
the measurements of Tq , and the determination of the energy 
position of the derivative photovoltage peak. 

For a quantitative analysis of the results given in Fig. 

1 1, it is convenient to rewrite Eq. (14a) for Ga, . Al^ As in 
an explicit form: 

Ax^AT^(k^^k iYi \ )/(«! + Oj Yix ). (17) 

where the suffixes 1 and 2 refer to A1 and As, respectively. 
For the substrate employed here. | J \ per unit current den- 

sity is 5 X 10 C/A cm ^ Accordingly, Ax/ A , corre- 
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FIG 1 1 The energy gap and corresponding compositional changes of elcc- 
iroepitaMally grown Ga, . Al, As as a function of current density (lower 
honzontal wale) and of Peltier effect induced change in interface tempera- 
ture A (upper scale) 
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tpondini to the dotted line in Fi|. II, iaO.SS ^ 0.IS 1/*C. 
Utilizing the phase diagram relationships'* of Oa, _ . AI, As 
and the above value of Ax/AT^, one obtains from Eq. (16) 
Xii * 3-2 ± Ii which, according to Eq. (10), yields 
Ha»E/v - 8 ± 2, taking - 3 - 3. •» 

in the growth velocity analysis. Thus, it is evident that the 
current-induced compositional changes and the growth ve- 
locity characteristics are satisfactorily accounted for on the 
basis of the present model utilizing a consistent set of param- 
eters (mobilities and diffusion constants) and the phase dia- 
gram data. 

V. SUMMARY 

b 

A theoretical model of electroepitaxial growth of multi- 
component semiconductor compounds was developed. It is 
based on diffusion, as well as electromigration mass trans- 
port in the liquid, and on phase diagram relationships. Mass 
transport by electromigration is dominant in the absence of 
convection, whereas in the presence of convection the contri- 
bution of diffusion-driven transport is enhanced. The result- 
ing growth velocity is proportional to the current density. 

The composition of the solid is controlled by the Peltier 
effect-induced change of the interface temperature, by the 
mobility and diffusion constants of the solute components, 
and by the growth velocity. A unique feature of electroepi- 
taxy is that the composition of the solid remains constant 
even for prolonged crystal growth, provided the current den- 
sity is kept constant. On the other hand, for a given solution 
composition, the composition of the solid can be varied by 
varying the current density, the variation of the composition 
depends on the phase diagram, and for Ga| . Alj, As it is 
more pronounced at higher temperatures and for lower val- 
ues of x(e.g., 880 'CxsO. 1-0.2). 

Electroepitaxial growth experiments were carried out 
employing the Ga, . ^ Al, As system. The results on growth 
velocity and solid composition were analyzed on the basis of 
the model and were found to be quantitatively consistent 
with all aspects of the model. 
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Minority-carrier mobility in p-type QaAs 
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Theoretical calculations of electron mobility in p*type OaAs were carried out taking 
into consideration the screening effects and all miyo'' scattering processes. Calculated 
values of mobility are presented as a function of carrier concentration, compensation 
ratio, and temperature. The basic differences between minority-carrier mobility in p- 
type GaAs and electron mobility in n -type GaAs are pointed out. A practical procedure 
is also pV^nted for the evaluation of minority-carrier mobility from available electron- 
mobility data. 

PACS numbers; 72.10. — d, 72.20.Fr, 72.20.Jv 


Although minority-carrier mobility is recognized as a 
key electronic parameter for device application, theoretical 
treatments of minority-carrier mobility in GaAs and related 
compounds have thus far been very limited. As a result, it 
has been a common practice to assume that the electron mo- 
bility in p-type GaAs is equal to the electron mobility in n- 
type GaAs with the same free-carrier concentration. Such 
practice ignores the presence of additional scattering centers 
(heavy holes) in p-type material.* Furthermore, it neglects 
differences in screening effects and electron statistics in n- 
type and p-type material. 

The present mobility calculations utilize the variational 
method *•* of solving the Boltzmann equation without in- 
volving the relaxation-time approximation. This method has 
recently been used to evaluate electron mobility in n-type 
GaAs as a fimction of carrier concentration and compensa- 
tion ratio. Comparison of these calculations with experimen- 
tal values of electron mobility has made possible the determi- 
nation of the compensation ratio in n-GaAs. It has also been 
shown’ that the compenution ratios obtained from mobility 
measurements are in good agreement with those determined 
from free-carrier absorption measurements. In the present 
paper the method used in Ref. 7 is extended to the calcula- 
tion of the minority-carrier mobility in p-type GaAs. This 
extension involves some modifications and additional as- 
sumptions. First it is assumed that the minority-carrier stat- 
istics in the conduction band are given by a nondegenerate 
distribution function described by the crystal-lattice tem- 
perature. In the case of photoexcited electrons this auump- 
tion requires that the electrons become thermalized within 
the conduction band prior to recombining with holes or be- 
ing trapped at localized centers. For GaAs in the tempera- 
ture range presently considered, this assumption is ju^ed 
since the effective relaxation time is several orden of magni- 
tude smaller than the minority carrier lifetime. 

In addition to the major scattering mechanisms consid- 
ered in Ref. 7, namely, ionized impurity, optical phonon. 
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piezoelectric scattering, and acoustical phonon, electron 
scattering by heavy holes will be taken into account. Because 
of their large effective mass, heavy holes will be regarded as 
fixed scattering centers.' Thus, the overall concentration of 
ionized scattering centers is given by + N f 

+ N D , wherepA is the concentration of heavy holes and 
N andfifo are the concentrations of ionized acceptors 
and donors, respectively. Neglecting the contribution from 
light holes and assuming that the minority-carrier concen- 
tration is much smaller thanpA one obtains 
p* = /V7 - SS andthusfV^2^4-. 

Another modification of the approach of Ref. 7 centers 
about screening effects. For a parabolic heavy-hole band, the 
screening length Lq is given by the following expression’: 



5.80X 10” 


(mA/mo)»^7-'^F 

^0 


( 1 ) 


where is the heavy-hole effective mass, do >* Oie static 
dielectric constant, F.(i 7 ^) is the nth-order Fermi-dirac inte- 
gral, and fff = Ef Tis the reduced-hole Fermi energy. 

The reduced screening energy, a, which describes screening 
of the electron-ionized impurity interaction is given by 


o= 

2m,Z< 


( 2 ) 


where m, is the conduction-band electron effective mass. 
The screening energy of the electron-optical phonon interac' 
tions has the following form'; 



where e. is the high-frequency dielectric constant. The 
screening effects arising from the conduction-band electrons 
and the valence-band light holes are negligible as the concen- 
trations of these carriers are much smaller than the concen- 
tration of the heavy holes. 

Compensation is included in the present calculations by 
introducing the compensation ratio defined as 

^ N p /S f . The hole concentrations enten into the cal- 
culation through the reduced-hole Fermi energy ; the 
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FIG. 1. Room-temperature electron mobility as a function of carrier con- 
centration. Solid curvet represent minority -carrier mobility in^type GaAs 
for different compensation ratios. Dashed curve shows the elwtron mobil- 
ity in a-type GaAs. 



FIG. 3. Temperature dependence of the electron mobility inp-type (solid 
curves) and a-type (dashed curvet) GaAs for two different carrier concen- 
trations: 10'* and 10" cm' and ^ « 0. 


compensation ratio affects the total concentration of ionized 
impurities as follows: 


The values of GaAs parameters used in Ref. 7 have also been 
adopted in the present case. The hole Fermi energy £ jf’ has 
been determined from the hole concentration taking the fol- 
lowing effective masses for the heavy (m^^) and light holes 
(ffi/): = 0.54mo and m/ s 0.089mo- 

The minority-carrier mobility as a function of hole 
concentration is shown in Fig. I, together with the electron- 
mobiiity,/i^”\ results reported in Ref. 7. It is seen that the 
values for|/|f ^ are lower than those for/i^''^ in n-type material 
with the same carrier concentration. Furthermore»/i^^ de- 
creases more rapidly with carrier concentration than/i^"\ 
The calculated room-temperature mobility ratios/i^V/i^”^ as 
a function of carrier concentration is given in Fig. 2. It is seen 



FIG 2. Ratio of electron mobility in p-type GtAt, to the electron 
mobility in n-type GaAs./i^”'. ei a function of free-carher concentration for 


that for high concentrations this ratio is as low as 0.6, which 
clearly indicates that the previously reported procedures 
for determining/i^’ from values of corresponding to the 
same free-carrier concentration, are not reliable for carrier 
concentrations exceeding 10'* cm~'. 

The calculated temperature dependence of/i|f^ andfi^"^ 
is given in Fig. 3. For the low concentrations the difference 
between and/i^”^ is more significant at low temperatures. 
This behavior is apparently due to an increased contribution 
of the electron scattering by heavy holes relative to that by 
the lattice. For higher carrier concentrations the screening 
effects become of greater significance and lead to a tempera- 
ture dependence of entirely different from that of 
Actually, for suiRciently low temperatures;/^^ becomes 
greater than 

The above behavior of the minority-carrier mobility 
can be understood if one considers the fact that with an in- 
crease in the hole concentration, there is (a) an increase of 
the total concentration of the scattering centers and (b) an 
increase of the heavy-hole screening. At room temperature 
the contribution from (a) is more significant than from (b), 
and the mobility decreases with increasing concentration. At 
lower temperatures the mean electron energy in the conduc- 
tion band becomes much smaller than the screening energy 
of the valence-band holes. Accordingly, screening is strongly 
enhanced and, consequently, (b) dominates (a). Thus, in 
contrast to n-type material, in p-type GaAs there is a tem- 
perature and hole-concentration region in which the elec- 
tron mobility increases rapidly .with decreasing temperature 
(Fig. 3) and also with increasing carrier concentration. 

The effect of compensation on the room-temperature 
electron mobility in p-type GaAs has been found to be slight- 
ly weaker than in n-type material. The computed values of 

with the compensation ratio as a parameter are given in 
Fig. 1 for hole concentrations ranging from 10'* to 10'* cm*'. 
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A comparison of the present oa/i^ with those oa 
(Ref. 7) showed that at room temperature can be ap- 
proximated with/i^"^ as follows: 

(5) 

taking n » and — (I d^/(3 — 6^ where d. and 0^ 
are the compensation ratios in n- and p-type material, re- 
spectively. Thus, for example in uncompensated p-type 
OaAs with a hole concentration p« is equal to the electron 
mobility in ti-type OaAs with an electron concentration 

and with compensation ratio 0, =• 0.33. This ap- 
proximation is pafticularly valid for low-carrier-concentra- 
tion material where screening effects are negligible. Howev- 
er, even for the higher concentration considered in this study 
the error associated with this approximation does not exceed 
10 %. 

From a practical point of view it should be noted that 
the validity of the above procedure is not limited to GaAs. 
Expression (S) can be applied to any p-type material with 
OsAs-like energy-band structure and scattering properties. 
For instance, in the case of indium phosphide the electron 
mobility values reported for n-type material can be uti- 
lized for obtaining the electron mobility values in p-type 
material. 

In summary, it has been shown that the minority-carri- 
er mobility in p-type GaAs is significantly affected by the 
presence of heavy holes which act as additional scattering 
centers. The present analysis shows that only at high tem- 


peratures and at low carrier concentrations the minority- 
carrier mobility in p*type GaAs equals the electron mobility 
in A-type mate^. For higher concentrations the room-lbm- 
perature minority-carrier mobility in p-type material can be 
significantly lower than the electron mobility in n-type mate- 
rial. The minority-carrier mobility wu found to exhibit an 
anomalous dependence on temperature and on carrier con- 
centration. A practical procedure was presented which per- 
mits the estimation of the minority-carrier mobility using 
available data on electron mobility in n-type material. 

The authors are grateful to the National Aeronautics 
and Space Administration for financial support. 
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Derivative surface photovoitage spectroscopy; a new approach to the study of 
absorption in semiconductors: GaAs 
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Derivative surface photovoitage spectroscopy was achieved with wavelength modulation; it was 
applied to GaAs and permitted, in a single experiment, the determination of changes in the 
ab^rption coefficient over an energy range of 0.5 to 4.3 eV. Photoionization characteristics of 
deep levels were determined. All known critical-point transitions up to 4.3 eV were clearly 
obtained from second-derivative spectra. In addition, previously unresolved transitions were 
observed- at about 2.6 eV. Oscillatory photovoitage in high magnetic fields was observed, and it 
was used to identify the transitions in the vicinity of the energy gap. 

PACS numbers: 78.20. — e, 78.50. — w, 71.25.Tn, 72.40. -h w 


This communication reports on the absorption changes 
in GaAs, over a wide energy range, involving photoioniza- 
tion of deep levels, interface states, transitions in the vicinity 
of the energy gap, and higher-energy critical points. Simulta- 
neous determination of these transitions was made possible 
with a simple experimental approach combining wavelength 
modulation and surface photovoitage spectroscopy with 
strikingly enhanced sensitivity through utilization of the 
MOS structure characteristics. 

Surface photovoitage spectroscopy has been employed 
primarily for surface state studies.' Its inherently low sensi- 
tivity in all standard experimental configurations does not 
permit derivative measurements. Wavelength modulation of 
photovoltaic effects associated with p-n junctions and metal- 
semiconductor barriers has been reported for indirect-gap 
semiconductors.^ ' However, />-/i junction studies' have been 
limited strictly to the indirect-absorption edge. The metal- 
semiconductor photovoitage spectroscopy has been ex- 
tended to the vicinity of direct edges.' Extension of this ap- 
proach well above or below the energy gap is not possible 
because of the high recombination velocity of minority carri- 
ers at the metal interface and the appreciable photoinjection 
of carriers from the metal to the semiconductor. None of the 
above limitations is encountered in derivative surface 
photovoitage. 

GaAs n-type single crystals and epitaxial layers with 
room-temperature carrier concentration l.Ox 10 to 
3.0 X 10 cm " ^ and with mobility 6900 to 3000 cm ^ A^sec 
were employed in the present study. The surface photovoi- 
tage (illumination-induced change of the surface barrier) 
was measureo with respect to a transparent reference gold 
electrode, evaporated on an oxide layer (800 to 1000 A thick) 
grown anodicall/ on the GaAs. With this insulating layer a 
sensitivity of 50 n V was achieved in the surface photovoitage 
measurements (i.e., two to four orders of magnitude higher 
than previously reported) permitting the determination of 
relative changes in the surface barrier as small as about 
10 


**00 leave from Institute of Physics, Polish Academy of Sciences. Warsaw. 
Poland 


The illumination system consisted of a quartz-halogen 
light source and a double-prism monochromator with a vi- 
brating slit. The illumination intensity and amplitude of the 
wavelength modulation were experimentally selected at lev- 
els (small) not affecting the measured spectra; typically 
AA Mo was of the order of 10 " Precautions were taken to 
eliminate the interference of spurious signals (arising, for 
example, from light-intensity modulation). The measure- 
ments were carried out in the temperature range 95-300 K. 

A typical sub-band-gap derivative surface photovoitage 
spectrum obtained with a melt -grown GaAs is shown in Fig. 
1 . The three-peak structure was found to be characteristic of 
all melt-grown GaAs crystals studied, although the relative 
amplitude of the peaks was found to vary with the supplying 
source of the crystals. In high-purity LPE layers, only the 
intermediate energy peak (-^0.9 eV) was observed, and its 
amplitude was about two orders of magnitude smaller. Ac- 
cordingly, the observed transitions are associated with bulk 
levels rather than with interface states. No photoionization 
transitions related to interface states were observed. Thus, 
the surface states in GaAs-oxide interfaces previously re- 
ported on the basis of capacitance measurements' apparent- 
ly are not optically active. 

By extending the surface photovoitage treatment of sur- 
face states' to the transitions involving bulk levels (in the 
surface space-charge region) it is shown that the derivative 
surface photovoitage (dV/dhv) is directly proportional to 
the derivative of the photoionization cross section (da/dhv) 
and thus to the derivative of the absorption coefficient. 

The basic features of the deep-level absorption charac- 
teristics (Fig. 1) coincide with the extensively studied emis- 
sion characteristics employing photoluminescence.^ Thus, 
the sensitivity of the magnitude of the peaks to the crystal- 
growth technique and the doping level are identical in both 
cases; similarly, in both cases the position of the peaks is not 
sensitive to the nature of the donor dopant. Furthermore, the 
position of the minimum ( 1 .02 eV) and that of the zero value 
of the derivative (1.23 eV) in Fig. 1 are similar to the ob- 
served energy positions of photoluminescence maxima. 
Thus, the two high-energy peaks (£, | 2 ) correspond to the 

same deep levels, previously observed by luminescence and 
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FIG. 1. Fint-dchvativc ipectnim of tubbandpp lurfacc phoCovoltife of 
( 100) GaAs(fnelt-|fX)wn) with n * 6x 10‘*cfn " ’.The contribution of the 
individual transitions is shown by dashed lines. Parameters used in fitting 
procedure: E, « MO. 0.905. and 0.83 eV below conduction band edge; 
r « 40. 1 S. and 40 meV, respectively; the Bohr radius of the shallow accep- 
tor 0 , * 3.S X 10 ’ ^cm; the ionization energy of shallow acceptor E, * 24.7 
meV was adopted from Ref. 10. 


attributed to complexes involving gallium vacancies and do- 
nor impurities; the similarity of the absorption and emission 
energies is significant, as it proves directly that there is no 
measurable Stokes shif: associated with these levels, con- 
trary to the literature report on Ge-doped GaAs.* In view of 
the above, the spectrum ofFig. 1 can be analyzed on the basis 
of the quantum-defect model of an impurity-band photoioni- 
zation cross section,* modified by the introduction of the 
Loreiitzian broadening F, 

Th^ results calculated from this model are shown in 
Fig. 1 . It is seen th?t they art in good agreement with experi- 
ment; furthermore, the present treatment enables the separa- 
tion of the individual contributions of the three deep ievel?< 
and the accurate determination of the binding energies, £, . 



Photon Energy iV 

FIG. 3. Second-derivative surface photovoltage of (100) GaAs epiuxial 
layer (a « 10 cm " ') taken in magnetic fields at 3(X) K. The direction of 
the magnetic field and that of the light propagation were perpendicular to 
the surface. 


It should be pointed out that with the same set of fitting 
parameters (given in the caption of Fig. 1) it was possible to 
obtain good agreement with experiment for all samples stud- 
ied in spite of the fact that the relative magnitudes of the 
experimental peaks were quite different. It is thus concluded 
that the optic^ly active deep centers in n-type GaAs can be 
approximated in the quantum-defect model as acceptor cen- 
ters, without invoking configurational changes. 

In the region of band-gap transitions second-derivative 
photovoltage spectra can be directly obtained. A typical 
spectrum is shown in Fig. 2. The photogenerated transition 
in the vicinity of the energy gap £o* spin-orbit splitting 
Eq -1- Jq* higher-energy structures £| and £, -f- J |, 

are clearly resolved. The energy positions and temperature 
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FIG. 2. Sccond-derivitive surface pbocovoltage tpactn of (100) GaAs (a « 6X 10 co ’ M^10'*cm. at nxxD tcmpcranire); d as h e d line 293 K; solid 
line 95 K. The amplitude of the peaks for -f EJ, £|, and £, -f J| has been multiplied by a factor of 50, 20, 3, and 3, rsapactivdy. 






dependence of these critical>puint trsnsitiont ate in good 
•greement with thoM obtained by reflectance meuure- 
ments."-" It it important to point out that the half-width of 
photovoltage structures £q ** ^ '* compara- 

ble to that obtained at 4.2 K in high-resolution electrorcflec- 
tance spectra.'* 

The structure observed in Fig. 2, designated as £ ? , 
(2.S9 eV at 93 K) coincides in energy with the transitioD 
originally reported in Ref. 1 1 but not oonflrmed in any later 
study.'* Since this structure, unlike the othen in Fig. 2, is 
characterized by a minimum-mazimum sequence, it can be 
tentatively identiflud with the oi My critical point.'* 

In order to establish directly the contribution from 
Py — Py transitions to the structure in the vicinity of £^, 
derivative surface photovoltage measurements were carried 
out in high magnetic fields. The results obtained at room 
termperalure are shown in Fig. 3. It is seen that (a) the domi- 
nant structure in the vicinity of £q is essentially unaffected 
by magnetic fields; accordingly, this structure is not associ- 
ated with P, — Py transitions or with free-exciton transi- 
tions, but rather with residual impurities and/or defects. 
Thus, the structure at £q cannot be identified with the ener- 
gy gap. (b) For higher photon energies magnetic field leads 
to an oscillatory magnetopliotovoltage with the minima and 
maxima shifting toward higher energies with increasing 
magnetic field. The quantitative account of the involved 
Landau-levels transitions is complicated in GaAs by the 
complex nature of its valence band. However, from the con- 
vergence of the magnetophotovoltage oscillations at zero 
magnetic field the energy gap at 300 K of high-quality GaAs 
is precisely obtained as = 1.440 ± 0.002 eV. 

In view of the above, it is not surprising that this value is 
greater than those previously obtained from the derivative 
spectral in the vicinity of £© (1.420-1 .427 eV).‘* Transitions 
in the vicinity of £q have been extensively studied in GaAs, 
as they are important in understanding laser action. Howev- 
er, due to uncertainties in the identification of transitions at 
room temperature, the energy-gap values obtained by var- 
ious methods range from 1.38 to 1.443." 

The present results (Figs. 2 and 3) clearly show that 
derivative surface pholovoltage is highly sensitive to band- 
structure optical transitions. This sensitivity can be qualita- 
tively explained by the theory of surface photovoluge aris- 
ing from trapping of photogenerated holes by surface (inter- 
face) states." The spectral response of surface photovoltage 
is determined by aL^{\ -f aL^) where a is the absorption 
coefficient. Thus, for small values of the minority-carrier 
diffusion length!, (eg., 10 “ ’em) the surface photovoluge 
is essentially proportional to the absorption coefficient; as 
!, increases the sensitivity of photovoluge to absorption 


coefficient decreases. Consistent with this model the ampli- 
tude of high-energy peaks in the spectra (as those shown in 
Fig. 2) wu found to decrease (with respect to the low-energy 
peaks) with increasing However, the high-energy 
(£, -f J ,) structure could be resolved even for high-quality 
(3aAs epitaxial layers (£, and a!, > 1); this result 

can be qualiutively explained considering that for high val- 
ues of a (high-energy region) the holes are generated in the 
space-charge region where an electric field (neglected in 
treatment of Ref. I6)ashighu 10’-I0*V/cm opposes their 
diffusion; accordingly, !, becomes a much shorter field- 
contracted diffusion length." 

In summary, it wu shown that derivative surface pho- 
tovoluge spectroscopy coiutitutes a unique means for de- 
rivative abMiption spectroscopy studies at energies below 
and well above the energy gap of GaAs. This technique 
should be applicable to unipolar semiconductors in general. 

The authors are grateful to the National Aeronautics 
and Space Administration for financial support. The authors 
are also grateful to Dr. R.E. Enstrom and Dr. D. Richman of 
RCA Laboratories for providing the high-quality epitaxial 
material and the Francis Bitter National Magnet Laborato- 
ry for providing the high magnetic field facilities. 


■H.C. Gatw md J. Lagowiki. J. Vac. Sci. Technol. 10. 130(1971). 

*T. Nithino and Y. Hamakawa, Pbyt. Sutut Solidi B M, 343 (1972). 

T. Nithino. M. Takeda. and Y. Hamakawa. Suif. Sci. 97, 404 (1973). 

*H. Hategawa and H.L. Hannagel. 3. Eleciroebem Soc. 129. 713 (1976) 
'L.A. Chctler and O.Y. Robimon. 3. Vac. Sci. Technol. IS. 1525 (I97S). 
*C.L. Balcttra. 3. Ugowtki, and H.C. Oaiot. Surf. Sci. M, 457 (1977). 
'O.P. Peka and V.A. Brodovoi. Sov. Pbyt Scfflicond. 7, 1 100 (1974); S. 
Meti and W. Friu. Inti. Pbyt Conf. Scr B 99a. 66(1976) 

•E.W Williamt and A.M. While. Solid Suie Commun. *. 279 (1971). 
’H.B Bebb. Pbyt Rev. lU. 1116(1969) 

"N.O. Lipardi and A. Baldaretchi. Pwttdintt Elevtuih tnumai'1. Cv\f 
Phj/t. Stmicond., Wartaw (Polith Scieniific PuMithen, Wanaw, 1972), p. 
1009 

"D.L. Orcenway. Pbyt. Rev. Leii. 9 . 97 (1962) 

"D.D. Sell and S.E. Slokowtki. hvct*dinp Tenth Initmat’l. Conf. A)y>. 
Semkond., Cambridfe. Matt (U.S. Alomic Ener|y Committion, Oak 
Ridge, Tenn.. 1970). p. 417, I P. Waller, R R.L. Zucca, M.L Cohen, and 
Y.R Shen, Pbyt Rev Leu 24. 102 (1970). Y.R. Shen. Surf Sci. 97, 522 
(1973) 

"D.E. Atpnet and A.A. Sludna. Surf. Sci 97, 631 (1973). 

“B. Btu. in Semiconduclon and SemimetoU, ediled by R.K. Willardton 
and Albert C Beer (Academic. New York. 1972), p 315. 

"See for etample, I Camatiel, D. Auvergne and H. Malhicu, J. Appl. Pbyt 
46. 2613 (1975) and referencct quoted therein. 

"N L Dimiiriuk, V I Lyathenko, A K. Terethenko, and S.A. Spektor, 
Pbyt. '' it Solidi A 20. 53 (1973). 

"N J. harrick, Solid-Suu Electron 1. 234 (I960) 


9001 


J. Appl Pttys . Vol 50. No 7. Jiiy 1070 


Commurucabona 


9061 


WAVELENGTH MODULATED PHOTOCAPACITANCE SPECTROSCOPY 
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Abstract 


Derivative deep level spectroscopy was achieved with wavelength 
modulated photocapacitance employing MOS structures and Schottky 
barriers. The energy position and photoioniratlon characteristics 
of deep levels of melt-grown GaAs and the Cr level In high resis- 
tivity GaAs were determined. The advantages of this method over 
existing methods for deep level spectroscopy are discussed. 


Journal of Applied Physics, In press. 
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Capaclcanctt measuraiMnts h«v« r«c«ncly r«c«lv«d • great deal of 

(1-3) 

atcancloa aa a maana for datarmlnlng daap lavala la aamlconductora. 

Tha moat affactlve and tanaltlva of thaaa mathoda raly almoat axclualvaly 

(3 4) 

on maaauremanta of tharmal amlaaion of carrlara trapped by daap lavala. ' * ' 
Tha photolonlaatlon of daap trapa haa bean utlllaad la conjunction vlth 
capacitance maaeuramanta* however » with much laaa auccaaa due to amblgultlaa 
In tha axperlmantal and theoretical datarmlnatlon of phot Ionization croaa- 
aactlona. 

In tha preaant paper It la ahown that a reliable datarmlnatlon of 
photolonlzatloa characterlatlca of daap lavala can be achieved by a almpla 
approach combining high aanaltlvlty photocapacltanca maaauramanta with wave- 
length modulation of tha Incident light. Tha advantagaa of wavelength modu- 
lation for atudylng daap lavala have recently baan damonatratad in derivative 
aurf aca photovoltage apactroacopy aa applied to GaAa MOS atructuraa. It 
%d.ll be praaantly ahown that the derivative wavelength modulated photocapacl- 
tanca apactroacopy la applicable to MOS aa wall aa to atructuraa almplar 
than MOS, a.g., matal-aamlconductor Schottky barrlara and p-n junctlona. 

Tha baalc capacitance maaaurementa ware performed with a 1 MHz commercial 
capacitance mater in a differential mode. Tha changea of capacitance Induced 
by wavelength modulation of tha incident light (vlth fraquanciea ranging from 
1 to 25 Hz) ware detected ualng a lock-in amplifier, aynchronlzad with the 
frequency and phaae of the light modulation. In thla way the flrat derivative 

of tha capacitance with raapact to tha photon energy vaa obtained. Tha ayatam 

-4 

employed permitted tha detection of capacitance changea aamllar than 10 pF 

—7 

and of relative changea In the total capacitance aa amall aa 10 -10 The 

Illumination ayatem conalated of a tungaten light aource and a double prlam 


■onochroBBtor with b vlbrBtlng silt. Th# illuminBtion IntBnalty vbs axparl- 
BtntBlly SBlBctBd Bt bobII IbvbIs not Bffactlng tha aaaaurad apactra. Tha 

aaplltuda of tha wavalangth BodulaClon* AX/X* was typically of tha ordar 

•2 -3 

of 10 . In cartaln BaBsuraaants a raductlon of 4X/X to about 10 waa 

found nadaaaary In ordar to ravaal a flna atructura of photoioiizatlon transitions. 

Pracautlons wars takan to allalnata tha Intarfaranca of spurious signals 

(arising, for axampla, froB light Intanslty Bodulatlon). 

Darlvatlva photocapacltanca aaasuraBants wars parfomad on n-typa GaAs 
utilizing MOS structuras and Schottky barrlars. Tha MOS structuras vara aada 
by avaporatlng a tranaparant gold alactroda on an oxlda layar (about 1000 A 
thick) grown anodlcally on CaAs. Schottky barrlars (MS) vara aada by avapor- 
atlng a tranaparant gold flln on GaAs. 

Typical subbandgap darlvatlva photocapacltanca spactra (dC/dhv) obtalnad 
with MOS and MS structuras ara shown In tha uppar part of Fig. 1 (the lotfar 
part of Fig. 1 will ba dlacussad balov); thasa structuras vara praparad on 
adjacant araas of tha sana aalt-grown GaAs wafar (rooB tanparatura carrlar con- 

1C 2 

cantratlon 6 x 10 cm and mobility 4.500 ca /Vsac). Both spactra exhibit 
tha sama vall-daflned thraa-pask atructura. Tha slallarlty of tha derivative 
photocapacitance spectra of MOS structure and Schottky barrier Is significant, 
as It constitutes a direct proof of the fact that the transitions involved are 
associated with bulk levels rather than with Interface states. Furthermore, 
tha positive sign of tha corresponding subbandgap photocapacltanca Identifies 
thasa transitions as photo-axci tat ions of electrons froa deep levels Into 
tha conduction band. 

Tha thraa-paak structure observed in derivative photocapactisnea was 
found previously to ba characteristic of all melt-grown n-cypa GaAs. 
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It is also of intortac to noca chat tha potlclon of tha 1.0 aV alniauo and 
chat of tha taro yalua of tha capaclcanca darivatlva (1.25 aV) in Pig. 1 
corraapond to tha paak anarglaa obaarvad In lualnaacanca. Thua* tha high 
anargy paaka, and £ 2 * corraapond to tha aama daap lavala pravloualy 
obaarvad by lumlnaacanca and attrlbutad to conplaxaa Involving gallium 
vacanclaa and donor Inpurlclaa. It ahould ba notad, howavar, that high 

a » 

raaoluclon maaourananta of darivatlva photocapacltanca hava ravaalad a flna 
thraa-paak atructura of paak £2 (aaa Inaart of Fig. 1) which hat not ba^ 
raaolvad In any pravloua atudlaa. Tha lowaac anargy paak, (dua to a 
laval at 0.8 aV balow tha conduction band adga) can ba attrlbutad to a laval 
commonly Idantlflad with oxygan, although chara la no unanblguoua proof of 
chit Idanclflcatlon. 

Tha aama thraa-paak atructura aa In darivatlva photocapacltanca (Fig. 1) 
waa obaarvad In darivatlva aurfaca photovolcaga maaturamanta parfomad on 
tha aama MOS davlcaa (Fig. 21 . Hovavar, tha darivatlva photovolcaga apactrum 
obtalnad with a Schottky barrlar (Fig. 2) doaa not axhlbic thla atructura. In 
chia caaa tha apactnm la domlnatad by Injacclon of alactront (cbraahold at 
hv s 0.8 aV) and of holaa (hv^ 1 aV) from tha natal Into tha aamlconducCor 
which maaka tha photolonizatlon tranaltlona from daap lavala. Tha advantaga 
of photocapacltanca maaauramanta with Schottky barrlara la thua apparent. 

The quantitative traataant of derivative photocapacltanca la almpllflad 
algniflcantly by tha fact that the Intarfaca charge la not modulated by aubbar'*- 
gap Ulumlnaclon (aa concluded from tha almllarlty between raaulta obtained 
with Schottky barrlara and MOS atructuraa) . By aztandlng aarllar photocapaci- 
tanca and photovolcaga traatmanta^^^*^^^ to tranaltlona Involving photolonizatlon 
of bulk lavala (In tha depletion apace charge region) It la readily abown chat 
anall algnal darivatlva photocapacltanca, dC/dhv, la directly proportional 
to tha derivative of tha photolonizatlon croaa-aactlon, do/dhv, 

i£_ . 1 do 

dhv 2 ''aff dhv 


( 1 ) 


- 5 - 

vh«rt X is ths Incldsnt photon flux, n^ la tbs concantratlon of occuplsd 

traps bslng probed, and 1* the net doping concentration, ■ C for 

2 

a Schettky barrier and ■ C[1-(C/C^^) ] for an HOS structure; Is 
the oxide capacitance. With the exception of very low teoperatures, n^ changes 
only slightly under low Intensity subbandgap Illumination and thus for the 
states located below the Fermi level, n^ practically equals the concentration 
of traps, N^. Similarly, the derivative photovoltege, dV/dhv, Is given by 

■ - V I A — ^ (2) 

The coefficient A In eqs. (1) and (2) depends on the thermal generation and 
on the recombination rates which characterize the Interaction of deep levels 
with conduction and valence bands. In general it Is also dependent on the 
rate of optical generation. 

For sufficiently high traperatures and low intensity of Incident illumina- 
tion, whereby a transient response of deep levels can be approximated by a 

2 2 1/2 

relaxation time, T, the coefficient A has a simple form A ” <t/ (1 u T ) >, 

where u) Is the angular frequency of wavelength modulation. A in Eqa. (1) and 
(2) represents an average value since t depends on the location of the traps 
In the space charge region. The approximation of small illumination Intensity 
was verified in the present experiments by the observed linear dependence 
of the signal on illumination intensity. The coefficient A was also found 
to be Inversely proportional to the frequency of light modulation, l.e. , 
wT>>l. At lower temperatures (<80*K) the derivative photocapacitance (and 
surface photovoltage) decreased significantly because under illumination the 
occupation of states decreases rapidly due to the increase of the relaxation 
time. No measurable derivative signals could be detected et temperature 
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50*K. The similarity of Eqs. (1) and (2) shows that the derivative photocapacl- 
tance spectra of Fig. 1 can be quantitatively analysed on the basis of the quan- 
tum defect photolonlzatlon model employed for the analysis of the derivative 
photovoltage spectra. Such a procedure enables the separation of the 
Individual contributions of the deep levels and the accurate determination of 
their binding energies* E^. The presently obtained values for peaks and E^ 
are 1.1 and 0.83 eV, respectively, l.e. , Identical to the values obtained from 
surface photovoltage measurements. The high resolution measurements performed 
with MOS structures showed, however, that the E^ peak in all Investigated 
samples exhibited a fine structure which cannot be accounted for by existing 
photlonlzatlon models. 

Since the derivative photocapacitance Is directly proportional to the 
derivative of the photolonlzatlon cross-section (Eq.. 1) It must also be pro- 
portional to the derivative of the absorption coefficient. A typical sub- 

( 12 ) 

bandgap absorption spectrum of GaAs Is shown In the lower part of Fig. 1. 

It Is seen that the derivative maxima indeed coincide with the highest absorp- 
tion slopes. Although the absorption spectrum corresponds to a lover tempera- 
ture (21*K), its comparison with the derivative photocapacitance spectra is 

valid In view of the weak temperature dependence of Che energy position of 

( 12 ) 

the derl'^atlve peaks and of Che corresponding thresholds. 

According to equations (1) end (2) , the ratio of Che derivative photo- 
capacitance to the photovoltage signal for MOS structures, (dC/dhv)/(dV/dhv) , 

Is equal to 1/2 ^ depletion layer, a decrease of the surface 

barrier height, V, corresponds to an Increase In capacitance C and thus 
(for C/C < 1). Accordingly, Che ratio (dC/dhv)/(dV/dhv) should Increase 

OX 

with decreasing the surface brv'ier. This marked difference between surface 
photcvoltage and photocapacitance has been verified by experiments in which 
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the surface barrier height was decreased by low Intensity background Illumina- 
tion with hv > E . As expected » both derivative spectra decreased In magnitude 

O 

due to the decrease In the occupation of the deep traps by electrons (capture 

of photoexclted holes). However, the ratio (dC/dhv)/(dV/dhv) Increased under 

backg'*oui^d Illumination; furthermore, the derivative photocapacitance spectrum 

could still be clearly observed when the photovoltage signal essentially 

vanished. This result Indicates again the advantage of photocapacitance 

measurements under conditions of low surface barrier and/or when Interaction 

of deep centers with minority carriers Is significant. 

In view of its high sensitivity to deep levels, derivative photocapacl- 

tance spectroscopy was applied to high resistivity Cr-doped GaAs. Typical 

derivative spectra obtained at room temperature and at 100 R are shown In 

the upper part of Fig. 3. In the lower part of Fig. 3 typical absorption 
n 3) 

spectra for similar material are shown. The zero crossing point of the 
derivative photocapacitance spectra at 0.9 eV clearly coincide with the absorp- 
tion peak observed at the low temperature. The shape of the derivative spectrum 
Is consistent with a transition between two discrete states, rather than between 

a localized state and the conduction band. This conclusion Is in agreement 

(13) 

with results of previous investigations. 

In summary. It was shown that derivative photocapacitance spectroscopy 
permits the detection of small changes of optical absorption associated with 
photoionization of deep levels In semiconductors. i technique can be used 
wit.. M0s> structures, Schottky barriers, and undoubtedly with p-n Junctions. 

It permits the detection of deep bulk levels in GaAs-MOS structures In the 
presence of a Mgh density GaAs-oxlde Interface state. According to a 
realistic estlu-.tlon, deep levels with a concentration two orders of magnitude 
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smaller chan shallow impurity doping can be detected. 
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FIGURE 'V,?TI0NS 

D«rlvativ« photocapacitance apectxa of nelt-grotm GaAs 
15 -3 

(n ■ 6 X 10 cm ) for MOS and MS atructures (upper portion 

»2 

of Fig.) obtained with wavelength modulation AX/X of about 10 . 

k 

The Insert shows a fine structure of peak E2 as revealed by 

high resolution measurements on MOS structures employing AX/X 

.3 

of about 10 . In lower portion of Fig. a typical absorption 

spectrum Is given (Ref. 12). 

Derivative surface photovoltage spectra for tbe same samples 

_2 

as In Fig. 1. The amplitude of wavelength modulation AX/X <• 10 . 

The MOS spectrum Is multiplied by a factor of 10. 

Derivative photocapacitance (upper portion of Fig.) and 

absorption spectra (see Ref. 13) of Cr-doped GaAs. The 

_2 

amplitude of wavelength modulation AX/X was about 10 . 
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. k ABSTRACT 

Interface states and bulk GaAs energy levels were siijiultaneously investi- 
gated in GaAs MOS structures prepared by anodic oxidation. These two types 
of energy levels were successfully distinguished by carrying out a comparative 
analysis of deep level transient capacitance spectra of the MOS structures 

and MS structures prepared on the same samples of epitaxially grown GaAs 

15 —3 

(n “ 10 cm ). The identification and study of the Interface states and 

bulk levels was also pursued by investigating the transient capacitance 

spectra as a function of the filling pulse magnitude. It was found that in 

the GaAs-anodlc oxide Interface there are states present with a discrete 

energy (0.65 eV below the conduction band) rather than with a continuous 

energy distribution as previously reported. The value of the capture cross- 

section of the Interface states was found to be 10 to 10 ^^cm^, which is 

~8 -9 2 

far more realistic than the extremely large values of 10 to 10 cm reported 
on the basis of conductance measurements. 



Proceedings of Seventh Annual Conference on Physics of Compound Semiconductor 
Interfaces, Estes Park, Colo., January 1980. 


I. INTRODUCTION 


The rapidly Increasing interest in the growth and properties of native 
oxide layers on GaAs stems primarily from the Importance of these layers 
In electronic device applications such as MOS devices and Integrated circuits. 
The electronic properties of GaAs-oxide Interfaces have so far been Investi- 
gated utilizing primarily C-V measurements. Commonly, GaAs MOS structures 
exhibit carrier-injection type hysteresis and large frequency dispersion which 
are considered to be due to the high density Interface states and/or to traps 
in the oxide. ^ * These itxtcrface characLerlstlcs are undesirable for active 

device applications and thus, the importance of the quantitative charac- 
terization of the Interface states and of the understanding of their origin 
is quite apparent. Precise determination of the Interface state parameters 
from the conventional C-V or -conductance methods Involves ambiguities related 
to the very broad time-constant dispersion of GaAs MOS structures. In 
the case of silicon MOS structures the time constant dispersion Is commonly 
attributed to surface potential fluctuations. It Is conceivable that 
in GaAs MOS structures also such fluctuations Influence significantly the 
results of C-V and conductance measurements. 

(8 9 ) 

In this investigation deep-level transient spectroscopy (DLTS)' * * was 
applied to study the parameters of GaAs-native oxide Interface states. Unlike 
C-V and conductance methods, transient capacitance measurements are Independent 
of the surface potential fluctuations. Recently DLTS has been applied to 
the Investigation of interface states in GaAs MOS structures having an 
oxide layer prepared by plasma oxidation; in* the present study the oxide layer 
was obtained by anodic oxidation in a solution of glycol and water. The experi- 
mental approach permitted the investigation of Interface states as well as 


residual bulk levels. 


U. EXPERIMENTAL 


GaAs LPE layers, n*type 45 yn thick, grown on Cr»doped semi-insulating 

substrates with an (100) orientation, were used In this study. At room 

IS -3 

temperature the layers exhibited an electron concentration of 1.1 x 10 cm 

2 

and a Hall mobility of 6900 cm /V-sec. A similar value of net impurity 

concentration ■ n) was obtained from C-V measurements. 

« 

Anodic oxidation of CaAs was carried out in two different electrolytes: 

(a) a mixture of propylene glyool and 3% aqueous solution of tartaric acid 

(liif£fcVo.il by ird^Ofl to pH of about 6.2) in a volume rur.io and 

(b) a mixture of diethylune glycol + 5% by weight water 0.5% by weight 

Oxidation was initiated under a constant current density of 
2 2 

1 mA/cm for electrolyte (a) and 0.1 mA/cm for electrolyte (b), until a pre- 
determined voltage was reached; than a constant voltage was maintained until 
the current decayed to about 10% of the initial value. During oxidation the 

GaAs samples were Illuminated. with 'a collimated beam of light from a tungsten 

( 12 ) 

lamp In order to generate holes required for the anodic reaction. * Prior 
to anodic oxidation the samples were etched in concentrated HCl to remove the 
as-grown oxide. The thickness of the anodic oxides employed was typically 
between 2000 and 2500 A. The HOS structure was prepared by vacuum deposition 
of gold through a metal mask. The MS structures were prepared from the same 
GaAs samples used for the MOS structures after dissolving the anodic oxide. 

In the C-V and dC/dV measurements an IMHs commercial capacitance meter 
was used. The DLTS measurements were performed in a standard experimental 
arrangement. The rate window was implemented with a double boxcar 
averager. The sampling times were selected so that t 2 * 2t^; t^ was varied 
from 0.5 to 25 ms. 
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IXX. RESULTS AND DISCUSSION 
3. 1 Capacltanf-Voltaftc M«asur«mgt >f 

Typical C-V character la tica of MOS atructuraa prepared by employing the 
propylene glycol electrolyte are ahown in Fig. 1. Very ainilar reaulta were 
obtained with MOS acructurea prepared by uaing diethylene glycol, although 
the resiativity of the oxidea prepared in the former electrolyte waa higher 
thnn the rcaiati.vlty of thoae prepared in the latter. Aa aeen in Fig. 1, 
r</.s MOS atructuros x^vcnl ccvcrcl distinct fcaturea conaistent with previous 
f im'inf.c: (1) CiiracitoncH-voltri;* chnr;jr*.. eristics f.ho.? certicr-injcction 

type hysteresis; (2) the 1~M1U capacitance of the MOS structures under forward 
bias is much lower than the oxide capacitance; (3) at high reverse bias the 
samples show a deep*depletion type characteristic. 

The deep>depletion type behavior is especially important for DLTS measure- 
ments and its validity was tested in more detail. Measurements of the total 
capacitance, C, and the derivative of the capacitance, dC/dV, (dV ■ 10 to 100 mV, 
frequency of voltage modulation 10 to 100 Hz) were carried out in the temperature 
range of 80 to 500 K. Provided that the MOS structure is in the deep depletion 
state these measurements allow the determination of the net doping concentration, 
Np - N^, according to the relation: 

( 1 ) 

where q is the electronic charge, is the dielectric permittivity of GaAs, 
and A is the area of the field electrode. 

It was established that the net concentration of ionized impurities, as 
determined from room temperature C-V results under deep-depletion conditions, 
were in excellent agreement with values obtained from Hall-effect, measurements. 

The doping concentration determined from eq. (1) changed only slightly with 
temperature and was consistent with that obtained from measurements with MS 
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atruccure>« Thesa results clearly show Chat in Che teaipersture range presently 
studied the MOS structures under high reverse bias were indeed in the deep- 
depletion state. 

3.2 PITS Meaaurcwents 

PITS measurements on MOS structures were carried out under deep-depletion 

conditions. Typical results are shown in Fig. 2. At low filling pulse voltages 

(up to 20 volts) only a high temperature (cbout 390 K) peak is observed. By 

iiicrcssin.'j, th---. fjll‘;<g pulse voltage rat additional peak is revealed at a lov;«r 

tc“.;i.'>ri;Mu‘c (aboci 3C0 K)* In gener.- X incrf.'.sin^ the- f^lliti;^ pulse VoJtap.u 

(8 9) 

should enhance the contribution from traps located close to the lnter£ace;_ * ^ 
thus, the lower temperature peak observed only for filling pulses exceeding 
20 volts can be attributed to Interface states. Consistent with this interpreta- 
tion, the lower temperature peak was not observed in MS structures prepared 
with the same samples; the high temperature peak (390 K) was clearly observed 
which indicates that this peak is associated with bulk traps. 

Thermal activation plots of T (where T is the temperature and « 

(tj-tj^) [ln(t2 /tj) l.e., the Inverse of the rate window^®^) are sho«m in 
Fig. 3. It is seen that the activation energy of the emission rate of bulk 
levels is the same for MOS and MS structures, and equal to 0.81 eV. The 
activation energy of the interface states is found to be 0.65 eV. The corres- 
ponding capture cross-sections determined from transient capacitance spectra 

*i3 2 

(assuming that they are independent of temperature) are about 10 cm and 

lO”^^ to 10**^^cm^ for bulk and Interface states, respectively. The present 

value of the capture cross-section of the interface states is far more realistic 

-8 -9 2 

than the extremely large values of 10 to 10 cm reported on the basis of 
conductance measurements. 


- 6 - 


The trenllenc cepacicancc spectra as a function of the filling pulse 
voltage in Fig. 2 show that the high temperature peak associated with 0.81 eV 
bulk states saturates at about 20 V. The lower temperature peak, associated 
with interfac<. states, is negligible at low filling pulse voltages, but it 
increases sharply for pulses exceeding 20 volts. This markedly different 
behavior of the bulk traps and Interface states permits the quantitative 
determination of their contributions to the total DLTS spectra. For a single 
energy level, tb^ t^’ ; eralurc dependence of the DLTS olgnal, ticasureci nt 
r.et..plin;\ t^'r^a t^ '• C(t 2 ) - ir. de; ibzd by equ'^tioi-.h,: 

AC(T) - exp (-tj^/i) (1 - exp (-t^/t)} 



*Vth"c 


( 2 > 


capture 

where tC is the DLTS signal at the peak maximum, is the electron/ cross- 
section, v^^ is the thermal velocity of electrons, is the density of states 
in the CaAs conduction band, !■ the activation energy of states, and k 

is the Boltzmann's constant. 

As shown in Fig. 4, aC(T) calculated fron eq. (2) assuming “ 0.807 cV, 

-14 2 

0 ■ 8.94 X 10 cm , and typical temperature dependence of v . and N is in 

n tn c 

excellent agreement with the experimentally determined transient capacitance 

spectrum obtained with a20V/20 p^ec pulse, associated with the bulk states. 

To separate the DLTS pes!; associated with the interface states, this saturation 

value of the bulk peak was subtracted from the total DLTS signal at higher 

filling pulse voltages. As shown in Fig. 5, the resulting interface peak is 

in good agreement with theoretical calculation assuming a discrete energy 

level at • 0.651 eV, characterized by an electron capture cross-section 

0 - 7.0 X lO'^^cm^. 

n 
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structures 


By extending earlier treatments of the capacitance transients in MOS 
(15) 


to DLTS measurements » it is readily shown that the concentration 


'to 


of the bulk states is given by 


K -6 AC ”d-”a 


while the density of interface states N is 

oS ' 


O) 






(A) 


v/i i *.rc C 5.n the tot:.! o2 t!i: t wi vtcJurc, C ^ id the oxide cepr:cJ.tancc, 

» on 

and is the net doping concentration. The concentration of the 0.81 eV 

• * 14 -3 

bulk states was found to be about 10 cm in all measured MOS and MS structures 

(in the case of MS structure in eq. 3 is equal to zero) and was Independent 

of depth. The density of the 0.65 eV interface states evaluated from the 

12 -2 

saturation value of the DLTS signal was of the order of 10 cm . 

The position of the discrete interface level reported in this paper is 


consistent with the energy-density distribution of the oxide-GaAs interface 

* and 
(16) 


states previously obtained by using Berglund's^^^ and Terman's^^^ methods as 


well as by the saturation photovoltage technique. As shown in Fig. 6, 

these results show an increase of the interface state density towards the 
energy corresponding to the position of the presently observed discrete inter- 
face level. It should be noted, however, that recent DLTS measurements of 
GaAs MOS structures prepared by oxidation in an oxygen plasma, induced by 

9 

high-frequency discharge, revealed interface states with continuous energy 

13 -2 -1 

distribution (about 10 cm eV ), exhibiting a peak at 0.43 eV below the con- 
duction band edge and an increase in density toward the valence band edge.^^^^ 

These results were obtained below room temperature. In the present study, transient 


C 


- 8 - 

capacitance measurements carried out under conditions similar to those 
reported In ref. 11 revealed no measurable DLTS signal below room temperature. 
Thus, It Is concluded that the Interface states at 0.43 eV are characteristic 
of MOS structuret. prepared by plasma oxidation rather than by anodic oxidation. 

IV. SUMMARY 

Interface states and bulk GaAs levels were simultaneously Investigated 

In GaAs MOS structures prepared by anodic oxidation in propylene and dlethylene 

jjlycol elcclroI-ytcE They \* 2 rc cuccccs:fully distinj^uished by carryirg cut a 

ct ; -.ierai i VP ;.r;; l.yyis of drop Icvi l tibiisfeiit c'.'prcitpnr.e f j-ecirc of and 

MS structures prepared on the same samples of epitaxially grown n-type GaAs. 

The transient capacitance spectra of both MOS and MS structures exhibited a 

peak corresponding to a bulk level at 0.81 eV below the conduction band and 

14 »3 electron 

with a concentration of about lO^^cm Independent of depth. The/capture 

-13 2 

cross-section of these states was found to be of the order of 10 cm , 
assuming to be Independent of temperature. The 0.81 eV bulk states can be 


attributed to a level commonly identified with oxygen,^ although there is 

(18) 


(17) 

An additional DLTS peak ob- 


no unambiguous proof of this identification, 
served only In the case of MOS structures was found to be very sensitive to the 
filling pulse voltage consistent with the behavior of Interface states. The 
results can be readily accounted for on the basis of Interface states with a 
discrete energy at 0.65 eV below the conduction band, with an electron capture 
cross-section of 10 to 10 ^cm and a density of about 10 cm . Thus, it 
is concluded that In the GaAs-anodlc oxide Interface near the mid-gap, there 
are Interface states with a discrete energy level rather than with a continuous 
energy distribution as previously reported. 
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FIGURE CAP’ilONS 


Fig. 1. Typical capacitance-voltage characteristics of GaAs NOS structures 
prepared by anodic oxidation in propylene glycol electrolyte for 
two sweep rates. 

Fig. 2. Dependence of transient capacitance spectra of NOS structures on 
filling pulse voltage; sampling times t^/t 2 ~ 5/10 msec. 

Flp. 3. Typical thermal activation plots for NOS and MS structures, 
rjj.. ^ 1 . Conpc risen of ncasuied (tj/l, 5/10 injec) uml calculated ti ancient 
cnp:.''.iti ncc spect i;r. o! bulk Pr.vctfc.iCT. a ur.ud iu the f i u i;ip, 

procedure: Er~E.. “ 0.807 eV and o_ ■ 0.9A x 10 *cm . 

V# t ^ la 

Fig. 5. Comparison of measured and calculated transient capacitance spectra. 

Total DLTS signal was obtained with sampling times **■ 5/10 msec 

and filling pulses 40 V/0.5 msec. The experimentally determined 

(see text) contribution of the bulk and interface traps Is shown 

by dashed lines. Parameters of interface states used in the fitting 

-15 2 

procedure are; E^-E^ ■ 0.651 eV and ■ 7.0 x 10 cm . 

Fig. 6. The interface state distribution in the energy gap as measured by 

Terman's method (after Sawada and Hasegawa - ref. 6) and saturation 
surface photovoltage technique (after Shlmano, Moritanl and Nakai - 
ref. 16). 
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AB8TEACT 




Thaoratleal and axparlaantal studlat of tha alaetron nobllltj and tha 
fraa^earrlar abaorptlon of n-typa ZnP vara earrlad out In th« tanparatura 
ranga of 77-30C*K. All najor aeattarlng procaaaaa and aeraanlng af facta vara 
takan Into eonaldaratlon. Zt vaa found that tha axparlaantal dapandanca of 
alaetron aobUltT and fraa-earrlar abaorptlon on tanparatura and/or on carrlar 
eoneantratlon can ba eonalatantly axplalnad only whan tha affact of conpanaa- 
tlon la quantltativaly takati Into account. Convanlant procaduraa 

ara praaantad for tha dataxniantion of 

tha eonpanaatlon ratio fron tha valuaa of alaetron nobility and fron tha fraa- 
carrlar abaorptlon coafflelant. Tha high contribution of optical phonon acat- 
taring In ZnP llnlta tha applicability of tha fraa-carrlar abaorptlon approach 
to alaetron concentration n > lO^^cn'^. Electron nobility, bowavar, can ba 
rallably aaployad for tha datamlnatlon of tha eonpanaatlon ratio for 
n > lO^^en"^ at 300*K and n > 10^*cn'^ at 77*K. 


Journal of Appllad Physlca, in praaa. 


I. mnoDucTioii 


TIm eoapMMCloD ratio la a vary laportant paraaatar la cha charactarlta- 

clOB of aaaleooductlag aatarlala. la raeaat yaara a au^ar of proeaduraa for 

avaluatlag thla paraaatar la a-typa aaBleoadueeora hava baaa propoaad. Soaa 

of thaaa*Mthodo raqulra axtaaalva lov taaparatura ■aaauravaata aad eharafora 

(1 2 ) 

ara of llaltad practical laportaaca. ' * ' Tha tlaplaat aad Boat vldaly used 
aathoda ara bated on tha fladlag that klaatlc phaaoBaaa such aa alactroa 
Boblllty aad/or fraa alactroa abaorptloa dapaad on tha total coacantratlon of 
loalaad lapurltlaa aad* thua» oa tha coapaaaatloa ratio. For Ill-V eonpound 
aaBleooductors» tha ralavaat thaoratlcal analysis la conpllcatad, slaca a 
ralaxatlOB tlaa cannot ba Introducad axcapt at vary low taaparaturas. Accor- 
dingly, oaa has to raly on auasrlcal calculatloaa^ lacludlag tha conpanaaclon 
ratio, aa a paraaatar. 

Thla type of thaoratlcal analysis of alactroa nob 111 ty and fraa-carrlar 
absorption In a-typa CaAs has baan carried out recently In coajuaetlon with 
aa axparlBsatal study. It has lad to tha foraulatloa of practical procedures 
for tha dataninatloB of tha coapaasatloa ratio over a wide range of fraa alac- 
troa coBcantratlona. la tha prasaat study this approach la axtandad to InP, 
which racaatly has baaa racalvlBg a groat daal of attention as a potential 
aubstltuta for CaAa la a auabar of nlcrowave aad optoalactronlc applications. 

Tha alactroa BOblUty la laP has baaa axtanalvaly studied, both thao- 
ratlcally aad axparlaaatally. Thus, It hat baaa soowa that theoretical 
aobUlty Halts agree with axparlaeatal data obtained with good quality un- 
coapaaaatad Bstarlal. Only a Halted effort has baan directed toward tha 
study of fraa-carrlar absorption in laP. There ara just two partlaant axparl- 
aaatal reports. Fkirtharaora, tha astarlals paraaatars used in the 


thtor«tic«l analysis of fraa-carriar abaorption' * hava alnca baan ravlsad. 

In tha present calculations of the electron mobility and the fran-carrlar 

absorption a sat of up-dated materials parameters la employed. The mobility 

(23) 

calculations are based on a variational procedure while the free-carrler 
absorptltm la calculated utilising the expression of ref. 18, derived with a 
perturbation theory. 

II. THEORETICAL CONSIDERATIONS 

Electron mobility and free-carrler absorption depend on the band structure 
parameters and the scattering processes which determine the ^«mentum relaxation 
of the free electrons. In III-V compound semiconductors the following scattering 
mechanisms play the most la^ortant role: electron-optical phonon, electron- 

acoustical phonon (Including plexoelectrlc scattering) and electron-ionized 
Impurity scattering. For InP the conduction band parameters and the electron- 
phonon coupling constants are fairly well known. Thus, InP Is a direct and 
wide gap semiconductor (Eg - 1.35 eV at 300*K) with the conduction band minimum 
located at the F point of the Brlllouln zone. The conduction band dispersion 
relation Is well described within the spherical and parabolic approximation 
with a constant effective mass m*. This approxlmxitlon allovs direct application 
of the relationships for the mobility and free-carrler absorption In the form 
given In ref. 4. 

All major scattering mechanisms, namely polar optical, deformation poten- 
tial acoustical and Ionized Impurity, are Included In these calculations. The 

piezoelectric acoustical phonon scattering Is neglected since In InP the 

2 

piezoelectric coupling parameter p Is -16 times smaller than that of GaAs 
In which the piezoelectric scattering contribution to the room teoperature 
mobility does not exceed 2Z. The deformation potential-optical phonon scat- 
tering (not considered In ref. 4) will also be neglected here, since in most 
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coapound Malconductora this scattarlng process Is efficiently nasksd by the 

h ‘ i _ * 

much stronger polar Interaction. 

The values of InP paraneters uaed In the computations are given in 
Table 1. It was assumed that all parameters other than tha effective t?ass 
are temperature Independent In the temperature range considered. 

A. Electron Mobility 

Due to the large optical phonon energy the electron-optical phonon scat- 
taring la an Inelastic process and thus the relaxation tlma approximation 
cannot be used In the mobility calculations. Accordingly, the varia- 
tional procedure in the form given In ref. 23 with modifications Introduced for 
free -‘electron screening effects Is used. This procedure docs not Invoke the 
Matthlessen rule (according to which average macroscopic mobilities resulting 

from different scattering mechanisms are summarized, 1/p ■ l^l/y^) since In 

(4 19) 

III-V compounds It leads to unacceptable errors. * Thus, In order to find 
the electron mobility as a function of free carrier concentration and compensa- 
tion ratio It Is necessary to perform numerical computations. 

Figure 1 shows the calculated temperature dependence of the total electron 
mobility In InP with two different free electron concentrations. In addition, 
the temperature dependence of the two most Important component mobilities is 
plotted. The Ionized iiq>urity contribution to the total mobility rapidly In- 
creases with Increasing* tepmerature. For the material with the lower dopant 
concentration the room temperature mobility Is practically Insensitive to 
Ionised Impurity scattering and its value Is determined mainly by optical 
phonon acatterlng. In such a case, more reliable results for the compensation 


I 


ratio can be obtained from the mobility value at the liquid nitrogen temperature. 

The calculated room and liquid nitrogen temperature mobilities as a function of 
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fr«*-carrl«r concentrations end compensation ratio are presented In Tables II 
and III* . v\L-- 

The results for the lattice mobility limits obtained In the present study* 

5 2 2 

l.e.* -3 X 10 cm /Vs at liquid nitrogen and ~4800 cm /Vs at room teoq>erature» 

are In gb^ agreement with thoae reported in ref. 19* l.e. * *>3.2 z 10^ and 

2 

>4700 cm /Va* respectively* where an Iterative procedure was used to obtain 
numerical solution of the Boltamaxin equation and the same set of InP parameters 
as In the present calculations. 

B. Free Electron Absorption 

In contrast to the mobility* the free-carrler absorption may be represented 
as a sum of absorption coefficients associated with various scattering processes 
cons! 'sred In the present paper; the total absorption coefficient Is: 


" ®«n + ®-« + 

1 op AC 


where a * a and a. are the absorption coefficients associated with electron- 
Op AC inp 

optical phonon* electron-acoustical phonon and electron-ionized Impurity Inter- 
action* respectively. 

Numerical calculations of absorption coefficients were carried out 
employing appropriate expressions for the component absorption coefficients 
as given In refs. 1 8 or 4. The results for room and liquid nitrogen temperatures 
at a wavelength of 10 ym are given In Tables II and III, respectively. 

As in the case of GaAs*^ ' It can be shown that the component absorption 
coefficients at any wavelength, X* can be obtained from the following approxi- 
mate relationships: 


“op*X " ®op*X^ " 

^’fimpfX " ®4mp,X 


®ac*X " ®ac,X^ (X/X^)^*^ 
o o 


(2c) 
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Thla approxlBatlon Is vslld for a-type InP with fres-carrlsr concsntrstlons 
18 ~3 

up to about 10 cm and for wavalengths X < 10 ym. Tha above expreaslons 
ara valid at larger wavalengths for lower dopant concentrations. 

The results of Tables II and III can be used to determine the absorption 
coefficient In saiq>les with a given Ionized Impurity concentration, namely 

®exp " n ®lmp ®ap * ®ac 

where la the total concentration of Ionized impurities and n la the free 

electron concentration. Using the above relation one can express the compensa- 
tion ratio In terms of the experimental and theoretical component absorption 
coefficients 

9 . ~ * V Vp^ 1 ^. 

V * %>■ * “•=) 

where 6 ■ ratio of the ionized acceptor to the Ionized donor 

concentrations. 

III. EXPERIMENTAL 

Hall-effect measurements, conductivity and Infrared absorption have been 
carried out on n-type InP samples cut from melt-grown single crystals with 
electron concentrations In the range of 3 x 10^^ to 2.5 x lO^^cm Assess- 
ment of the saiq>le homogeneity was found necessary in order to obtain consistent 

and meaningful results. The measurements carried out for this assessment, 

(28 29) 

nasiely the carrier concentration profiling. with IR laser scanning^ * and 
the Ball-constant measurements as a function of magnetic field up to 150 KOe, 
are presented In Appendix I. In considerable portions of melt-grown crystals 
the spatial variation of the absorption coefficient, at X ■ 10.6 ym was found 
to exceed 30Z of the average value. In such Inhomogeneous samples the Hall 


constant exhibited an anomalous magnetic field dependence consistent with 
the Hall-effect analysis of Ref. 25. A reliable interpretation of electron 
mobility and free carrier absorption on the basis of average macroscopic values 
is not possible for inhomogeneous material. Accordingly » only selected samples 

s 

s 

with carrier concentration Inhomogeneitles not exceeding lOZ and exhibiting 
the standard Hall constant dependence on magnetic field are considered in this 
study. For these samples, the electron siobillty and electron concentration 
were obtained from conductivity and Hall-effect measurements. The saturation 
value of the Hall constant, in high magnetic field, at which the Hall factor, 
r ■ 1, was used to determine the electron concentration. In this way common 
ambiguities associated with the value of the Hall factor, r, were eliminated. 
The measured parameters of the samples employed are listed in Table IV. 

A. Electron Mobility 

In the present theoretical treatment, the electron drift mobility was 
taken as p ■ o/en, where o is the low DC electric field conductivity. Thus, 
a comparison between theory and experiment requires experimental data of 
drift mobility rather than of the commonly measured Hall-effect mobility. 

Mg “ ru. Since the present measurements of the Hall constant, Rg, were deter- 
mined in a high magnetic field, where r ■ 1, the electron concentration 
(n ■ r/eRg) and the drift mobility can be reliably determined. The results 
on the room temperature mobility a« a function of electron concentration are 
presented in Fig. 2, together with the theoretically calculated values for 
different compensation ratios, 6. It is seen that the presently determined 
mobility values are below the theoretical limit for uncompensated InP (6 ■ 0). 
Utilizing Table III, the values of compensation ratios were determined for 
these samples and are given in Table IV. 
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Th« llttraturt nobility plottad In Fig. 2 rapraaant Hall 

mobllltlaa» and thua It la not aurprlalng that In tha raglon of low alactron 
concantratlona* aona valuaa azcaad allghtly tha llnlt of tha drift mobility. 
According to theoretical conalderatlona preaanted in ref. 26, tha Hall factor, 
r. In pura InP at room temperature la about 1.25. With Increasing electron 
concentration the value of r decreases, approaching r ■ 1 for degenerate 
material. Thus, for high electron concentrations, the Hall mobility approaches 
the drift mobility value, and the literature results fall below the theoretical 
curve corresponding to uncompensated material. 

Mobility data at 77*K are given In Fig. 3. At this temperature the dif- 
ference between the Hall-effect and drift mobility does not exceed 
l.e.. It is within the experimental accuracy of Hall-effect measurements. Thus, 
the present results and the literature data fall below the theoretical curve 
of 9 "0 over the entire electron concentration range. The compensation ratios 
determined from Table III are given In Table IV. A comparison with the results 
obtained at 300*K shows that there Is a noticeable Increase in the compensation 
ratio at 77*K for the samples In which the electron concentration decreases at 
77*R (samples 4, 5 and 7). For the other samples the differences in compensation 
ratio do not exceed the experimental error. However, It should be noted that 
for electron concentrations below 10^^, the mobility at 77 *K is much more sen- 
sitive to ionized liq>urities than at 300*K. 

The necessity to Include the compensation ratio In the theoretical analysis 
of the mobility Is especially evident In the dependence of mobility on temperature. 
Thus, In Fig. 4 the experimentally determined mobility Is shown as a function of 
temperature for two InP samples with different free electron concentrations. The 
theoretical curves for uncompensated samples lie well above the experimental 
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value* and exhibit a distinctly different temperature dependence. On the 
other hand, as shotm In the same figure, there Is a vary reasonable agree- 
ment between theory and experiment If the compensation ratio Is taken Into 
account In the theoretical calculations. 

It Is of Interest to note that, according to the present analysis, melt- 
grown InP exhibits much lower compensation than that encountered In melt- 
grown GaAs. This difference Is evident In Fig. 5 where the normal range of 
Nimp/n Is plotted as a function of electron concentration for commercially 
available high-quality GaAs' ' and for InP. In the low electron-concentration 
region the melt-grown GaAs Is highly compensated, with a total concentration 
of Ionized Impurities exceeding by an order of magnitude the free carrier con- 
centration. InP, on the other hand, even at low concentrations can be obtained 
by melt groirth with a concentration of Ionized impurities of the same order of 
magnitude as the electron concentration. This marked difference Is probably 
related, In part, to the lower melting point of InP (1062*C for InP and 
1238*C for GaAs) and thus to the lower concentration of native point defects; 
point defects In conq>ound semiconductors are known to participate In the 
formation of compensation centers. Furthermore, the covalent radii of In and 
P differ by about 40Z, whereas those of Ga and As are about the same; thus, 
a dopant Impurity Is more likely to occupy preferentially In or P sites In 
InP than Ga or As sites In GaAs (amphoteric or compensating doping). 

B. Free Carrier Absorption 

Room temperature IR absorption was determined from transmittance 
measurements employing a Fourier spectrometer. The absorption coefficient, 
a. Is obtained from the relationship: 
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X • (1-R) expC-gd) 
1-R^ exp(-2cd) 


( 5 ) 


where T le the transmit tence, d Is the thickness of the ssiq>le, and R Is 
the reflectance. Values of R as a function of wavelength* X* for InP are 
given in Ref. 27. 

The present experimental results on the absorption coefficient at 
X ■ 10 ym together with those in Ref. 16 are presented In Fig. 6 as a function 
of electron concentration. The theoretical values of the total absorption 
coefficient and the absorption components associated with different scattering 
mechanisms are also shown In this figure. It Is seen that the major contri- 
bution to free carrier absorption In InP Is provided by optical phonon scat- 
tering. The Ionized Is^urlty scattering Is of significance only at high con- 
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centratlons* n > 10 cm . 

As pointed out In Ref. 4* a comparison between the experimental and 

theoretical dependence of free carrier absorption on the electron concentration 

(for samples having different compensation ratios) can be achieved by reducing 

the experimental absorption coefficient* a * to an absorption coefficient 

* ^ N. 

corresponding to zero compensation, o = o - (-^ - Da. ; which according 
to ref. .(4) Is; h exp n imp 


o* 

exp 




( 6 ) 


given In Table II. In Fig. 7, a* Is plotted as a function of 


The values of g*^ were calculated using the values of compensation ratio deter- 
mined for each sample from Its electron mobility (Table IV) and from the computed 
values of 

carrier concentration together with the theoretical dependence oi^(n). It Is 
seen that good agreement la obtained between theoretical and experimental results. 

Using Eq. (4) and the theoretical values of the absorption coefficient 
given In Table II* values of the compensation ration for the samples Investigated 
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v«r« obtained. The reaulta are aumaariaad In Table IV. It la aaan that 
the valuaa obtained by the optical oathod are in aatiafaetbry agraamant irith 
thoaa datarainad from electron nobilitiaa. 

The theoretically calculated valuaa of the abaorption coefficient at 
77 *K ara^ ahoun in Fig. 8 aa a function of carrier concentration. A coapariaon 
with the room temperature reaulta in Fig. 6 ahowa that the optical phonon 
contribution to the total absorption coefficient does not dacraasa signifi- 
cantly with decreasing teaparatura. Thus, lov *’eoT>erature IR absorption is 
not appreciably aore sensitive to the ionised impurity scattering than room 
taaperatura absorption. This behavior la significantly different from that 
of the electron mobility discussed in the previous section. The optical phonon 
contribution to the total aobillty dacraases rapidly with decreasing tamparattire 
and thus the low temperature mobility is very sensitive to the ionized impurity 
scattering. This different behavior of two closely related processes can be 
explained as follows: At low temperatures when the optical phonon occupation 

number is very small the main contribution to electron optical phonon scattering 
consists of transitions in which optical phonons are emitted by electrons. 

In low DC fields phonon emission requires a transition of sn electron from 
an occupied state at a given energy E to an empty state at an energy lower by 
the phonon energy, Such a process is hardly possible at low temperatures 

(kT«lkJ ) whereby the density of final unoccupied states is essentially negll- 
o 

gib la. Free carrier absorption Involves additional intaractlon with an incident 
photon. Thus, even at low tamperaturas the photon esiisslon is energetically 
posaibla providad that tha anergy of the absorbad photon, hv, is higher than 
the energy of the emitted phonon 


IV. SUMMAKT 


Th« •Icctron aobillty «nd IR fr«« alcctron abaorptlon in n-typ« InP hnv« 
bacn InvaatlgaCad theortcicallj and axparimantally. Thaoratlcal calculationa 
of mobility basad on a variational procadura agraa vary wall with pravloualy 
publiahad reaulta. The room tamperatura mobility In pura In? la domlnatad by 
optical phonon acattarlng, a contribution of laaa than lOX mada by acouatlc 
phonon acattaring. On tha othar hand, tha room taoparatura abaorption coafflcl- 
ant of pura InP la antlraly datarmlnad by alactron-optlcal phonon Intaractlona; 
tha contribution from tha othar acattaring procaaaaa la nagllglbla. With In- 
craasing doping tha lonlzad liiq>urlty contribution to tha total acattaring 
Incraasaa, which parmlta tha lonlzad l^urlty concantratlon to ba datarmlnad 
from tha mobility and fraa carrlar abaorption maaauraments. At lowar tampara- 
tures the optical phonon acattaring contribution to the mobility dacraaaas; 
thus, the Interaction between alactrona and lonlzad Impurltlaa bacomea tha 
main process controlling the electron mobility over a wide range of free 
electron concentrations . 

In contrast, in fraa carrlar absorption, the absorption coefficient re- 
sulting from the optical phonon scattering decreases slowly with decreasing 
temperature, and Its value at liquid nitrogen temperature is smaller onl'* by 
a factor of two than Its room temperature value. A comparison of theory with 
experimental data showed that the experimental values of the room temperature 
mobility and of the free carrier absorption coefficient, as well as the 
temperature dependence of the mobility can be consistently explained if the 
compensation ratio of the Investigated material Is taken Into account. 

Tha tabulated values of electron mobilities and free carrier abaorption 
coefficients provide a practical means for the determination of the compenaation 
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ratio In n-typa InP, with n > 10 cm from room tonparatura maaouramanta 
of althar tha nobllltj or tha fraa carrier absorption coafflclant. For fraa 
carrier concantratlona ranging fron 10^^ to lO^^cm ^ tha companaatlon ratio 
la datamlned from nobility neaauraaanta at tha liquid nitrogen tonparatura. 

Aa shown recently » under certain condltlona the nobility calculatlona 
for n-type GaAs can be utilised to eatimate nlnorlty carrier nobility In p~type 
GaAa.^^^^ This nethod can be directly applied to InP. 

Finally, It waa found that data on mobility in InP are vary cloae to 
theoretical values obtained for unconpensated material; in contraat, for 
GaAa with aimilar free carrier concentratlona the free carrier nobilltlea 
Indicate compenaatlon ratioa In the range of 0.6 to 0.8. 
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Appmdlx X 

ASSESSMENT OP HOMOGENEITT OF InP 

TIm «l«ccron conctntratlon profllas of InP oamplot voro sMourod with 

(28 29) 

tho IE Xdsor acannlng tochnlquo. * * Tho oxporliMntnl arrangtaont conalatad 

of a CO 2 laaar (tunabla from 9.2-11.8 ym), a ataga vlth aoclont an optical 
ayataa for aonltorlng cha powar and vavalangth of tha laaar baaa, and an IR 
datactor. Tha optically flat vafar vaa nountad on tha ataga and tha Intan- 
city of a colllaatad baas (-30 ya In diaaater) paoalng through tha vafar vaa 

aaaaurad as a function of tha poalt‘^m along tha aampla. Tha II tranaalaalon 

an 

acana vara convartad Into /abaorpt Ion coafflclant utilising Eq. (3). Tha aaasura- 
aants vara rapaatad for dlffarant vavalangtha of tha laaar radiation In ordar 
to aaparata posslbla Intarfaranca affacta rasultlng froa alight variations 
In vafar thicknass. 

Tha alactron concantratlon profllas vara datamlnad from abaorptlon coaf- 
flclant profllaa utilizing tha tabulatad thaoratlcal valuaa of tha fraa car- 
rlar absorption coafflclant. Typical raaulta obtained vlth nalt-grovn InP 
ara glvan In Pig. 9. It la aaan that local changaa In alactron concantratlon 
ara as high as -391 from tha averaga value. Similar large fluctuations of 
carrier concantratlon vara found In about 201 of the nalt-grovn InP crystals 
obtained from various sources. Reliable analysis of such Inhomoganaous 
material on tha basis of tha average optical and electrical properties Is not 
posslbla. 

As shovn In Pig. 10, Inhomoganaous InP samples exhibit anomalous Hall 
constant dapandanca on the magnetic field; In samples 1 and 2 tha Hall 
constant Incraasas In high magnetic fields and In sample 3 It even passes 
through a minimum. In n-typa homogeneous samples tha Hall factor r (and 
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I 

thus ths Hall constant K ■ r/sn) should dscrssss nonotonlcslly to ths vslus 

of r ■ 1 (or K ■ 1/sn) with Incrssslng oagnstle field. 

This anomalous Ball constant behavior is in excellent agreement vlth 

( 25 ) 

a recent treatment of the Ball affact in Inhomogeneous semiconductors 

a 

where it is shown that the standard analysis of DC conductivity and low nag- 
netic fiald Ball constant applied to Inhomogeneous material can lead to false 
mobility values noticeably exceeding theoretical llalta. It should also ba 
noted that the standard analysis of the average free carrier ebsorption and 
low field Ball constant leads for inhomogeneous samples to an III absorption 
coefficient which is appreciably smaller than the theoretical value. 

In the present study cf the electron mobility and the free carrier absorp- 
tion only selected samples were employed which exhibited standard Hall constant 
dependence on the magnetic field » and fluctuations of electron concentration 
not exceeding -lOZ. 
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Table I* InP parametera used In the 


present computations 


' Low-freqency dielectric 

i ' constant 

I 0 

; High-frequency dielectric 

i constant e_ 

• Optical-phonon 

energy hu^ 

S Deformation 

I potential, 

Longitudinal 
elastic constant, 

Effective mass, m*/m 

o 


(after Ref. 7 and 10) 

12.38 ' 

,.55 

42.8 meV 
6.8 eV 

12.1 X 10^^ dyne/cm^ | 

0.082 (77*K) 0.078 (300*K) j 

] 

j 



1.0x10^^ 

4710 

4690 

4660 

4620 

4570 

4510 

4410 

4270 

4160 

4020 

3820 

3490 

2900 

1.5 

4670 

4630 

4590 

4540 

4480 

4390 

4270 

4090 

3970 

3800 

3560 

3210 

2610 

2 

4640 

4600 

4550 

4490 

4410 

4300 

4160 

3960 

3810 

3630 

3380 

3010 

2420 

3 

4570 

4510 

•4450 

4370 

4270 

4140 

3970 

3730 

3570 

3360 

3090 

2720 

2150 

4 

4520 

4450 

4370 

4280 

4160 

4010 

3820 

3560 

3380 

3170 

2900 

2530 

1980 

5 

4460 

4390 

4300 

4190 

4060 

3900 

3690 

3410 

3230 

3020 

2750 

2380 

1840 

6 

4420 

4340 

4240 

4120 

3980 

3810 

3590 

3300 

3120 

2900 

2620 

2270 

1740 

7 

4380 

4290 

4180 

4060 

3910 

3730 

3500 

3200 

3010 

2790 

2520 

2170 

1650 

8 

4330 

4230 

4120 

3990 

3830 

3640 

3410 

3110 

2920 

2700 

2440 

2090 

1570 

9 

4300 

4200 

4080 

3940 

3780 

3580 

3340 

3030 

2850 

2630 

2360 

2020 

1510 

1.0x10^^ 

4250 

4140 

4010 

3870 

3700 

' 3500 

3250 

2950 

2760 

2550 

2290 

1950 

1450 

1.5 

4080 

3950 

3800 

3640 

3460 

3240 

2990 

2690 

2500 

2300 

2044 

1720 

1240 

2 

3960 

3820 

3660 

3480 

3290 

3070 

2810 

2510 

2330 

2120 

1880 

1560 

1100 

3 

3770 

3610 

3440 

3250 

3050 

2820 

2570 

2270 

2100 

1900 

1660 

1360 

918 

4 

3630 

3460 

3280 

3090 

2880 

2660 

2400 

2110 

1940 

1750 

1510 

1220 

800 

5 

3530 

3350 

3170 

2970 

2770 

2540 

2290 

2000 

1830 

1640 

1410 

1120 

718 

6 

3440 

3270 

3080 

2880 

2670 

2440 

2190 

1900 

1740 

1550 

1320 

1040 

654 

7 

3370 

3190 

3000 

2800 

2590 

2360 

2120 

1830 

1660 

1480 

1250 

973 

604 

8 

3310 

3130 

2940 

2730 

2530 

2300 

2050 

1770 

1600 

1410 

1190 

920 

563 

9 

3260 

3070 

2880 

2680 

2470 

2250 

2000 

1710 

1550 

1360 

1140 

875 

530 

1.0x10^^ 

3220 

3030 

2830 

2630 

2420 

2200 

1950 

1660 

1500 

1320 

1100 

835 

501 

1.5 

3090 

2900 

2700 

2490 

2280 

2060 

1800 

1520 

1360 

1180 

970 

723 

423 

2 

2960 

2760 

2560 

2350 

2140 

1910 

1660 

1370 

1210 

1040 

839 

611 

345 

3 

2840 

2640 

2430 

2210 

1990 

1750 

1500 

1220 

1060 

896 

712 

507 

277 

4 

2770 

2560 

2340 

2120 

1890 

1650 

1400 

1120 

967 

806 

633 

444 

239 

5 

2720 

2500 

2270 

2050 

1820 

1570 

1320 

1040 

896 

741 

577 

400 

213 

6 

268C 

2450 

2220 

1990 

1760 

1510 

1260 

987 

844 

694 

536 

369 

194 

7 

2650 

2410 

2180 

1940 

1710 

1460 

1210 

942 

802 

657 

505 

346 

181 

8 

2620 

2380 

2140 

1900 

1660 

1420 

1160 

900 

764 

623 

477 

325 

169 

9 

2590 

2350 

2110 

1860 

1630 

1380 

1130 

867 

733 

596 

455 

309 

160 

1.0x10^® 

2570 

2320 

2080 

1830 

1590 

135P 

1100 

839 

708 

573 

436 

295 

152 

1.5 

2470 

2210 

1960 

1710 

1460 

1220 

979 

737 

617 

495 

373 

250 

128 

2 

2390 

2130 

1870 

1620 

1380 

1140 

904 

674 

561 

448 

336 

224 

113 

3 

2290 

2010 

1750 

1500 

1260 

1030 

811 

598 

495 

393 

293 

194 

98 

4 

2210 

1930 

1670 

1420 

1190 

965 

755 

553 

456 

361 

268 

177 

89 

5 

2150 

1880 

1610 

1360 

1140 

920 

717 

523 

431 

340 

252 

166 

83 

6 

2100 

1830 

1560 

1320 

HOC 

886 

687 

501 

411 

325 

241 

156 

79 

7 

2070 

1790 

1530 

1290 

1070 

861 

667 

484 

398 

314 

232 

152 

76 

8 

2040 

1760 

1500 

1260 

1040 

838 

647 

470 

386 

304 

225 

148 

74 

9 

2010 

1740 

1480 

1240 

1020 

821 

635 

462 

378 

296 

219 

144 

72 


19 


1.0 xlO 


1990 1710 1460 1230 1010 809 624 452 370 289 215 141 


71 
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Table Ilb. Free carrier absorption 


Concentration 

cn”J 


300*K, 
k • 0.078 

CM 

00 

• o 
1 

0 


“ac 

‘V»p 

®lnp 

‘Vc 

®OP 

1.0x10^® 

0.004 

0.034 

0.623 

0.003 

0.005 

0.320 

1.5 

0.008 

0.052 

0.932 

0.007 

0.008 

0.476 

2 

0.014 

0.069 

1.239 

0.012 

0.011 

0.629 

3 

0.031 

0.104 

1.850 

0.028 

0.016 

0.928 

4 

0.056 

0.139 

2.456 

0.048 

0.022 

1.217 

5 

0.086 

0.173 

3.051 

0.075 

0.028 

1.503 

6 

0.123 

0.208 

3.646 

0.106 

0.033 

1.780 

7 

0.167 

0.243 

4.240 

0.143 

0.039 

2.053 

8 

0.217 

0.278 

4.815 

0.185 

0.045 

2.321 

9 

0.273 

0.313 

5.397 

0.232 

0.050 

2.588 

1.0x10^^ 

0.314 

0.325 

5.578 

0.284 

0.056 

2.849 

1.5 

0.690 

0.491 

8.227 

0.617 

0.086 

4.127 

2 

1.201 

0.660 

10.79 

1.066 

0.117 

5.366 

3 

2.602 

1.005 

15.75 

2.296 

0.181 

7.793 

4 

4.474 

i;360 

20.52 

3.949 

0.250 

10.19 

5 

6.790 

1.726 

25.16 

6.009 

0.322 

12.60 

6 

9.510 

2.100 

29.65 

8.460 

0.398 

15.02 

7 

12.64 

2.488 

34. U 

11.27 

0.476 

17.41 

8 

16.13 

2.879 

38.44 

14.49 

0.560 

19.90 

9 

20.00 

3.285 

42.75 

18.05 

0.646 

22.38 

1.0x10^® 

24.22 

3.699 

47.01 

21.96 

0.736 

24.89 

1.5 

50.28 

5.912 

67.80 

46.69 

1.229 

37.69 

2 

83.91 

8.354 

88.02 

78.85 

1.798 

50.91 

3 

170.3 

13.87 

127.0 

165.13 

3.147 

78.60 

4 

276.7 

20.12 

164.1 

276.8 

4.751 

98.46 

5 

396.1 

26.98 

199.6 

405.0 

6.57 

118.6 

6 

522.8 

34.34 

233.4 

533.3 

8.527 

137.7 


Table III. Coniputed values of electron nobility 
In n-type InP at 77*K 
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Table IV. CoapensaClon ratloa detenilned froa nobility 
and IR absorption neasurenents 
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Flgur* 1 


Figure 2 


Figure 3 


Figure 4 


Figure 5 


Figure 6 


FinuRZ CAPTIONS 


Theoretical temperature dependence of electron mobility in 
uncompensated In? (see text). 

Theoretical (solid lines) and experimental values of mobility 
at 300*K as a function of electron concentration in n-type InP 
for various compensation ratios; + - present results (drift 
mobility): other experimental results correspond to measured 
Hall mobility; I - ref. 5; o - Ref. 10; Ref. 11; x - Ref. 12 
^ - Ref. 13; •- Ref. 14; A- Ref. 15. 

Theoretical (solid lines) and experimental values of mobility 
at 77 *R as a function of electron concentration In n-type InP 
for various compensation ratios; + - present results (drift 
mobility); other experimental points correspond to reported 
Hall mobility values: I - Ref. 5; ■ - Ref. 9, 

A - Ref. 11; ^ - Ref. 13; #- Ref. 14; Ref. 15. 

Experimental (+) and theoretical (lines) values of electron 
mobility in n-type InP as a function of temperature. Dashed 
lines correspond to uncompensated material. Solid lines repre- 
sent electron mobilities calculated using a compensation ratio 
3 Independent of temperature. 

Normal range of N^^p/n as a function of electron concentration 
in available melt-gro%m InP and GaAs (see text). 

Theoretical and experimental values of the room temperature 
(300*K) absorption coefficient as a function of electron concen- 
tration in InP (for X ■ 10 urn. # - present results; x - Ref. 16. 
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Figure 7. 


Figure 8. 


Figure 9. 


Figure 10. 


Theoretlcel (eolld line) end experlaentel valuee) obtained 
In preaent etudy) of the absorption coefficient reduced to 
sero conpensetlon (see text). 

Theoretlcel and experimental (Ref. 16) values of low temperature 
(77*K) free carrier absorption coefficient In n-type InP. 

Typical electron concentration mlcroproflle of InP obtained 
with scanning IR absorption. 

Hall constant dependence on magnetic field for Inhomogeneous 


n-type InP. 


Electron Mobility (cmVv-s) 










0*O.On 

\ 

\ 

\ 

\ InP 

\ + experiment 

\ ” theory 






Absorption coefficient (cm ) 



Absorption coefficient orexp (cm*^ 




Absorption coefficient (cm ) 




Rh(H)/R 



DETERMINATION OF CARRIER CONCENTRATION AND 
COMPENSATION MICROPROFILES IN GaAs 


L. Jascrzebskl» ^ J. Lagowski, W. Waluklavicz^ ^ and H. C. Gatoa 
Department of Materials Science and Engineering 
Massachusetts Institute of Technology 
Cambridge, Massachusetts 02139 

Simultaneous microprof lllng of semiconductor free carrier, donor 
and acceptor concentrations was achieved for the first time from the 
absolute value of the free carrier absorption coefficient and its wave- 
length dependence determined by IR absorption In a scanning mode. Em- 
ploying Ge- and Si-doped melt-grown GaAs, striking differences were 
found between the variations of electron concentration and those of 
Ionized impurity concentrations. These results showed clearly that 
the electronic characteristics of this material are controlled by 
amphoteric doping and deviations from stoichiometry rather than by 
impurity segregation. 


Journal of Applied Physics, In press 
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The slnultaneous determination of the donor and acceptor concentrationa 

la of particular Importance In III-V aemiconductor compounds which exhibit 

an appreciable degree of compensation. A nuotber of techniques utilizing 

spreading resistance, ^ ^ Schottky barriers, EBIC,*^*^ electroreflectance® ^ 

(6*9 ) dOA 

and IR absorption or topography have been developed for profiling 

the carrier concentration In semiconductors. However, these techniques 

essentially fall to distinguish between changes in impurity concentration 

and In compensation ratio. The present method Is an application of the 

02 ) 

recent quantitative studies of the free carrier absorption in GaAs and 
InP^^^ as a means of determining the compensation ratio. The principle of 
the method is understood by considering that the absolute value and the 
wavelength dependence of the free carrier absorption coefficient are experi- 
mentally measured quantities which contain only two unknown parameters, 
namely the electron concentration n • compensation ratio 

6 ■ concentrations of ionized donors and acceptors, 

respectively.) 

The experimental arrangement for microprofile determinations consists 
of a CO^ tunable laser (with about 90 lines in the spectral range of 9.16-11.02 ym) , 
a stage %rlth x-y scanning motion, an optical system for monitoring the power 
and wavelength of the laser beam, and an IR detector. Optically flat samples 
were mounted on the stage, and the intensity of a collimated bean (- 30 ym in 
diameter) transmitted through the wafer was measured as a function of its 
position along the aample. 

From the IR transmission scans the absorption coefficient, a, was obtained 
utilizing the atandard relationship: 



1-R‘«xp( - 2 od) 


1 where T le the trenemlttence, d Is the thickness of the sample, and R la 

j' 

the reflectance (values of R are tabulated in Ref. (14)). Measurements 
I I %wre performed with different wavelengths of the laser radiation in order 

I Ij to determine the wavelength dapandanca of the absorption coefficient. 

Common difficulties encountered in scanning IR absorption measurements 

f ( 

J } arise from Interferenc e affects. The large refractive index of semiconductors 

^ ’ in the Infrared region (usually for X « 10 ym, the refractive index exceeds 3) 

i 

enhances the interference due to alight variations of the wafer thickness. 

For thicker wafers, the temperature dependence of the refractive index can 
also lead to interference effects in the presence of temperature gradients. 

I These effects were essentially eliminated by careful polishing of the wafers, 

by utilizing a email power incident laser beam (to prevent heating of the 
wafer) and by performing the measurements under conditions of low transmittance 
(od > 2, eq. 1). The negligible role of Interference effects was verified by 
carrying out scanning IR absorption measurements at different wavelengths and 
for different values of incident beam power. 


In ref. (12) the study of free carrier absorption was primarily addressed 
to the determination of the compensation ratio, 6, from maasured values of the 


absorption coefficient, o . The simultaneous determination of the electron 


• d£ncx 

concentration p n ■ *^d“*^X* obtained If the logarithmic slope p = 

ie also available. 


For n-type GaAs at room temperature, a can be expressed as 


(15,12) 


a * 0 + 0 + a. N. / 
exp op ac imp imp/n 


( 2 ) 


where a and a are the absorption coefficients corresponding to electron 
op ac 


interaction with screened optical phonons and acouston phonons, respectively; 
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0 ^^ !• tht absorption coafflciant dua to the interaction of alactrona with 
acraaned ionitad impurities in uncompansatad material, and is the actual 

total concentration of ionised impurities + N^. For a given temperature 
and wavelength X the component absorption coefficients, and are 

determined* knowing the electron concentration, n. Values of these cofflecients 
for X ■ 10 pm are tabulated in ref. 12 as a function of n. 

The component absorption coefficients are characterised by different 


wavelength dependences. In the vicinity of X^ ■ 10 pm they can be expressed 


as 


( 12 ) 




X 


1 

a 

op 

■ a 

X 

* • 

^0 0 

a 

AC 

>* 

■ 

o 




1.5 


and 


°iap L " 


^ 0 


(3) 


By differentiating equations 2 and 3 it is shown that the slope 
p 5 ^ in the vicinity of X^ satisfies the relationship 


(3.5 - p) 


o +2o. 
op ac 

a 

exp 


(4) 


Computed value» cf a 4-2a as a function of electron concentration are 

op AC 

shown in’ Fig. 1. In the range of n considered, exhibits a slightly 
sublinaar dependence on n (screening effects) while Increases stronger 
than linearly with increasing n (effect of degeneracy). However, as seen 


in Fig. 1, a *f2a - 3.33x10 
op ac • 

convenient form 


,-17 


n. Thus eq. (4) can be rewritten in the 


n[cm“^] - 3.10^^(3.5-p)a _[cn"^l 

exp 


(5) 


which permits the direct determination of n from the measured values of p 


and a 
exp 


1 • 
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Aa shown In Fig. 2, tht valuta of p In n-typt GaAa can vary batvaan 

2.5 and 3.5, depending on Che electron concantracion and conpenaatlon ratio 

(or total concentration of lonlxad Inpurltiaa) . These changea of p with n 

(6*9) 

and/or 6 were neglected In previous creatDents of scanning IR absorption' ' 

n was 

assumed; such simplified relationship Is justified only when scattering by 


and topography ^^^*^^^where direct proportionality batvaan a and 

exp 


Ionized impurities is negligible. However, it Is evident from eq. 5 and Fig. 2 
that for n > 10^^ the reliable determination of the electron concentration in 
GaAs requires the measurement of the absorption coefficient and its wavelength 
dependence. When the electron concentration, n, is determined, the total con- 
centration of ionized impurities, can be obtained from expression 2, 

using the values of snd tabulated in ref. 12. Knowledge of 

n and permits the determination of the compensation ratio 6 and/or the 
concentration of ionized donors Np and acceptors N^. 

Mlcroproflles of free carrier, donor and acceptor concentrations obtained 
with melt-grown (horizontal Bridgman) Ge- and Si-doped GaAs are presented in 
figures 3 and 4, respectively. The Inhomogeneities seen In these figures, 
although of an amplitude larger than average, are typical of all crystals ob- 
tained from several sources. Actually, as a general rule, all presently 
svallable melt-grown GaAs is highly compensated and exhibits pronounced varia- 
tions in free carrier concentration. 

A striking feature of the resulta -of Fig. 3 is that the electron concen- 
tration exhibits large fluctuations of a magnitude comparable to its average 
value, while the total concentration of ionized impurities remains essentially 
constant. The average compensation ratio of this material was found to be 

about 0.4, which is a typical value for GaAs with an electron concentration 
18 —3 fl21 

in the range of 10 cm . ' ^ It is also seen in Fig. 3 that the spatial 


Ri 




J 


It is also seen In Fig. 3 that the spatial 


6 


variation* of donor and accaptor concantration ara out of phaao* i.a., 

maxina in coincida with ninima of 

Tha variation* of th* donor and accaptor variation* *hown in Fig. 3 

cannot b* axplainad on tha baai* of atandard impurity aagragation behavior 

controlldf by th* cryatal growth velocity or th* diffuaion boundary layer 

thicknea*. If it i* aaauaed that th* donor concentration variation* 

raflect atandard aegregation bahavior, th* naxlna in donor concentration auat 

correapond to naxina in growth rate or diffuaion boundary layer thickneaa, 

(18) 

ainc* th* diatribution coefficient of Ge la amaller than one. In that 

caae the acceptor Isipurity muat be an unintentionally added lapurlty with a 

diatribution coafflcient graater than ona, ainc* Ita concentration m a rla a 

coincide with th* Ge concentration ainlna. Howevar, there are no known lapuri- 

(18) 

tie* with a diatribution coafflcient greater than on* in the GaA* ayatea. ' 

Mor ia it likely that an unintentionally added. laq>urity 1* preaent at aa high 
a level aa 10 cn . Thua, both donor and acceptor level* auat be attributed 
to C* (amphoteric doning). Accordingly, th* aegregation of Ge (Nj^ + N^) under- 
goea only aaall fluctuetiona, and the free carrier variation* are due to the 
variation* of the relative Ge concentration in Ga (donor) aitea and in Aa 
(acceptor) aitea. Such varletiona (which are conaiatenc with fra* carrier con- 
centration fluctuation*) auat be aaaociated with factor* affecting th* relative 
concentratlona of Ga and Aa vacanciea or atoichloaatry of GaAa. Mo correlation 
between point defect formation and growth parameter* ia available. 

Similarly, th* fra* carrier concentration profile* in Si-doped GaAa ca nn ot 
be accounted for by atandard aegregation behavior. Thua, in Fig. 4a the donor 
and acceptor concentration change* arc not nearly a* pronounced. Th* result* 
of Fig. 4b were obtained from the same Si-doped crystal as those of Fig. 4a, 
but a different segment. In this case, the fluctuations of the free carrier. 
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th« donor and tht total lonlaad lapurlty 4 N^) concantratlona art of 
almilar natura and nagnituda, but tha lonlaad accaptor concantratlon undar- 
goal ralatlvaly anall variations. 

It la apparant fron flguraa 3 and 4 that fluctuations in tha Individual 
profllas of carrlar, donor* accaptor and total lonlaad lapurlty concantratlon 
can readily ba olstakan aa ralatad to standard Inpurlty aagragatlon behavior 
as they rasai^la wall estabUshad lapurlty aagragatlon variations, particularly 
In aleoantal semiconductors. Thus, all of thasa mlcroprolflas must ba sisul* 
tanaously considered in assessing the origin of Inhomoganaltlas in compound 
semiconductors . 

In suianary, a method is developed for tha simultaneous datarminatlon of 
free carrier, donor and acceptor microprof lias. It is alao shown that in 6a- 
and Si-dopad melt-grown CaAs dopant Inhomoganaltlas are not controlled by 
standard impurity segregation behavior, but rather by growth factors affecting ■ 
amphoteric doping and stoichiometry, which cannot ba as yat unambiguously 
defined. 
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FICUU CAPTIONS 

Plgurt 1. D«p«nd«nc« of a . ^ 2a . on oloctron concontration (••• tone). 

op oc 

Figuro 2. Calculntod logarithmic slope, p, of cha fraa carrier absorption 
^ cosfficiant of n-typa GaAs as a function of conpanaation ratio, 
6, la n*typa GaAs at room taaparatura. 

Figure 3. Electron concentration and ionised impurity nicroprofilea of 
Ge-doped melt-grown GaAs. 

Figure 4. El* ^tron concentration and ionised impurity nicroprofilea of 
Si-doped melt-grown GaAs. (a) and (b) correspond to different 
segments of the same crystal. 
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